0HHT1> 


4 


Q  ANGLE  AT  NOSEBOOH-HOUNTED  VANES 


KENNETH  RAWLINGS  III 


JUNE  19B1 


*  i  fC 

ELZCTE 
OCT  2  6  mC 


'*'*v 


THIS  DOCUMENT  HAS  BEEN  APPROVEO  FOR  PUBLIC  RELEASE  AND  RESALE 
ITS  DISTRIBUTION  IS  UNLIMITED 


AIR  PORCI  PL’OHT  TltT  CINTKR 
EDWARDS  AIR  PORCI  BASE,  CALIFORNIA 
—  AIR  FORCE  SYSTEMS  COMMAND  — 
UNITED  STATES  AIR  FORCE 


81  10  26  ©IT 


This  Technical  Information  Memorandum  (AFFTC-TIM-81-1) ,  A  Method  of 
Estimating  Upwash  Angle  at  Noseboom-Mounted  Vanes,  was  submitted  under 
Job  Order  Number  SC6342  by  the  Commander,  6520  Test  Group,  Edwards  AFD, 
California  93523. 

This  report  has  been  reviewed  and  cleared  for  open  publication  and/or 
public  release  by  the  AFFTC  Office  of  Public  Affairs  in  accordance  with 
AFR  190-17  and  DODD  5230.9.  There  is  no  objection  to  unlimited  distribu¬ 
tion  of  this  report  to  the  oublic  at  large,  or  by  DDC  to  the  National  Tech¬ 
nical  Information  Service  (NTIS) .  At  NTIS,  it  will  be  available  to  the 
general  public  including  foreign  nationals. 


Prepared  by: 


AFFTC,  Edwards  AFB  CA 


This  report  has  been  reviewed  and  is 
approved  for  publication: 


EDWARD  B.  RUSSELL,  Colonel, 
Commander,  6520  Test  Group 
AFFTC,  Edwards  AFB  CA 


USAF 


When  U.S.  Government  drawinns,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  a  definitely  related  government  procurement  opera¬ 
tion,  the  government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  government  may  have  formulated,  furnished, 
or  in  any  way  supplied  the  said  drawings,  snecif ications ,  or  any  other 
data  is  not  to  be  regarded  by  implication  or  otherwise,  as  in  any  manner 
licensing  the  holder  or  any  other  person  or  corporation  to  convevinn  any 
rights  or  permission  to  manufacture,  use,  or  s^ll  any  patented  invention 
that  may  in  any  way  be  related  thereto. 


Do  not  return  this  copy;  retain  or  destroy 


*  % 


V, 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  THIS  PAGE  (Vh«n  Dele  Entered) 


0 


10 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

L  REPORT  HUMBER 

AFFTC-TIM- 81-1  ^ 

2.  GOVT  ACCESSION  NO. 

A  h-  A^oL 

3.  RECIPIENT'S  CATALOG  NUMBER 

% 00 

A*  TlTi-E  (and  Subtitle 

A  METHOD  OF  ESTIMATING  UPWASH  ANGLE  /ft 

^T^NOSEBOOM- MOUNTED  VA^ES,  ^  1  j 

fy  *  ^ 

3.  TYPE  OF  REPORT  &  PERIOD  COVEREO 

1  FINAL  /■*'£  bit* 

*.  ^RFORMmO  OR«. 

7.  AUTHORS 

KENNETH /RAWLINGS .III 

8.  CONTRACT  OR  GRANT  NUMSERft) 

9  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Air  Force  Flight  Test  Center 

6520  Test  Group/ENDT,  Stop  239 

Edwards  AFB ,  California  93523 

10  PROGRAM  El  LMENT.  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 

JON  SC6342 

1  t.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Air  Force  Flight  Test  Center  Jh 

6520  Test  Group/ENDT,  Stop  239 

Edwards  AFB,  California  93523 

12.  REPORT  DATE 

Juna  Ifrftl  / 

>»■  NU»*«l*OF  paces  '  /'ll  .  .  .  u  O 

14  MONITORING  AGENCY  NAME  A  ADDRESSiT/  dllloront  from  Controlling  Olllce) 

IS.  SECURITY  CLASS,  (of  thtrfSport} 

Unclassified 

IS*  DECLASSIFICATION/ DOWNGRAOING 

schedule 

-i 


l«.  DISTRIBUTION  STATEMENT  (ol  tble  Report) 

THIS  DOCUMENT  HAS  BEEN  APPROVED  FOR  PUBLIC  RELEASE  AND  RESALE; 
ITS  DISTRIBUTION  IS  UNLIMITED. 


17.  DISTRIBUTION  STATEMENT  (of  the  mb, tract  entered  In  Block  20,  II  dllloront  Irom  Report) 


i*  supplementary  notes 


It  KEY  WORDS  fContln ua  on  rovmree  aide  II  neceeeery  and  Identity  by  block  number) 

upwash  fuselage  upwash  angle  of  attack 

upwash  calculation  wing  upwash  angle-of-attack  error 

upwash  calibration  theoretical  upwash  angle-of-attack  calibration 
upwash  estimation 


20.  ABSTRACT  (Continue  on  fri'ani  elde  It  nmcemernry  and  Identity  by  block  number) 

This  technical  information  memorandum  documents  a  theoretically- 
based  method  of  estimating  the  upwash  angle  at  no seboom- mounted 
vanes  of  an  arbitrary  aircraft.  The  method  is  commonly  called 
the  Yaggy-Rogallo  technique  and  has  been  used  successfully  for 
many  years  at  the  NASA/Dryden  Flight  Research  Center  (NASA/DFRC) 
and  the  Air  Force  Flight  Test  Center  (AFFTC; .  The  technique 
divides  any  aircraft  into  a  collection  of  bodies  of  revolution 


/ 


2 


00  ,'STn  1473 


EDITION  OF  I  NOV  f>  IS  OBSOLETE 


UNCLASSIFIED 


ft 

S 


•  t*r 


1 

4 


m 


4 


‘I  I 

m  I 

I 

■S  ! 

I  I 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  Cate  Motored) 


StCUt»  TV  CLASSIFICATION  OF  THII  PAOIflWln  Pit  Kmfn4) 

20.  and  thin  airfoils,  makes  upwash  angle  estimates  for  each 

component,  and  combines  the  individual  estimates  into  a  single 
estimate.  Development  of  the  theory  for  upwash  angle  estimates 
of  both  bodies  of  revolution  and  thin  airfoils  are  presented. 
Recent  improvements  in  solution  and  use  of  the  Yaggy-Rogallo 
equations  are  also  presented.  Guidelines  for  using  the 
upwash  angle  estimates  in  data  reduction  and  analysis  are 
presented  along  with  comparisons  to  flight  test  upwash  calibra¬ 
tion  data  showing  excellent  results.  FORTRAN  V  computer 
programs  were  written  to  simplify  and  expedite  making  upwash 
angle  estimates.  Listings  of  the  program  code,  a  user's 
guide,  aid  a  programmer's  guide  are  included  in  the  memorandum. 


Accer.3i.cn  For 

"NTIS  CRAfcl 
D7IC  TAB 
Unannounc  ct 
Justifies l :  t 


By - - 

Distribution/ 

Availability  Codes 
Avail  and/or 
Dist  I  Special 


UNCLASSIFIED 

SCCUNITV  CLASSIFICATION  of  THIS  FAOIfWN**  Dmim  Kntm4) 


PREFACE 


The  determination  of  upwash  angle  induced  by  the  aircraft  at 
noseboom-mounted  angle-of -attack  vanes  is  a  problem  the  flight  test  com¬ 
munity  faces  on  a  recurring  basis.  The  Air  Force  Flight  Test  Center  (AFFTC) 
r  and  the  NASA  Dryden  Flight  Research  Center  (NASA/DFRC)  use  variations  of  an 

upwash  estimation  technique  developed  at  NASA/DFRC  which  is  commonly  called 
the  Yaggy-Rogallo  technique.  This  technique  has  been  applied  successfully 
to  numerous,  varied  aircraft  by  both  organizations.  Confusion  on  the  part 
of  "first-time"  users  in  applying  the  technique  has  resulted  because  there 
was  little  formal  documentation  of  the  technique  and  most  information  was 
passed  by  informal  contact  with  previous  users.  Since  no  proven  computer 
program  existed,  numerous  procedures  to  implement  the  technique  were 
developed.  The  procedures  ranged  from  hand  calculation  to  use  of  small 
programmable  calculators  to  use  of  desk-top  and  large  mainframe  computers. 
Most  of  the  procedures  were  insufficiently  documented  to  have  any  retained 
value. 

Donald  R.  Bellman  of  NASA/DFRC  was  the  key  individual  in  development 
of  the  Yaggy-Rogallo  upwash  estimation  technique.  It  was  his  insight 
into  the  work  of  Yaggy  and  Rogallo  and  his  expertise  in  translating 
their  work  into  an  easily  used  format  that  led  to  creation  of  this 
technique,  lie  refined  the  technique  and  applied  it  to  many  projects 
only  one  of  which  is  referenced  in  this  memorandum.  The  technique  would 
have  undoubtably  been  lost  due  to  lack  of  use  had  it  not  been  for  Edwin 
J.  Saltzman,  also  of  NASA/DFRC.  His  strong  advocacy  of  the  technique 
and  consultation  with  potential  users  resulted  in  its  continued  use  at 
NASA/DFRC  and  its  introduction  to  AFFTC.  Although  no  original  work  by 
either  Bellman  or  Saltzman  was  available  for  use  directly  in  writing 
this  memorandum,  the  consultation  they  provided  to  other  users  and  myself 
was  indispensable  to  the  creation  of  this  document. 

The  present  effort  was  undertaken  to  consolidate  existing  memos, 
notes,  computer  programs,  and  personal  knowledge  before  they  were  lost.  . 

This  document  is  really  a  collection  of  work  by  numerous  individuals  which 
was  edited  and  significantly  expanded  as  an  expedient  method  of  disseminating 
information.  The  list  of  individuals  whose  work  was  incorporated  in  the 
text  are  recognized  below.  To  provide  a  coherent  and  understandable  doc¬ 
ument,  the  author  modified  and  reorganized  material,  changed  notations, 
and  expanded  many  sections.  Errors,  inconsistencies,  and  misunderstandings 
which  were  introduced  during  this  process  are  totally  the  fault  of  the 
author,  and  in  no  way  reflect  on  the  individuals  who  contributed  to  the 
document.  Those  who  contributed  were: 

Glenn  M.  Sakamoto 
Aerospace  Engineer 
NASA/DFRC,  Edwards  CA 

Steven  J.  Grier 

Cooperative  Engineering  Student 
NASA/DFRC,  Edwards  CA 

Arthur  E.  Sheridan 
Cooperative  Engineering  Student 
NASA/DFRC,  Edwards  CA 


John  W.  Hicks 
Aerospace  Engineer 
AFFTC,  Edwards  AFB  CA 

Patricia  Juenemann 
Data  Technician 
AFFTC,  Edwards  AFB  CA 

Mark  T.  Korsmo 

Cooperative  Engineering  Student 
AFFTC,  Edwards  AFB  CA 
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The  accurate  determination  of  angle  of  attack  is  a  significant  problem 
to  much  of  the  flight-test  community.  Since  angle  of  attack  is  a  primary 
parameter  to  performance,  stability-and-control ,  and  aerodynamic  analyses, 
angle-of -attack  errors  affect  much  flight-test  data.  The  usual  approach  to 
minimizing  these  errors  is  to  obtain  indicated  angle  of  attack  from  the 
angle-of-attack  portion  of  a  noseboom  instrumentation  unit  (NBIU)  mounted 
on  a  flight-test  noseboom.  The  NBIU  normally  has  a  calibration  from  a  wind 
tunnel  which  accounts  for  angle-of-attack  errors  due  to  NBIU  local  flow 
effects.  At  subsonic  Mach  numbers,  however,  there  are  errors  in  indi¬ 
cated  angle  of  attack  due  to  the  fuselage  and  wing  upwash  which  are  present 
even  with  noseboom-mounted  vanes.  Angle-of -attack  errors  caused  by  the 
aerodynamic  influence  of  bodies  results  from  two  -ic.uces ;  firs';,  effects 
of  local  flow  about  the  NBIU  such  as  aerodynamic  interference  boundary 
layer  effects,  and  shock  interaction  and,  second,  upwash  effects  of  air¬ 
craft  components  such  as  the  fuselage  and  wing.  .  >.e  NBIU  local  flow  effects 
can  be  adequately  defined  only  by  wind  tunnel  testing.  The  NBIU  wind- 
tunnel  calibration  curves,  however,  have  an  upwash  component  due  to  the 
noseboom  used  in  the  tunnel  tests  which  must  be  adjusted  to  account  for 
different  noseboom  configurations.  The  upwash  adjustment  for  the  NBIU 
noseboom  as  well  as  upwash  of  the  wing  and  fuselage  can  be  estimated  using 
the  techniques  presented  in  the  memorandum.  An  in-flight  calibration,  how¬ 
ever,  must  be  performed  to  establish  the  angle-of-attack  errors  and  allow 
empirical  adjustment  of  the  upwash  estimates. 

A  good  upwash-estimation  technique  is  important  both  during  applica¬ 
tion  of  the  NBIU  wind-tunnel  calibration  to  flight-test  data  reduction 
and  during  analysis  of  the  in-flight  upwash  calibration.  Adjustments  to 
the  NBIU  calibrations  can  be  adequately  predicted  by  calculation  and 
new  curves  developed  for  use  in  data  reduction  for  removal  of  NBIU  effects 
prior  to  any  analysis  of  fuselage  and  wing  upwash.  Upwash  angle  estimates 
for  the  fuselage  and  wing  are  available  well  before  flight-test  results 
and  are,  therefore,  available  for  initial  prediction.  Characteristics 
of  subsequent  in-flight  calibrations  should  be  consistent  with  the  predicted 
characteristics.  Separation  of  Mach  number  and  lift  effects  on  upwash  angle 
should  be  accomplished  based  on  empirical  adjustment  of  the  level  of  initial 
curves  rather  than  attempts  to  define  characteristics  based  on  the  flight 
test  data.  It  should  be  emphasized  that  the  empirical  adjustment  is  made 
only  to  upwash  estimation  curves  and  not  to  the  NBIU  wind-tunnel  calibra¬ 
tion  curves. 

An  upwash  estimation  technique  which  has  consistently  proven  to  be 
simply  used  and  give  acceptable  results  is  based  on  work  by  Paul  F.  Yaggy 
and  Vernon  L.  Rogallo  and  is  commonly  called  the  Yaggy-Rogallo  technique. 

The  Yagqy-Rogallo  technique  is  based  largely  on  theoretical  techniques 
used  by  Yaggy  and  Rogallo  in  analyzing  their  experimental  investigations 
with  some  empirical  modification  by  Yaggy,  Rogallo  and  Donald  R.  Bellman. 

The  work  of  Yaggy  and  Rogallo  was  refined  and  simplified  by  Donald  R. 

Bellman  of  NASA/Dryden  Flight  Research  Center  (NASA/DFRC)  for  application 
to  flight-test  data.  Although  the  technique  has  been  in  use  for  many  years 
at  NASA/DFRC  and  the  Air  Force  Fliqht  Test  Center  (AFFTC) ,  the  theory, 
procedures,  and  computer  programs  were  never  adequately  documented.  The 
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equations  and  solutions  developed  to  implement  the  technique  were  often 
inadequate  or  had  unique  characteristics  which  restricted  them  from  general 
use.  Terminology  and  notation  varied  enough  to  cause  considerable  confusion 
in  adapting  and  applying  an  existing  implementation  to  subsequent  aircraft. 

The  procedures  for  applying  the  technique  varied  from  hand  calculation  up 
through  a  FORTRAN  IV  computer  program  for  a  large  mainframe  computer.  Docu¬ 
mentation  for  much  of  the  derivation  and  many  of  the  computer  programs  was 
incomplete  or  nonexistent.  The  Background  section  of  this  memorandum  * 

describes  the  development  o*  the  Yaggy-Rogallo  technique  and  reviews  appli¬ 
cable,  previous  work  both  documented  and  undocumented. 

The  Yaggy-Rogallo  technique  treats  any  vehicle  as  a  collection  of 
aerodynamic  bodies  whose  upwash  may  be  calculated  in  one  of  two  ways.  Cyl¬ 
inder-like  bodies  which  are  basically  nonlifting  such  as  nosebooms,  fuse¬ 
lages,  and  stores  are  treated  as  bodies  of  revolution.  Primary  lift-producing 
surfaces  such  as  wings,  canard’*,  and  tails  are  treated  as  thin  airfoils.  A 
concise  method  of  estimating  upwash  for  each  of  the  two  types  of  bodies  has 
been  developed  from  the  work  of  Yaggy  and  Rogallo.  In  the  case  of  both  types 
of  bodies,  assumptions  and  refinements  were  made  to  simplify  application  of 
the  equations  and  expedite  making  upwash  estimates.  The  upwash  for  the 
vehicle  is  then  determined  by  combining  the  estimates  for  the  individual 
parts.  The  Yaggy-Rogallo  Analysis  section  of  this  memorandum  describes  the 
procedure  for  dividing  a  vehicle  into  components  and  determining  equivalent 
body  shapes.  The  equations  for  calculating  upwash  of  the  components  are 
described  as  is  the  procedure  for  combining  the  component  upwash  estimate's 
into  a  single  estimate  for  the  vehicle.  Detailed  development  of  the 
equations  and  notation  for  upwash  angle  estimates  of  bodies  of  revolution 
and  thin  airfoils  are  presented  in  Appendices  A  and  B  respectively.  For 
users  who  are  familiar  with  the  Yaggy-Rogallo  technique,  a  procedural 
outline  is  presented  in  Appendix  C  for  use  as  a  quick  reference  and  check¬ 
list  in  making  upwash  angle  estimates.  The  equations  and  procedures  were  , 

implemented  in  two  FORTRAN  V  computer  programs  which  are  well  documented 
in  a  user's  guide  and  programmer's  guide,  Appendices  D  and  E  respectively. 

Practical  application  of  the  current  implementation  of  the  Yaggy-Rogallo 
technique  and  corresponding  computer  programs  is  -rhown  in  the  Procedural  Ap¬ 
plications  section  of  this  memorandum.  Two  examples  of  obtaining  upwash 
angle  estimates  from  Yaggy-Rogallo  analysis  are  presented.  The  examples 
show  how  the  aircraft  is  broken  into  components  and  data  is  obtained  for 
entry  into  the  production  software.  The  card  input  and  printed  output  are 
presented  to  show  use  of  the  software  and  the  plotted  results  are  presented 
to  show  expected  characteristics.  The  results  are  applied  to  actual  flight- 
test  data  to  demonstrate  the  excellent  results  which  can  be  expected  from 
use  of  the  Yaggy-Rogallo  technique  in  flight-test  data  reduction. 

BACKGROUND 

The  present  upwatih  estimation  technique  evolved  through  successive 
refinements  by  numerous  persons  and  organizations.  The  original  technique 
was  developed  by  Paul  F.  Yaqgy  in  Reference  1  and  greatly  generalized  by 
Vernon  L.  Rogcllo  in  Reference  2.  The  technique,  which  became  known  as 
the  Yaggy-Rogallo  technique,  represented  any  vehicle  as  a  collection  of 
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1”^  fo?S.  «“f*cUnt  data  for  airship  hulls  from  potential  How 

/oiutions  for  b-i^ 

a  simple  method  ot  “J“^w^;oi»«lyai5iuiSrLlSMinTN' 1  s  1476,  1491, 
erences  4  and  5  (Reference  a  tiv*lv) .  Using  a  modified  liftinq-line 

^  £.»”?'»!».«  «*»•>>  «e.u.»t 

Paul  F.  Yaggy  used  Von  Kansan's  work  and  basic  liftinq-line  theory  in 
a  study  of  upwash  angle  induced  at  the  propeller  planes  of  multiengine, 

winq-mounted^engines  in  R.f ar.no.  1.  In  order  to  optimrre  J**" 

/*  , _ .,r.Seac  it  was  required  to  estimate  tne  upwasn  anyic 

induced0 by  theming  and  nacelle.  Yaqqy  modified  Von  Karman;s  development 
tests  fro^NACA/Ames11  Aeronautical1  Laboratory3  ( later  NASA/Ames  Research 

i;  S 

inq  airflow  didn't  materially  effect  results  even  though  the  assumptions 
£f9potential  flow  were  violated.  Yaggy's  work  was  confined  to  unswept  wings 
and^incompressible  flow  so  that  basic  liftinq-line  theory  gaYe  g°°d  **®^tS* 

He  did,  however,  recognize  the  deficiencies  and  recommended  lnc°rP°™^ng 
John  ne  Young's  work  if  if  forts  were  expanded  to  swept  wings  and  higher  Mach 

numbers. 

Vernon  L.  Rogallo  expanded  Yaggy's  work  to  propeller  planes  of  aircraf^ t 
with  sweet  winqs  and  applied  Mach  corrections  to  wing  upwash  effects. 

Ref erenc^ 2,  Rogallo  described  a  means  of  calculating  Mach  effects  by  increas¬ 
ing  the  effective  distance  of  the  wing  aft  of  the  propeller  plans  and  in 
creasing  the  effective  wing  sweeps.  Although  Rogallo  added  Ma°J utioJ^for 
the  wing  component  of  upwash,  he  used  Yaggy  s  incompressible  solution  for 
nacelles.  Rogallo  comnared  his  theoretical  results  with  experimental  data 
with  good  agreement. 

The  concepts  used  by  Yaqqy  and  Rogallo  were  expanded  to  upwash. 
tion  at  noseboom-mounted  angle-of-attack  vanes  at  NA5A/DFRC.  Donald  R.  Be  li¬ 
man  and  Edwin  J.  S»ltzman  modified  the  technique  to  account  for  differen  e 
between  the  propeller  'nd  noseboom  applications.  Reference  9  documents  that 

compressibility  corrections  were  added  to  the  body-of mnlif^and 
tions.  The  technique  for  lifting  surfaces  was  also  refined  to 
expedite  application  with  little  reduction  m  accuracy.  By  certain  simpli¬ 
fying  assumptions,  the  independent  parameters  used  were  changed  to  geometric 
parameters  only  and  the  curves  were  expanded  to  a  range  of  distances  ahead 
of  the  wing  which  would  include  angle-of-attack  vanes.  The  modified  tech¬ 
niques  were  applied  to  a  noseboom  calibration  in  Reference  10  and  to  the 
M2-F1,  HL-10 ,  and  M2-F2  aircraft  in  References  11,  12,  and  13,  * 

with  good  comparison  with  full  scale  test  of  the  vehicles  m  NASA/ARC  wind 

tunnels . 
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LIST  OF  SYMBOLS  AND  ABBREVIATIONS 


ITEM 

DESCRIPTION 

USCU1 

UNITS 

SI1 

UNITS 

SYMBOLS 

A 

body  cross-sectional  area 

.  2 
m 

2 

m 

avn 

Weissinger  influence  coefficient  for 
symmetric  loading 

N-D 

N-D 

AR 

2 

aspect  ratio  of  a  lifting  surface,  AR=b  /S 

N-D 

N-D 

b 

span  of  a  lifting  surface 

ft 

m 

c 

local  chord  length  of  a  lifting  surface 

ft 

m 

c 

mean  aerodynamic  chord  length  of  a  lifting 
surface 

ft 

m 

CL 

lift  coefficient  of  a  lifting  surface, 
CL=L/qS 

N-D 

N-D 

section  lift  coefficient  of  a  lifting 
surface,  C^VqS 

N-D 

N-D 

% 

lift-curve  slope 

/deg 

/deg 

C^c 

span  loading  coefficient 

N-D 

N-D 

V 

d 

longitudinal  distance  ahead  of  angle-of- 
attack  vane 

in 

m 

dx 

element  of  the  longitudinal  axis 

in 

m 

FS 

fuselage  station 

in 

m 

G 

n 

dimensionless  circulation  about  the  wing 
at  span  station  n 

N-D 

N-D 

i 

incidence  angle 

deg 

deg 

constant,  see  equation  (A32) 

in 

in 

*U.S.  Customary  Units  (USCU)  and  International  System  of  Units  (SI)  were 
used  interchangeably  during  analyses  in  this  report.  Analyses  were  based 
on  ratios  and  non-dimensional  parameters  as  much  as  possible  so  that  con¬ 
sistency  within  systems  is  much  more  important  than  the  system  itself. 
Conversion  factors  between  systems  were  obtained  from  Reference  14  when 
required. 
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LIST  OF  SYMBOLS  AND  ABBREVIATIONS  (CONTINUED) 


CTEM 

DESCRIPTION 

USCU 

UNITS 

SI 

UNITS 

SYMBOLS 

K2 

constant,  see  equation  (A32) 

in 

m 

constant,  see  equation  (A40) 

in2 

2 

m 

constant,  see  equation  (A40) 

in2 

2 

m 

kav 

average  of  ratios  of  actual  sectional  lift- 
curve  slopes  to  theoretical  values 

N-D 

N-D 

L 

lifting  surface  lift 

lbf 

n 

l 

section  lift 

lbf 

n 

l 

general  unit  of  length  (see  footnote  3, 
page  89) 

— 

— 

M 

Mach  number 

N-D 

N-D 

N 

number  of  points 

— 

— 

P 

arbitrary  point  in  flowfield  of  a  body 
revolution 

— 

— 

q 

dynamic  pressure 

lb/ft2 

n/m2 

R 

equivalent  radius  of  a  body,  R  =  /a/ti 

in 

m 

r 

perpendicular  distance  from  longitudinal 
axis  to  center  of  pressure  of  angle-of- 
attack  vane 

in 

m 

S 

reference  area  of  a  lifting  surface 

ft2 

2 

m 

U 

longitudinal  velocity 

ft/sec 

m/sec 

V 

airspeed 

ft/sec 

m/sec 

w 

upwash  velocity 

ft/sec 

m/sec 

X,Y.  Z 

general  rectangular  coordinate  system 

— 

— 

i 

angle  of  attack 

deg 

deg 

B 

compressibility  parameter, 

N-l) 

N-D 

r 

c;rculation  about  the  wing 

,„2  . 
ft  /sec 

2  . 

m  ,  sec 

i 

upv.’ash  angle,  *  =sin  ^  (W/V^) 

deg 

deg 

it 


LIST  OF  SYMBOLS  AND  ABBREVIATIONS  (CONTINUED) 


ITEM 

DESCRIPTION 

uscu 

UNITS 

SI 

UNITS 

SYMBOLS 

0 

polar  coordinate  angle,  0  =  cot  ^(d/|r|) 

rad 

rad 

A 

sweep  angle  of  a  lifting  surface 

deg 

deg 

V 

moment  per  unit  length  of  a  doublet  element 

lbf 

n 

n 

constant,  3.14159 

N-D 

N-D 

T 

dimensionless  distance  from  quarter  chord 
to  center  of  pressure  of  angle-of-attack 
vane,  t  =  d/(b/2) 

N-D 

N-D 

4 

potential  function 

N-D 

N-D 

n 

angular  location  of  the  center  of  pressure 
of  the  angle-of-attack  vane 

deg 

deg 

SUPERSCRIPTS 

~ 

value  corrected  for  compressibility 

— 

— 

SUBSCRIPTS 

0 

zero  -  lift 

— 

— 

C/4 

quarter  chord 

— 

— 

F 

fuselage 

— 

— 

i 

the  ith  point  along  a  body 

— 

— 

LE 

leading  edge 

— 

— 

LF 

local  flow 

— 

— 

NB 

noseboom 

— 

— 

n 

spanwise  station  influencing  upwash  at 
station  v 

— 

— 

r 

radial  direction 

— 

— 

V 

center  of  pressure  of  angle-of-attack  vane 

— 

— 

V 

spanwise  station  where  upwash  estimates 
can  be  calculated 

— 

— 

T 

true 

— 

— 

TE 

trailing  edge 

12 

— 

— 

i 

i 


ifirirtfch* 


V 


H 


, 

LIST 

OF  SYMBOLS  AND  ABBREVIATIONS 

(CONCLUDED) 

% 

USCU 

SI 

'  h 

t. 

ITEM 

DESCRIPTION 

UNITS 

UNITS 

> 

SUBSCRIPTS 

r 

W 

wing 

— 

— 

ft 
'  1 

6 

value  corrected  for  compressibility 

— 

— 

i 

3if 

12 

tangential  direction 

— 

— 

i 

ABBREVIATIONS 

pi 

i 

AFFTC 

Air  Force  Flight  Test  Center 

— 

— 

1 

AWACS 

Airborne  Warning  and  Control  System 

— 

— 

jj 

& 

FORTRAN 

FORmula  TRANslation  programming  language 

— 

— 

NASA/ARC 

NASA/Ames  Research  Center 

— 

— 

i* 

-3 

NASA/DFRC 

NASA/Dryden  Flight  Research  Center 

— 

— 

NBIU 

noseboom  instrumentation  unit 

— 

— 

Saab 

Saab-Scania  AB 

— 

— 

SI 

International  System  of  Units 

— 

— 

TACT 

Transonic  Aircraft  Technology  Program 

— 

— 

USCU 

U.  S.  Customary  Units 

— 

— 

n 

| 

WINGLETS 

USAF/NASA  KC- 1 35/WINGLETS  Program 

— 

— 

i  i 

3 

i 


’i 
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YAGGY-RRGRLLO  ANALYSIS 


The  application  of  the  Yaggy-Rogallo  technique  to  flight  test  vehicles 
requires  a  functional  and  geometric  analysis  of  the  aircraft.  The  Yaggy- 
Rogallo  technique  allows  analysis  of  two  functional  types  of  bodies.  The 
actual  test  aircraft  must  be  represented  by  a  collection  of  bodies  of  rev¬ 
olution  and  thin-airfoil  type  lifting  surfaces.  Bodies  which  are  not  pri¬ 
marily  lifting  surfaces  such  as  nosebooms,  fuselages,  wing-mounted  nacelles, 
and  external  stores  are  treated  as  bodies  of  revolution.  Conventional  wings, 
canards,  and  horizontal  tails  (if  treated  at  all)  are  treated  by  thin  air¬ 
foil  theory.  The  analysis  of  the  actual  vehicle  and  division  into  compo¬ 
nents  is  the  first  and  most  important  step  in  Yaggy-Rogallo  analysis. 

Most  vehicles  are  easily  analyzed  and  the  definition  of  the  components 
rather  straightforward.  There  are,  however,  many  new  and  unique  aircraft 
that  are  not  easily  analyzed.  Bodies  such  as  the  Airborne  Warning  and  Con¬ 
trol  System  (AWACS)  fuselage-mounted  radome,  for  instance,  are  not  easily 
categorized.  More  common  are  aircraft  such  as  the  F-16  and  F-18  which  have 
strakes  extending  forward  along  the  side  of  the  fuselage.  The  strakes  are 
not  efficient  lifters  at  low  angles  of  attack  and  have  such  a  low  aspect 
ratio  that  they  are  not  easily  handled  by  thin-airfoil  theory.  If  the 
strakes  are  included  in  the  fuselage  analysis,  they  cause  much  greater 
apparent  fuselage  area  than  would  be  obtained  from  the  equivalent-circular- 
body  analysis  described  in  the  following  sections.  The  "chine"-type  fuse¬ 
lage  such  as  the  SR-71  present  the  same  problem  to  a  greater  degree.  Al¬ 
though  no  quantitative  guidance  could  be  provided,  it  appears  that  some 
combination  of  increasing  the  effective  diameter  of  the  fuselage  and  treat¬ 
ing  them  as  a  very  low  aspect  ratio  canard  would  be  appropriate.  Another 
problem  arises  from  wings  such  as  those  of  the  Saab  JA37  "Viggen"  which  have 
changes  in  leading  edge  sweep  with  spanwise  station.  Since  the  theory  dis¬ 
cussed  in  later  sections  assumes  a  straight  quarter-chord  line,  it  is 
necessary  to  define  an  equivalent  wing  with  slightly  different  geometry  and 
straight  quarter-chord.  These  examples  demonstrate  the  need  for  a  careful 
analysis  of  the  vehicle  and  the  possibility  of  the  need  to  change  the  geo¬ 
metric  characteristics  to  get  an  equivalent  aerodynamic  representation. 

The  upwash  of  each  component  obtained  from  the  previous  analysis  is 
obtained  from  an  independent  analysis  of  upwash  and  the  results  summed. 

In  most  cases  the  results  are  summed  with  no  regard  for  interference 
effects.  Where  required,  there  are  rather  simple  means  of  accounting  for 
first-order  interference  effects  which  have  been  proven  with  test  data. 
However,  for  analysis  of  numerous  bodies  which  are  in  close  proximity  and 
with  larqe  interference  effects,  there  is  no  adequate  theory  for  inter¬ 
ference  effects.  Care  should  be  taken  in  extending  the  following  analyses 
to  collections  of  bodies  where  interference  would  be  a  dominate  effect. 


CYLINORICRL  COHPOHCMTS 

The  upwash  estimates  for  non-lifting  bodies  which  the  functional  and 
geometric  analysis  has  shown  must  be  treated  as  bodies  of  revolution  are 
obtained  from  equations  developed  in  Appendix  A.  An  equivalent  body  of  rev¬ 
olution,  as  shown  in  Fiqure  1,  must  be  developed  for  each  separate  component 
which  must  be  analyzed.  Bodies  whose  longitudinal  axes  are  not  concurrent, 
such  as  "drooped"  nosebooms,  must  be  analyzed  separately  since  the  vane 
position  relative  to  the  axis  is  different  and  the  bodies  have  different 
local  angles  of  attack.  Also  bodies  having  different  local  angle  of  attack 
due  to  induced  angle  of  attack  from  nearby  bodies  must  be  analyzed  separately. 
A  coordinate  transformation  must  also  be  made  so  that  distances  are 

M 


B)  EQUIVALENT  BODY  OF  REVOLUTION 

FIGURE  1«  REPRESENTATION  OF  NON-LIFTINO  BODIES 

J 


measured  forward  and  aft  of  the  angle-of-attack  vane  as  shown  in  Figure  2. 
Once  the  equivalent  geometry  has  been  developed,  compressibility  effects 
can  be  accounted  for  by  further  changes  to  an  "effective"  geometric  config¬ 
uration.  These  changes  allow  upwash  estimates  to  be  calculated  using 
equations  developed  in  Appendix  A. 

Experience  has  shown  that  bodies  with  significantly  different  cross- 
sectional  shapes  can  be  represented  by  an  equivalent-body-of-revolution 
analysis  with  little  effect  on  upwash  estimates.  The  body  of  revolution 
is  obtained  by  calculating  a  radius  for  a  circle  with  equal  area  to  the 
actual  body.  The  radius  of  the  equivalent  circle,  R,  is  related  to  the 
cross  sectional  area,  A,  by  equation  (1).  By  calculating  the  equivalent 
radius  at  enough  sections  from  plotted  or  tabulated  data  the  body  of  rev¬ 
olution  can  be  determined.  The  centers  of  the  circles  lie  on  the  original 
longitudinal  axis  and  the  radii  are  perpendicular  to  it.  In  determining 
the  area  of  a  section,  the  area  attributable  to  lifting  surfaces  such  as 
wings  or  tails  is  ignored  except  in  unusual  cases  such  as  strakes  which 
were  discussed  previously. 
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The  location  of  the  center  of  pressure  of  the  nos eboom-mounted  vane 
relative  to  the  centerline  of  each  component  must  be  specified.  The  radial 
distance,  r,  and  the  angular  position,  Q,  are  a.  shown  F^ure 

2.  It  is  important  in  establishing  the  value  of  Q  to  realise  that  the  tech 
nique  is  extremely  sensitive  to  the  value  of  R  and  ^hat,  in  many  cases, 
there  is  significant  uncertainty  in  its  value.  In  most  cases,  the  vertical 
location  of  the  axis  of  the  equivalent-circular  body  representing  a  fuselage 
is  SiknSwn  by  several  inches  or  more.  For  the  AFFTC  standard  noseboom 
instrumentation  unit,  a  +7  inch  uncertainty  could  vary  the  value  of  n  from 
458  £nd  135^  and  vary  the  upwash  estimate  from  zero  to  its  maximum  and  back 
to  ,“o.  Wh»  th.rAs  uncertainty,  it  i.  beet  to  ».k.  initial' 
with  an  a  value  of  90f  and  adjust  the  value  based  on  an  inflight  calibration. 
However ,  where  the  body  actually  is  cylindrical  and  tjj 

well  known  such  as  a  probe  refueling  receptical,  the  calculated  value  or 
R  should  be  used.  Distances  along  the  longitudinal  axis,  d^  are  “®aJ“  d 
positive  forward  of  a  plane  perpendicular  to  the  longitudinal  axis  whicn 
contains  the  center  of  pressure  of  the  angle  of  attack  vane  as  stances 

Figure  2.  To  account  for  compressibility,  these  actual 

must  be  adjusted  by  a  compressibility  factor,  B,  defined  by  equation  (2). 

The  effect  of  compressibility  is  to  create  an  effective  geometry  such 
that  sections  forward  of  the  vane  appear  closer  to  the  vane  and  ®®^°ns 
aft  of  the  vane  appear  farther  away.  This  continues  until  Mach  number,  M, 
of  1.0.  Beyond  »toch  1.0  the  body  has  practically  no  effect  and  “P"a®Jf*®tive 
considered  to  go  to  zero.  For  Mach  numbers  between  0.0  and  1.0  the  effecti 
distances  are  defined  by  equations  (3).  The  Yaggy-Rogallo  equation  shown 

as  equation  (4)  was  developed  in  Appendix  A  and  is  used  f ji^b^equation 
upwash  angle,  e,  per  angle  of  attack,  a.  The  angle  B  is  defined  by* “JJ^ion 

(5)  as  explained  in  footnote  2  (following  page). 

formulation  for  B  into  equation  (4) ,  it  can  be  completely  defined  and 
solved. 


d  =  d»B 


d  >  o  (forward) 
d  <  o  (aft) 


c  sin  r. 


2„  f  6TF  , 

/  H2si 

2rz  J  : 


sin6  d0 


9  =  cot  ( d/ 1 r ' ) 


O  <G  <  IT 
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FIGURE  3*  IgmgWTRTIOH  Of  BOOIEI  OF  REVOLUTION 


Several  solutions  to  equation  (4)  are  discussed  in  Appendix  A  and  the 
most  practical  was  determined.  The  recommended  procedure  is  determination 
of  an  upwash  estimate  through  summation  of  upwash  estimates  for  individual 
segments  of  the  body  when  summed  from  leading  to  trailing  edge.  The  body 
is  divided  into  N-l  segments  by  N  points  where  point  1  is  the  leading  edge 
and  point  N  is  the  trailing  edge  as  shown  in  Figure  3.  The  length  and  num¬ 
ber  of  segments  should  be  selected  such  that  the  area  between  points  varies 
linearily,  or  approximately  so,  between  points.  Normally,  this  is  accom¬ 
plished  by  selecting  points  from  the  fuselage  area  versus  fuselage  station 
plot  in  such  a  manner  that  a  straight  line  between  points  closely  approx¬ 
imates  the  curve.  The  fuselage  stations  must  be  transformed  into  distances 


2  —1 

The  equation  0  =  cot  (d/|r|)  is  the  necessary  representation  of  the  angle  0 
As  the  value  of  d  moves  from  far  upstream  (d  =  +«)  to  far  downstream  (d  =  -») 
the  angle  8  moves  from  0  to  »  as  shown  in  the  sketch  of  the  cotangent  func¬ 
tion.  However,  most  computers  do  not  allow  the  inverse  cotangent  function 
and  the_^nverse  tangent  function  must  be  used  instead.  The  equation 
0  =  tan  x(]r|/d)  for  0  £  0  £  it  is  equally  valid  as  the  cotangent  function. 
Most  computers,  however,  return  inverse  tangent  values  between  — ti/2  and 
tt/2  since  this  avoids  the  discontinuity  at  0  =  tt/2  as  shown  on  the  sketch 
of  the  tangent  function.  In  this  case  the  equation  must  be  implemented  as 


0  =  tan  * ( | r |/d) 

O  =  tan  1  ( |  r  |  /d )  +  it 
so  that  0  will  move  from  0  to  tt  as  d 


Y=C0T( X ) 


(r/d)>0 
( | r |/d) ^0 
from  +°°  to  -«. 
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forward  or  aft  of  the  center  of  pressure  of  the  vane  using  equation  (6). 

The  corresponding  fuselage  area  must  be  transformed  to  the  radius  of  the 
equivalent  circle  using  equation  (7) .  The  distances  must  be  converted  to 
effective  distances  for  each  Mach  number  where  upwash  estimates  are  desired 
using  8  calculated  using  equation  (2).  The  effective  distances  and  cor¬ 
responding  effective  angles  are  defined  in  equations  (8)  and  (9)  respec¬ 
tively.  Upwash  estimates  in  terms  of  e/a  can  then  be  calculated  using 
equation  (10).  The  calculations  in  equations  (8),  (9),  and  (10)  must  be 
repeated  for  each  Mach  number,  thus  8  value,  where  upwash  estimates  are 
desired.  This  process  is  demonstrated  fully  in  the  Procedural  Applications 
section  of  this  memorandum. 


d.  =*  FS  -FS. 
l  vi 


(6) 


df  =  d^*8 


d .  =  d .  /  3 
l  r 


(7) 


d  >  o  (forward) 
d  <  o  (aft) 


(8) 


0  . 

l 


cot-1 (di t | r| ) 


(9) 


2  2 

-  =  sin  F-cos  71 

a  “  2 

2r 


(Kj (sin(K-sinO^_^)  +  K.,  (cos0^_^-cos0^)  ] 


i=2 


(10) 


"  2  2 

where  =  *Ri-l“Ri*  /  (cotS  ^-cotfl ^ ) 


Kj  *  -K^cotC ^ 


The  summation  procedure  is  very  repetitive  and  becomes  tedious  when 
many  upwash  estimates  at  several  Mach  numbers  must  be  made.  A  FORTRAN  V 
computer  program  has  been  written  which  implements  this  procedure.  It  is 
described  fully  and  in  detail  in  Appendix  D. 


LIFTING  SURFACE  COMPONENTS 


The  upwash  estimate*  for  components  to  be  treated  as  thin-air foils, 
as  shown  by  functional  and  geometric  analysis,  are  obtained  from  equa¬ 
tions  developed  in  Appendix  B.  The  lifting  surfaces  must  be  analyzed  and, 
if  necessary,  equivalent  planforms  must  be  developed  and  equivalent  geometric 
characteristics  such  as  area,  quarter-chord  sweep,  aspect  ratio,  and  span 
be  determined.  From  these  geometric  values  the  upwash  estimates  for  the 
lifting  surfaces  can  be  determined. 

Experience  has  shown  that  unless  the  planform  of  the  lifting  surface 
differs  significantly  from  usual  configurations,  the  upwash  can  be  deter¬ 
mined  from  curves  or  the  equation  developed  in  Appendix  B.  The  develop¬ 
ment  assumes  that  all  lifting  surfaces  be  in  the  horizontal  plane  of  the 
angle-of-attack  vanes  and  that  the  vanes  are  located  midspan  of  the  lift¬ 
ing  surface.  Thus,  the  location  of  the  vanes  relative  to  the  lifting 
surface  are  specified  by  the  nondimensional  distance,  t,  defined  by 
equation  (11).  The  effective  wing  sweep,  Ag,  can  be  calculated  by  using 
equation  (12).  These  values  can  be  used  inpequation  (13)  or  the  curve 
in  Figure  4  to  obtain  the  upwash  parameter  eAR/C..  The  upwash  parameter  for 
lifting  surfaces  e/ CL  can  then  be  determined  by  dividing  eAR/C  by  aspect 
ratio.  It  is  highly  recommended  that  the  estimates  be  implemented  in  terms 
of  c/CL  for  data  reduction  and  analysis.  When  this  is  not  possible,  it 
may  be^necessary  to  accept  increased  uncertainty  and  modify  the  form  of  the 
correction  as  described  in  the  Combination  of  Component  Estimates  section. 


T  -  <FSv-FSc/4)  /  (b/2) 

! 


(11) 


tan  1  (tanAc^4/R) 


(12) 


f.  AR/C.  -  10lMog<T/?)  *  +  C] 

L 

where  A  ■  -1.488973010 
b  =  -0.008447888 
C  *  -1 .099388684 


(13) 


COMBINATION  OF  COMPONENT  ESTIMATE# 

The  upwash  estimates  from  various  components  must  be  combined  during 
data  reduction  and  analysis  to  give  a  total  estimate  of  the  upwash  angle. 
The  combination  of  upwash  estimates  should  normally  be  a  straightforward 
process  of  combining  the  estimates  for  the  non-lifting  surfaces,  combin¬ 
ing  the  estimates  for  lifting  surfaces,  and  adding  the  results.  The 

21 


o  o  ooo  o—  n  r»  ■» u>« 
T/e  (SEN [SPANS) 


o  o  ooc  oo 
m  n  «io«  • 


FI DUPE  4 i  UPMASH  EST [NATION  CUAVES  FOA  UFTINO  SURFACES 


process  can  become  more  complicated  if  it  is  determined  to  be  necessary 
to  include  the  first  order  effects  due  to  interference.  The  process  becomes 
still  more  complicated  and  significantly  less  accurate  if  the  user  chooses 
to  modify  the  form  of  the  lifting  surface  estimates  from  values  of  e/C* 
to  values  of  e/a.  The  values  of  lift-curve  slope  and  zero-lift  angle  Of  at¬ 
tack  vary  greatly  across  the  lift  coefficient  range  and  may  introduce  large 
errors  when  used  during  the  modification.  Cars  should  be  taken  during  the 
combining  process  to  insure  the  best,  practical  result  is  obtained. 

COm  1  NATION  OF  MQN-LIFTIMO  COgONCW  EBTIWATIS: 

The  combination  of  estimates  for  non-lifting  bodies  which  are  treated 
as  bodies  of  revolution  can  normally  be  combined  without  correction  for 
interference  effects.  Estimates  for  a  noseboom  (NB)  and  fuselage  (F) ,  for 
instance,  can  be  combined  using  equation  (14).  The  wing  angle  of  attack,  a, 
must  be  corrected  for  the  wing  and  noseboom  incidence  angles,  i  and  i 
respectively.  This  is  the  recommended  procedure  for  most  applications. 


c 


(c/a).w(a-iNB)  +  (c/a)p(a-iw) 


(14) 


The  1st  order  effects  of  interference  can  be  accounted  for  by  a  simple 
correction  discussed  in  Reference  1.  The  angle  of  attack  of  the  various  com¬ 
ponents  can  be  adjusted  for  upwash  of  components  downstream  of  the  body  of 
interest.  The  angle  of  attack  of  the  noseboom,  for  instance,  could  be  in¬ 
creased  by  the  induced  angle  of  the  wing  and  fuselage  at  the  noseboom.  Al¬ 
though  the  correction  is  rather  simple,  it  becomes  complicated  to  implement 
if  there  are  several  bodies  and  adds  little  to  the  accuracy  of  the  method 
since  first  order  interference  effects  are  second  or  third  order  effects 
on  upwash  angle  estimates.  Although  the  method  has  proven  acceptable  for 
a  few  bodies,  the  correction  may  be  inadequate  for  large  interference 
effects  of  multiple  bodies  such  as  stores  and  external  tanks. 

COMBINATION  OF  LIFTING  SURFACE  EBTIMATEI: 

The  combination  of  estimates  for  lifting  surfaces  is  far  more  difficult 
than  for  non-lifting  bodies.  If  only  the  wing  is  involved,  the  entire  air¬ 
craft  lift  can  be  attributed  to  the  wing  and  the  upwash  is  simply  determined 
from  equation  (15).  If  two  lifting  surfaces  are  involved  such  as  a  wing  (w) 
and  canard  (c) ,  the  equation  would  be  equation  (16).  Normally  it  is  diffi¬ 
cult  to  accurately  divide  the  lift  between  the  wing  and  canard.  If  it  is 
possible,  it  should  be  remembered  that  the  lift  coefficients  are  based  on 
the  area  and  dynamic  pressure  at  the  surface. 


E  -  <C/CL)W.CL 


(15) 


c 


+  (e/CL,C*% 


(16) 
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A  simpler,  but  significantly  less  accurate,  method  is  to  convert  the 
lifting  surfaces  estimates  to  an  c/a  format.  The  values  of  lift-curve 
slope  and  zero-lift  angle  of  attack  vary  greatly  across  the  lift  coeffi¬ 
cient  range  and  may  introduce  large  errors  in  upwash  angle  estimates. 

This  is  especially  true  at  high  angles  of  attack  where  upwash  angles  are 
large  and  the  lift-curve  slope  is  significantly  different  than  an  average 
value  for  moderate  lift  coefficients.  The  upwash  estimate  in  terms  of  c/a 
can,  however,  be  determined  from  e/CT  and  the  lift  curve  slope,  CT  ,  at  the 

Lo 

appropriate  Mach  number.  If  the  angle  of  attack  at  zero  lift  is  a.,  the 
upwash  angle  is  given  by  equation  (17).  This  method  is  not  recommended  and 
should  be  avoided  if  possible. 


t  =  (£/Cl)*Cl 


<a-aQ) 


((e/Cl)*Cl  a)  -  ((e/C  ) 
a 


CL  V 
a 


(17) 


PROCEDURAL  APPLICATIONS 


The  Yaggy-Rogallo  technique  was  used  to  make  upwash  estimates  for  two 
common  flight- test  applications.  The  applications  were  presented  primarily 
to  demonstrate  the  procedures  to  be  followed  in  developing  upwash  estimates 
and  applying  them  to  flight-test  data  reduction.  The  results,  however,  also 
demonstrated  the  excellent  quality  of  the  upwash  estimates  and  the  versatility 
of  the  Yaggy-Rogallo  technique.  The  first  application  discussed  is  develop¬ 
ment  of  total  upvash  estimates  for  use  in  reduction  of  data  from  the  F-lllA/ 
TACT  aircraft.  Development  of  the  F-lllA/TACT  upwash  estimates  was  a  complete 
demonstration  of  the  Yaggy-Rogallo  technique  and  showed  the  excellent  abso¬ 
lute  level  of  curves  obtained  using  the  technique.  The  second  application 
discussed  is  development  of  corrections  to  upwash  estimates  to  account  for 
changing  noseboom  configuration.  This  application  demonstrated  the  quality 
of  relative  values  between  upwash  estimates  for  two  configurations  and  showed 
the  versatility  of  the  technique.  Both  applications  are  documented  in  the 
following  subsections. 

F-lllR/TACT  UPWASH  ESTIMATES 

The  Yaggy-Rogallo  technique  was  applied  to  the  supercritical  wing  con¬ 
figuration  of  the  F-lllA/TACT  aircraft  to  completely  demonstrate  a  typical 
application  of  the  technique.  This  aircraft  was  selected  because  the  vari¬ 
able  wing  sweep  and  large  Mach  number  range  of  the  aircraft  allowed  maximum 
demonstration  and  verification  cf  the  technique.  Also,  a  large  quantity  of 
high  quality,  angle-of-attack-calibration  data  was  available  for  comparison. 
Upwash  estimates  were  developed  for  three  wing  sweep*  and  for  the  Mach 
range  from  0.0  to  1.0.  The  development  followed  procedures  recommended  in 
the  Yaggy-Rogallo  Analysis  section  of  this  memorandum  entirely  as  well  as 
using  some  alternative  procedures  necessary  to  convert  the  wing  upwash 
estimates  from  c/C  to  c/a  format. 

Li 
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The  process  of  making  comprehensive  upwash  estimates  began  with  a 
preliminary  analysis  of  the  test  aircraft.  As  discussed  in  the  Yaggy-Rogallo 
Analysis  section  of  this  report,  it  was  necessary  to  separate  the  test 
aircraft  into  a  collection  of  bodies  to  be  represented  either  as  bodies  of 
revolution  or  thin-airfoil  lifting  surfaces.  Analysis  of  the  F-lllA/TACT 
aircraft  revealed  that  only  one  body  of  revolution  and  three  lifting  sur¬ 
faces  were  necessary  to  represent  the  aircraft  adequately  for  making  up¬ 
wash  estimates.  This  was  due  in  part  to  choosing  an  aircraft  equipped  with 
the  AFFTC  standard  noseboom  instrumentation  unit  (NBIU) .  The  NBIU  has 
unique  angle-of-attack  corrections  dominated  by  local  flow  effects  which 
can't  possibly  be  estimated  using  the  Yaggy-Rogallo  technique)  the  wind 
tunnel  calibration  of  the  NBIU  must  be  used  to  correct  for  local  flow  ef¬ 
fects.  The  wind  tunnel  calibration,  however,  also  contains  correction  for 


the  upwash  of  the  noseboom  used  in  the  wind  tunnel  tests  which  normally 
must  be  corrected  as  demonstrated  in  he  following  subsection  entitled 
Noseboom  Correction  Estimates.  The  F-lllA/TACT  aircraft  was  equipped  with 
a  noseboom  which  exactly  matched  the  wind  tunnel  configuration  so  that  anal¬ 
ysis  of  the  noseboom  was  unnecessary.  As  with  most  conventional  tails, 
the  upwash  of  the  tail  was  insignificant  and  could  be  ignored.  Representing 
the  aircraft  was  thus  reduced  to  representing  the  fuselage  by  an  ejuivalent 
body  of  revolution  and  each  of  the  three  wingsweeps  by  a  thin  airfoil.  Up¬ 
wash  estimates  were  prepared  for  the  body  and  each  airfoil  using  procedures 
recommended  in  the  Yaggy-Rogallo  Analysis  section  of  this  memorandum.  The 
upwash  estimates  were  prepared  using  the  production  software  described  in 
Appendix  D  and  are  presented  as  production-software  output  although  the 
same  results  could  have  been  obtained  by  hand  calculation  if  desired.  The 
resulting  upwash  estimates  were  compared  with  flight  test  angle-of-attack 
calibrations. 


FIGURE  6 <  F-l 1 IR/TRCT  FU6ELR0E  CROSS  SECTIONAL  AREA  POINT  SELECTION 
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FUtCLROE  Ufjjjgw  E8TIMITER: 

Upwash  eatimates  for  the  F-lllA/TACT  fuselage  were  obtained  bf  an  equi- 
valent-body-of-revolution  analysis  and  the  Yaggy-Rogallo  production  software. 
In  order  to  define  the  equivalent  body  of  revolution,  a  graph  of  fuselage 
area  versus  f  oc=lage  station  was  obtained  from  the  physical  description  data 
contained  in  reference  17.  The  graph  which  is  shown  in  Figure  5  had  curves 
for  total  cross  sectional  area  as  well  as  fuselage  alone.  All  curves  except 
for  the  fuselage  alone  were  ignored,  and  that  curve  was  estimated  by  a  series 
of  points  chosen  along  the  curve  such  that  straight  lines  connecting  the 
points  accurately  represented  the  curve.  In  this  case  twenty-one  points  were 
adequate  to  accurately  define  the  curve  as  shown  in  Figure  6.  The  location 
of  the  center  of  pressure  of  the  angle-of-attack  vane  was  located  at  fuselage 
station  of  -68.45  inches.  As  previously  discussed,  the  radial  distance  and 
angular  location  were  not  easily  determined.  The  radial  distance  was 
set  to  the  lateral  location  of  7.875  inches  and  the  angular  location  was 
set  to  90.0  degrees  with  the  intent  of  adjusting  the  level  of  the  upwash 
estimation  curves  based  on  flight  test  data.  The  data  were  entered  on 
punched  cards  using  the  usr's  Guide  for  the  Yaggy-Rogallo  Production 
Software  as  a  reference  to  obtain  the  deck  shown  in  Figure  7.  The 
input  data  deck  was  processed  through  program  BODY  to  obtain  upwash 
estimates. 
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The  output  of  program  BODY  consisted  of  two  types  of  lineprinter  list* 
ings  describing  the  actual  geometric  data,  the  effective  geometric  data,  and 
the  upwash  estimates.  The  first  page  of  output  shown  in  Figure  8  printed 
the  input  and  calculated  actual  geometric  data.  The  output  is  described  in 
detail  in  Appendix  C.  The  remaining  output  from  the  input  data  cards  in  Fig¬ 
ure  7  were  pages  of  the  second  type  of  output.  One  page  was  printed  for  each 
Mach  number  where  an  upwash  estimate  was  requested.  Figure  9  gives  an  exam¬ 
ple  of  a  single  page  for  a  Mach  number  value  of  0.8.  Th«*  effective  geometry 
and  incremental  upwash  estimates  are  described  in  detail  in  Appendix  C  but 
the  important  value  is  the  total  upwash  estimate  in  terms  of  epsilon/alpha 
which  is  printed  at  the  bottom  of  the  page.  These  estimates  were  plotted 
v  versus  Mach  number  as  shown  in  Figure  10  to  give  a  concise  upwash  estimate 

for  the  fuselage. 
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HIMa  UPHMH  E8TjHRTC8: 

Upwash  estimates  for  the  F-111A/TACT  supercritical  wing  wars  obtained 
at  three  wing  sweep  angles  using  a  thin-airfoil  analysis  and  the  Yaggy- 
Rogallo  production  software.  Leading  edge  sweep  angles  of  26°,  35°,  and 
58°  were  selected  because  they  gave  a  wide  range  of  wing  sweep  angles  and 
because  the  best  flight  test  data  were  available.  The  geometry  of  the  wing 
was  obtained  from  physical  description  data  contained  in  reference  17.  The 
data  were  for  a  "theoretical  trapezoidal”  wing  so  that  it  was  not  necessary 
to  develop  an  equivalent  geometry.  The  physical  data  used  in  the  analysis 
are  presented  in  Table  1.  These  data  were  entered  on  punched  cards  using 
the  User's  Guide  for  the  Yaggy-Rogallo  Production  Software  as  a  reference  to 
obtain  the  input  data  deck,  shown  in  Figure  11  for  26°  wing  sweep,  and  two 
other  decks  not  shown.  The  input  data  decks  were  processed  through  program 
WING  to  obtain  upwash  estimates  for  the  three  wing  sweep  angles. 


TABLE  It  F-lllA/TRCT  MINO  ORTA 


HI 

NINO 

SPAN 

ASPECT 

RATIO 

-A-C/4 

PS  c/4 

MEM" 

684.02 

5.07 

BEK 

457.28 

Wtm 

617.34 

4.39 

ESS’ 

442.96 

B9 

471.64 

2.56 

m 

384.56 

The  output  of  program  WING  consisted  of  a  single  pagg  of  lineprinter 
listing  for  each  wing  sweep  angle.  The  output  for  the  26  wing  sweep  as 
shown  in  Figure  12  corresponds  to  the  input  data  deck  in  Figure  11.  The  out¬ 
put  which  is  described  fully  in  Appendix  C  were  obtained  for  the  three  wing 
sweep  angles  and  plotted  as  shown  in  Figure  13.  This  is  the  recommended  form 
for  using  the  data  in  data  reduction  and  analysis  and  normally  the  data 
should  be  left  in  this  form.  It  is  not,  however,  compatible  with  many  exist¬ 
ing  data  reduction  procedures  so  efforts  were  made  to  change  the  form  of  the 
data  and  combine  the  data  with  the  fuselage  estimates  to  obtain  a  single  up¬ 
wash  estimate. 
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FIOURE  lit  "NINO"  INPUT  ORTA  FOX  F-lllA/TACT  28°  MINO  8MEEF  CONFIGURATION 
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Combination  of  upwaah  estimates  for  the  F-111A/TACT  aircraft  was  easi¬ 
ly  accomplished  using  the  procedures  recommended  in  the  Yaggy-Rogallo  Analysis 
section  of  this  memorandum.  However,  the  c/CL  format  was  incompatible  with 
the  programs  and  data  base  which  were  used  to  compare  Yaggy-Rogallo  results 
with  flight  test  data.  The  estimates  were,  therefore,  converted  to  an  e/ a 
format  using  the  lift-curve  slopes  contained  in  reference  18  and  presented 
in  Figure  14  and  zero-lift  angles  of  attack  (not  presented) .  The  results 
were  satisfactory  for  an  angle-of-attack  calibration  comparison  and  are 
discussed  in  the  following  subsection.  At  th-  lower  lift  coefficients  (0.1 
to  0.2),  however,  a  50%  difference  in  estimated  upwash  angle  was  noted 
between  using  the  average  lift  curve  slope  obtained  at  moderate  lift 
coefficients  (near  0.465)  and  the  correct  one.  Although  changing  from 
c/C.  to  c/a  yielded  acceptable  results  for  this  study,  it  is  not  recommended 
for  normal  use. 

COHPWK I SOM  WITH  FLIOHT  TEST  OUTHt 

The  upwash  angle  estimates  generated  for  the  F-111A/TACT  aircraft  were 
ooaipared  with  the  actual  flight  test  data  obtained  during  the  angle-of-attack 
calibration.  The  indicated  angle-of-attack  values  were  obtained  from  the 
AFFTC  NBIU  on  the  aircraft  and  the  true  angle  of  attack  values  were  obtained 
from  the  inertial  navigation  and  pitot-static  systems.  The  data  ware  obtained 


YAGGY-ROGALLO  UPWASH  CALCULATION  PACE 

UPWASH  ESTlMATfc  FOR  LIFTING  SURFACE  COMPONENTS  1/  i 


E -1 1 IA/T  ACT  AIRCRAFT 
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LIFTING  SURFACE  GEOMETRIC  OtTAIL: 
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LIFTING  SURFACE  UPWASH  ESTIMATES: 
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0.0000  l.OOOC  1.5835  23.3400  .025482  .005026  .287976 

.1000  .9950  1.5911  73.4449  .025241  .004978  .285247 

.2000  .9798  1.6161  23.7684  .024515  .004835  .277036 

.3000  .9539  1.6599  24.3387  .023298  .004595  .263265 

.4000  .9165  1.7277  25.2111  .021561  .004257  .243860 

.5000  .8660  1.6284  26.4847  .019349  .003616  .218663 

.6000  . 80CC  1.9793  28.3411  .016584  .003271  .187416 

.7000  .7141  2.2173  31.1410  .013263  .002616  .1*9860 

.8000  .8000  2.6391  35. 7222  .009361  .001546  .105766 

.9000  .4359  3.6327  44.7098  .004864  .000963  .055191 

.4900  .1411  11.2250  71.8961  .000536  .000106  .006063 
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from  numerous  constant-altitude  accelerations  and  decelerations  accom¬ 
plished  throughout  the  flight  tests.  The  indicated  values  were  corrected 
to  remove  local  flow  effects  determined  from  a  wind  tunnel  test  documented 
in  Reference  19.  The  curve  of  local  flow  errors  presented  in  Figure  15  was 
used  to  change  the  indicated  angle-of -attack  values  to  "corrected  angle  of 
attack"  values.  The  upwash  angle  estimates  presented  in  Figures  10  and 
13  were  then  used  to  obtain  "calculated  angle  of  attack"  from  the  corrected 
values.  The  data  a.'e  presented  for  26,  35,  and  58  degree  wing  sweeps  in 
Figures  16(a),  16(b),  and  16(c),  respectively. 

The  plots  of  true  angle  of  attack  versus  indicated  angle  of  attack  show 
errors  in  indicated  angle  of  attack  which  are  normally  attributed  totally 
to  upwash.  However,  when  the  NBIU  calibration  was  applied  to  obtain  the 
corrected  angles  of  attack,  about  half  of  the  error  was  removed.  The  upwash 
angle  estimates  obtained  from  the  Yaggy-Rogallo  technique  were  applied  di¬ 
rectly  without  any  modification  to  obtain  the  calculated  angle  of  attack  val¬ 
ues.  As  can  be  noted  from  Figure  16,  these  estimates  were  excellent  in  pre¬ 
dicting  the  upwash  angles.  Close  examination  indicates  that  the  Yaggy-Rogallo 
technique  slightly  under-estimated  upwash  angle  for  the  26  degree  wing 
sweep  and  slightly  over-estimated  upwash  angle  for  the  58  degree  wing  sweep 
upwash,  but  the  errors  could  be  easily  removed  by  only  a  small  change  in 
level  of  the  upwash  estimation  curves  based  on  this  data.  On  the  58 
degree  wing  sweep  plot,  several  points  which  lay  outside  the  data  grouping 
were  removed  during  analysis.  The  points  were  not  random  scatter  but 
were  caused  by  the  inability  to  accurately  calculate  rate  of  climb  (thus 
flight  path  angle)  from  the  pitot  static  system  as  the  aircraft  acceler¬ 
ated  or  decelerated  through  a  Mach  number  of  1.0.  These  points  were 
actually  the  result  of  inability  to  accurately  calculate  true  angle  of 
attack  and  were  not  related  to  the  NBIU  or  the  Yaggy-Rogallo  technique. 
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2.  CURVE  OBTRINEO  FROM  REFERENCE  19. 
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NOSEBOOM  CORRECTION  ESTIMATES 


The  correction  of  the  NBIU  calibration  for  chaagea  in  nosebooa  upwash 
is  another  important  application  for  Yaggy-Rogallo  technique.  The  NBIU  cali¬ 
bration  contains  an  error  in  angle  of  attack  due  to  upwash  of  the  wind  tunnel 
noseboom  which  changes  when  a  different  noseboom  is  used.  Rather  than  sub¬ 
tract  the  whole  upwash  estimate  predicted  by  Yaggy-Rogallo  from  the  NBIU 
calibration  data,  the  difference  in  upwash  estimates  between  the  wind  tunnel 
boom  configuration  and  the  test  boom  configuration  is  added.  This  is  consis¬ 
tent  with  the  AFFTC  philosophy  of  making  only  incremental  corrections  to 
test  data  and  doesn't  tie  the  user  to  the  Y&ggy-Rogallo  technique  if  a  new 
and  better  technique  is  developed. 

The  wind  tunnel  noseboom  configuration  shown  in  Figure  17  was  processed 
through  the  Yaggy-Rogallo  production  software  as  was  the  KC-135/WINGLBTS 
noseboom  configuration  shown  in  Figure  18.  The  results  of  both  runs  are 
displayed  in  Figure  19.  The  difference  curve  was  entered  as  a  correction 
to  the  KC-1 3 5/WINGLETS  upwash  estimate  and  combined  with  estimates  for 
other  components. 


FI  (HIRE  17 1  NINO  TUNNEL  N08E900H  CONFIGURATION 


FIGURE  18s  KC-1SS/NIN0LET8  NOSEBOOM  CONFIGURATION 
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FIGURE  19i  KC-136/HINIH.ETS  CORRECTION  FOR  NOSEBOOK  CONFIOURRTION 


SUMMARY 


The  Yaggy-Rogallo  technique  for  estimating  upwash  angles  at  noseboom- 
mounted  vanes  has  been  documented  in  this  memorandum.  The  origin  of  the  tech- 
nique  has  been  investigated  and  referenced  to  establish  confidence  in  the 
technique.  The  limitations  and  assumptions  have  been  discussed  and  recom- 
mendations  on  use  of  the  technique  have  been  made.  The  source  of  all  curves 
and  equations  have  been  shown  and  an  audit  trail  of  all  data  has  been  esta¬ 
blished.  From  these  investigations,  standard  software  and  procedures  have 
been  developed  and  documented  for  future  use. 

The  use  of  the  Yaggy-Rogallo  technique  in  data  reduction  and  analysis 
significantly  increases  the  accuracy  ct  angle-of -attack  data  and  reduces 
the  work  in  analyzing  upwash  calibration  data.  By  predicting  characteristics 
and,  in  most  cases,  the  level  of  curves,  the  technique  allows  development 
of  data  reduction  curves  with  minimum  effort.  The  newly  developed  production 
software  produces  preliminary  curves  based  on  aircraft  geometry  alone. 

When  used  in  conjunction  with  an  NBIU  wind-tunnel  calibration  to  remove 
local-flow  effects  and  an  upwash  calibration  to  adjust  level  of  the 
upwash  estimation  curves,  the  technique  offers  a  uniform,  yet  flexible, 
ate t hod  of  obtaining  accurate  angle-of -attack  data. 
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THE  EQUATION  FOR  UPWASH 
BY  A  BODY  OF  REVOLUTION 
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INTRODUCTION 

The  basic  concepts  for  calculating  the  upwash  angle  induced  by  a  body 
of  revolution  were  devr'.oped  by  Theodor  Von  Karman  in  Reference  3.  Although 
the  goal  of  the  memorandum  was  calculation  of  pressure  distribution  around 
airship  hulls,  part  of  the  solution  involved  calculation  of  a  potential 
function  for  flow  around  the  body  of  revolution  and  hence  velocity  distribu¬ 
tion.  Paul  F.  Yaggy  applied  this  theory  to  determination  of  upwash  angle 
induced  by  nacelles  on  the  flow  at  the  propeller  plane  in  Reference  1.  In 
his  development  Yaggy  narrowed  the  development  to  determination  of  induced 
angle  only  and  validated  the  technique  by  comparison  with  test  data  on  sev¬ 
eral  nacelle  configurations. 

Yaggy 's  technique  was  later  applied  to  nacelles  on  swept  wing  aircraft 
by  Vernon  L.  Rogallo  in  Reference  2.  Although  Rogallo  corrected  other  parts 
of  his  analysis  for  compressibility  effects,  he  applied  Yaggy's  nacelle  equa¬ 
tions  without  correction  even  at  high  subsonic  Mach  numbers.  Violation  of 
potential  flow  assumptions  at  the  higher  Mach  numbers  lowered  confidence  in 
the  calculation.  Donald  R.  Bellman  developed  a  compressibility  correction 
to  the  Yaggy  equation  which  was  used  in  Reference  9.  The  correction,  which 
was  based  on  determining  an  "effective"  geometry  which  varied  with  Mach  num¬ 
ber,  allowed  acceptable  estimation  of  upwash  at  all  Mach  numbers.  Develop¬ 
ment  of  the  Yaggy  equation  and  the  subsequent  compressibility  correction 
are  outlined  in  the  following  sections. 

YROOY  EQUATION  DEVELOPMENT 

The  derivation  for  upwash  induced  by  a  nacelle  presented  in  Reference 
1  is  suitable  for  general  application  except  for  its  restriction  to  the  hor¬ 
izontal  meridian  plane.  The  following  derivation  of  the  Yaggy  equation 
closely  follows  Yaggy's  original  work  except  for  minor  terminology  and  sign 
changes  necessary  to  make  it  compatible  with  currently  accepted  aircraft 
coordinate  systems.  The  derivation  was  also  expanded  to  allow  varying  rad¬ 
ial  position  around  the  body. 

OEOWETR I C/C00R0 I NATE  SPECIFICATION : 

The  geometry  of  an  arbitrary  body  of  revolution  as  specified  ir  rectan*- 
gular  coordinates  is  shown  in  Figure  Al,  The  longitudinal  body  axis,  X  ax¬ 
is,  is  the  axis  of  revolution  (axis  of  symmetry)  and  is  positive  forward. 

The  lateral  body  axis,  Y  axis,  and  vertical  body  axis,  Z  axis,  are  positive 
as  shown.  The  body  is  at  angle  of  attack,  a,  such  that  a  transverse  flow 
in  the  direction  of  the  -Z  axis  of  velocity  W  is  induced.  The  angle  of 
attack,  transverse  velocity,  and  uniform  freestream  velocity,  V_,  are  rela¬ 
ted  by  equation  (Al) .  The  longitudinal  velocity  component,  U,  Is  defined  by 
U  -  V  cos  (a) .  The  lateral  velocity  component  is  considered  identically 
zero  for  all  cases  so  that  angles  of  sideslip  other  than  zero  are  never  con¬ 
sidered. 


a  =  sin"1 (W/VT) 


(Al) 
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FIGURE  Alt  BODY  OF  REVOLUTION  IN  RECTRNBULAR  COORDINATES 


The  potential  function  for  doublet  flow  is  more  simply  developed  in 
polar  coordinates  so  a  coordinate  transformation  was  performed.  The  angu¬ 
lar  position  around  the  longitudinal  axis,  ft,  and  the  radial  distance  from 
the  longitudinal  axis,  r,  of  an  arbitrary  point,  P,  are  defined  and  rela¬ 
ted  to  the  rectangular  coordinate  system  in  Figure  A2.  Distances  along  the 
longitudinal  axis  are  defined  by  a  new  variable,  d,  such  that  d  ■  0  where 
the  perpendicular  piane  containing  point  P  intersects  the  longitudinal  axis. 
The  sign  of  distances  are  determined  to  be  positive  when  the  point  P  is 
downstream  of  a  doublet  element  located  on  an  element,  dx,  of  the  longitud¬ 
inal  axis.  Thus,  distances  forward  of  d  -  0  are  positive  and  distances  aft 
of  d  -  0  are  negative.  The  definition  of  adjusted  longitudinal  distances 
and  of  the  remaining  angular  coordinate,  8,  of  the  polar  coordinate  system 
are  shown  in  Figure  A3.  The  angle  0  as  defined  in  equation  (A2)  goes  from 
6  =•  0  an  infinite  distance  upstream  (forward'  to  8  -  it  an  infinite  distance 
downstream  (aft) . 

0  =  cot'1  (d/I  r|  )  0  <  0  <  it  (A2) 


180* 


FIOURE  A2»  DEFINITION  OF  A  ANO  r 
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POTENT I RL  FUNCTItW  OCVCLOfWiHT  : 

An  infinitely  long  cylinder  of  constant  radius  in  a  transverse  flow 
nay  be  represented  by  covering  the  longitudinal  axis  with  doublets  of  moment 
per  unit  length,  p.  If  the  transverse  flow  is  along  the  -Z  axis  and  the 
doublets  are  oriented  with  their  axes  along  the  Z  axis,  the  potential  func¬ 
tion  of  a  doublet  element  of  strength  pdx  in  polar  coordinates  is  defined 
in  equation  (A3) .  The  potential  function  for  the  infinite  cylinder  can  then 
be  obtained  by  integrating  equation  (A3)  from  infinitely  upstream  /here  0 
>  0  to  infinitely  downstream  where  0  -  it.  The  potential  function  -s  then 
defined  by  equation  (A4).  The  potential  function  for  a  body  of  finite 
length  may  be  obtained  by  integrating  from  0  at  the  leading  edge,,  eLE»  to 
S  at  the  trailing  edge,  8^,  of  the  body  as  in  equation  (A5). 


cos!)  sin6  dO 

4irr 


(A3) 


♦ 


sinO  dO 


p  sinC 


dO 


(A4) 


(A5) 


Potential  functions  for  bodies  of  varying  radius,  R,  can  be  developed 
in  a  similar  manner  by  varying  the  strength  of  the  doublet  elements  along 
the  longitudinal  axis.  The  strength  of  the  elements  must  be  related  to 
the  radius  of  the  body  and  the  transverse  velocity  by  equation  (A6).  Sub¬ 
stituting  equation  (A6)  into  equation  (A5)  yields  the  potential  function 
for  a  finite  body  of  varying  radius  shown  in  equation  (A7) . 
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Substituting  equations  (A10)  and  (All)  into  equation  (A12)  and  rearranging 
yields  equation  (A13).  But  W  -  V  sin  a,  or  for  small  angles  of  attack  w 
-  V  a,  so  that  equation  (A13)  becomes  equation  (A14) .  However,  since  the 
induced  anqle  is  small,  c  =  vz/vt  and  ®flua*''ion  (*14)  becomes  (A15)  . 
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2  2  f  ®TE 

v„  =  w  /  r2  3ine  de 

J  0LE 


(A13) 


V.  -  v_o  /  9t!V  sine  ae 

^  '  »,E 


(A14) 


|  -  ainWa  f  \>  sl„8  d8 

2r^  J  9 

LE 


(A15) 


CWTHMIBIUTY  COWItCTI— : 

Violation  of  potential  flow  assumptions  used  in  formulating  tha  Yaggy 
equation  lowered  confidence  in  the  results  obtained  at  high  subsonic  Mach 
numbers.  Donald  R.  Bellman  developed  a  compressibility  correction  based 
on  determining  an  "effective"  geometry  which  varied  with  Mach  number.  To 
account  for  compressibility,  these  actual  geometric  distances  must  be 
adjusted  by  a  compressibility  factor,  B,  defined  by  B  *  /i-M* .  The  effect 
of  compressibility  is  to  create  an  effective  geometry  such  that  sections 
forward  of  the  vane  appear  closer  to  the  vane  and  sections  aft  of  the  vane 
appear  further  away.  This  continues  until  at  Mach  number  of  1.0  and  beyond 
the  body  has  practically  no  effect  and  upwash  is  considered  to  go  to  zero. 

For  Mach  numbers  between  0.0  and  1.0  the  effective  distances  are  defined 
by  equation  (A16) .  The  effective  angles  are  similarly  redefined  in  equation 
(A17) .  Equation  (A15)  then  is  transformed  to  equation  (A18)  when  corrections 
for  compressibility  are  included. 


d  =  d‘B 
d  =  d/B 


d  >  0  (for/ard) 
d  <  0  (aft) 


(A16) 


0  *  cot”* (d/| r | ) 


0  <  0  <  it 


(A17) 
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2^-C032»  / 0TFd2 
r  J  9 


sin9  d0 


(A18) 


moor  EQUATION  SOLUTION 


The  solution  of  the  Yaggy  aquation  was  accomplished  by  several  tech¬ 
niques  to  determine  the  best  approach  for  practical  application.  The  radius 
of  the  equivalent  body  of  revolution  as  a  function  of  6  is  rarely  an  analy¬ 
tic  function,  but  is  usually  a  graphical  or  tabular  representation.  The 
variation  of  the  radius  is  such  that  assumptions  about  the  relationship  of 
R  to  0  are  valid  only  for  small  ranges  of  0.  Equation  (Al8) ,  therefore,  was 
converted  to  a  summation  equation  so  that  the  integration  could  be  accomp¬ 
lished  over  small  segments  by  assuming  a  relationship  between  R  and  6.  The 
integral  portion  of  equation  (A18)  could  be  changed  to  a  sum  of  integrals 
as  shown  in  equation  (A19)  or  more  corcisely  in  equation  (A20) .  Thus,  if 
the  body  is  divided  into  N-l  segments  by  N  points  where  point  1  is  the  lead¬ 
ing  edge  and  point  N  is  the  trailing  edge,  equation  (A20)  can  be  substi¬ 
tuted  into  equation  (A18)  to  yield  equation  (A21) . 


R  sin0  d0 


sin0  d0 


sinO  d0 
(A19) 


f  0TE 


R  sinO  clO 


(A20) 


ft 

a 


sin2ft-cos2ft 


2r 


(A21) 


Further  simplification  of  equation  (A21)  depended  on  the  method  chosen 
to  perform  the  integration.  The  term  (sin*ft  -  cos’ft)  is  independent  of  the 
integration  so  ft  was  assumed  to  be  90  degrees  so  that  equation  (A21)  was 
transformed  to  equation  (A22).  Since  the  dominant  method  in  past  imple¬ 
mentations  had  been  numerical  integration,  several  numerical  algorithms  were 
tried  first. 


e 

a  = 


sin0 


d0 


(A22) 
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NUMERICAL  APPROXIMATIONS; 


The  initial  algorithm  investigated  for  solution  of  equation  (A22)  was 
one  that  had  been  used  extensively  in  past  application  of  the  Yaggy-Rogallo 
technique.  The  algorithm  is  presented  in  equation  (A23)  and  assumes  that 
the  radius  of  the  segment  is  constant  along  its  length  with  a  value  equal 
to  the  aft  radius  of  the  actual  segment.  During  a  sensitivity  study  which 
is  discussed  in  the  Sensitivity  portion  of  this  subsection,  the  algorithm 
was  found  to  have  a  high  sensitivity  to  the  number  of  segments  into  which 
the  body  is  divided.  Because  of  the  unacceptable  sensitivity,  three  other 
numerical  algorithms  were  tried. 


e 

a 


<v 


ei-l> 


(A23) 


The  three  numerical  methods  which  were  investigated  were  slight  varia¬ 
tions  of  the  same  algorithm.  In  each  case  the  segment  lengths  were  assumed 
to  be  short  enougn  that  R  and  0  could  be  assumed  constant  along  the  segment. 

Rear  Radius.  The  first  numerical  method  assumed  that  each  segment  of 
the  body  was  a  cylinder  with  constant  radius  equal  to  the  rear  radius  of  the 
segment  (R.)  and  angle  equal  to  the  average  value  for  the  rear  of  the  segment 
(0.)  and  ti*le  front  of  the  segment  0.  ,  so  that  equation  (A22)  was  transformed 
into  equation  (A24) . 


(A24) 


Forward  Radius.  Similarly,  the  second  numerical  method  assumed  that 
each  segment  of  the  body  was  a  cylinder  with  constant  radius  equal  to  the 
forward  radius  of  the  segment  { i ) •  The  angle  equal  to  the  average  value 
for  the  segment  was  again  used  on  this  method.  Equation  (A22)  was  trans¬ 
formed  into  equation  (A25)  using  the  second  parameters. 


c 

a 


sin 


O.+0.  , 

l  l-l 


(0i-0i-l> 


(A25) 


Average  Radius.  The  third  numerical  method  assumed  that  each  segment 
of  the  body  was  a  cylinder  of  constant  radius  equal  to  the  average  of  the 
front  and  rear  radii  of  the  segment.  Using  these  parameters,  equation  (A22) 
was  transformed  into  equation  (A26). 


(A26) 
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Equations  (A24) ,  (A25) ,  and  (A26)  were  used  to  evaluate  the  sensitivity 
of  numerical-integration  techniques  to  the  method  of  estimating  radius ,  to 
the  choice  of  number  of  points,  and  to  round-off  and  recording  errors  during 
hand  calculation. 

Sensitivity.  The  sensitivity  study  was  conducted  on  a  typical  fighter 
fuselage  by  picking  various  numbers  of  segments,  choosing  points  to  define 
the  segments,  and  calculating  upwash  estimates  using  equations  (A23) ,  (A24) , 
(A25) ,  ana  (A26) .  The  numbers  of  points  chosen  were  9,  21,  42,  and  83  points 
along  an  approximately  900  inch  fuselage.  For  the  9  and  21  point  tests,  the 
points  were  selected  to  best  represent  the  area  versus  fuselage  station  plot 
with  straight  line  segments.  The  42  point  test  essentially  defined  the 
curve  exactly  and  the  83  points  were  selected  by  linear  interpolation  be¬ 
tween  the  42  points.  The  sensitivity  of  the  primary  method  used  in  the  past 
to  the  number  of  points  selected  is  shown  in  Figure  A4  and  the  sensitivity 
to  the  number  of  points  and  to  the  method  of  estimating  radius  of  the  three 
numerical-integration  techniques  are  shown  in  Figure  A5.  The  number  of 
points  significantly  affected  the  accuracy  of  the  upwash  estimates  as  shown 
by  comparison  with  the  optimum  solution  obtained  from  exact  integration.  If 
either  the  forward  or  aft  radius  of  the  segment  were  chosen  as  in  Figure  A5, 
the  technique  did  not  converge  to  the  optimum  solution  with  83  points.  If 
the  average  radius  were  chosen  as  in  Figure  A5,  at  least  42  points  were 
needed  and  83  points  preferred  although  83  points  was  well  within  past  guide¬ 
lines  for  selecting  the  number  of  points.  The  sensitivity  of  the  techniques 
to  round-off  errors  incurred  during  hand  calculation  and  recording  were  in¬ 
significant. 


STHBOL  DESCRIPTION 

□  9  POINT  SOLUTION 
O  21  POINT  SOLUTION 
A  42  POINT  SOLUTION 
+  83  POINT  SOLUTION 


FIOURE  R4i  SENSITIVITY  OF  ACCEPTED  METHOD 
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FIGURE  RSt  COMPARISON  OF  SENSITIVITY  OF  THREE  NUHER I CAL- INTEORAT I ON 

METH008  OF  UPNA8H  ESTIMATION  FOR  A  TYPICAL  FIOHTER  FU8ELAOE 


The  sensitivity  study  of  the  numerical  techniques  indicated  deficiencies 
in  the  numerical  method;  the  prime  one  being  the  extreme  sensitivity  to  the 
number  of  points  selected.  A  user  could  unknowingly  select  too  few  points 
and,  coupled  with  a  poor  choice  of  radius,  get  highly  erroneous  results.  In¬ 
creasing  the  number  of  points  and  making  other  assumptions  about  variations 
of  R  with  6  could  improve  the  accuracy  of  the  methods  to  acceptable  levels, 
but  make  the  whole  technique  ungainly  and  difficult  to  implement.  Exact 
integration  of  the  integral  portion  of  equation  (A22)  proved  a  better  alter¬ 
native  than  improved  numerical  techniques. 

EjjjCT  INTEggtftTIOH  TECHMIOUCfc: 

Exact  integration  of  the  integral  portion  of  equation  (A22)  is  a  prac¬ 
tical  and  far  more  accurate  means  of  solving  the  equation.  As  with  the  nu¬ 
merical  methods,  assumptions  must  be  made  about  the  variation  of  R  with  6. 

It  is  not,  however,  necessary  to  assume  that  the  segments  approach  infini¬ 
tesimal  length  and  finite-length  segments  may  be  used.  The  length  of  the 
segments  may  be  as  long  as  desired  if  the  assumed  function  of  R  with  6  is 
valid.  Accuracy  of  the  solution  is  totally  dependent  on  proper  selection 
of  the  function  of  radius  along  the  length  of  the  body.  Three  functions 
of  R  with  0  were  investigated  to  determine  the  most  practical  method  for  fu¬ 
ture  application. 

Finite  Cylinder.  The  finite  cylinder  was  selected  as  the  first  estima¬ 
tion  of  body  shape  both  because  of  its  simplicity  and  because  of  its  corres¬ 
pondence  to  the  third  numerical  technique.  The  radius  of  the  segment  was 
taken  as  the  average  of  the  front  and  rear  radii  as  shown  in  equation  (A27). 
Using  this  relationship  which  makes  R  independent  of  6,  equation  (A22)  can 
be  transformed  to  yield  equation  (A28) .  The  integral  can  be  easily  evaluated 
to  yield  equation  (A29) . 


R 


VRi-i 


(A27) 


(A28) 


(A29) 


Finite  Cone.  The  finite  cone  solution  was  similarly  obtained  by  assum¬ 
ing  that  the  ratlius  of  each  segment  varies  linearly  along  the  segment.  The 
slope  of  the  cone  edge  is  defined  by  the  expression  (R,  ,  ~  Ri^^i-1  “  3.) 
so  that  the  radius  is  defined  by  aquation  (A30) .  Using  the  relationship  of 
3  to  B  defined  in  equation  (A17) ,  equation  (A30)  was  transformed  to  equation 
(A31) .  Equation  (A31)  was  rearranged  to  yield  equation  (A32) .  Equation  (A32) 
can  then  be  substituted  into  equation  (A22)  to  yield  equation  (A33). 
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Ri-rRi 


di-rdi 


(d-di)  +  R. 
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cot6i_l_cot0i 


(cotO-cotCL)  +  R^ 


(A31) 


R  =  KlCotO  +  K2 

where  «*  (R^_j-R^)/(cote^_  j,-cot0^) 

K2  “  Ri  "  Kicot0i 


(A32) 
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Linear  Area.  The  third  exact  integration  solution  was  obtained  by  assu¬ 
ming  that  the  cross-sectional  area  of  the  segments  varied  linearly  along  the 
length  of  the  segment.  The  area  is  then  defined  by  the  equation  (A35) .  The 
area,  however,  is  defined  by  the  equation  A  ■  *R2  and  the  relationship  of 
d  to  5  is  defined  in  equation  (Al7)  so  that  equation  (A35)  was  transformed 
into  equation  (A36) .  Equation  (A36)  was  then  transformed  into  equation  (A37) 
and  substituted  into  equation  (A22)  to  yield  equation  (A38) .  Equation  (A38) 
was  integrated  to  yield  equation  (A39) . 


A 


(d-d.) 


+  A. 
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(A3  5) 


Sensitivity.  The  sensitivity  of  the  exact  integration  techniques  to 
the  method  o£  estimating  radius,  to  the  choice  of  number  of  points,  and  to 
round-off  and  recording  errors  was  evaluated  with  the  same  data  points  used 
on  the  numerical-integration  techniques.  As  expected,  all  three  methods 
were  less  sensitive  to  the  number  of  points  selected  than  the  numerical  tech¬ 
niques  as  shown  in  Figure  A6.  The  finite  cone  showed  some  sensitivity,  but 
the  finite  cylinder  and  linear  area  methods  were  essentially  insensitive 
from  21  through  83  points  as  shown  in  Figure  A6.  Check  cases  were  hand  cal¬ 
culated  for  all  three  methods  and  the  finite  cone  proved  to  be  difficult  to 
hand  calculate  and  extremely  sensitive  to  round-off  and  recording  errors. 

No  particular  problems  were  noted  with  the  other  methods. 

The  analysis  of  both  the  numerical-  and  exact-integration  techniques 
indicated  that  the  most  practical  method  of  solving  the  Yaggy  equation  was 
the  linear-area-variation  method.  Only  on  bodies  which  were  actually  coni¬ 
cal  did  the  finite  cone  provide  better  results,  and  its  added  complexity 
and  sensitivity  were  unwarranted.  The  linear-area-variation  method  provided 
the  best  upwash  estimates  in  practically  all  cases,  was  the  least  sensitive 
to  the  number  of  points  selected,  and  was  almost  as  simple  and  adaptable  to 
hand  calculation  as  the  numerical  techniques.  Therefore,  equation  (A21) 
was  solved  using  the  linear-area-variation  method  to  evaluate  the  inte¬ 
gral  portion  and  resulted  in  equation  (A40) .  The  recommended  equation 
for  general  use  and  the  equation  implemented  in  the  production  software 
was  equation  (A40) . 
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APPENDIX  B 

DERIVATION  OF  THE  EQUATION  FOR  UPWASH  ANGLE 
INDUCED  BY  AN  ARBITRARY  WING 


INTRODUCTION 

The  estimation  of  induced  velocities  and  angles  in  the  flow  field  of 
arbitrary  wings  is  an  extensive  and  complicated  subject  which  received  much 
attention  by  early  theoretical  and  experimental  aerodynamists.  The  methods 
and  theories  developed  by  early  researchers  were  compared  extensively  with 
test  results  and  experimental  data.  Although  individual  methods  compared 
more  or  less  favorably  with  experimental  results,  general  conclusions  and 
trends  could  be  drawn.  The  most  important  conclusion  was  that  most  aero¬ 
dynamic  characteristics  of  a  vehicle  could  be  obtained  accurately  by  inde¬ 
pendent  treatment  of  the  wing,  fuselage,  and  nacelle  components.  Except 
for  flow  in  the  immediate  vicinity  of  the  junction  of  bodies,  the  aerodyna¬ 
mic  characteristics  could  be  determined  based  on  geometry  of  individual 
components  and  combined  with  similar  characteristics  from  the  other  compo¬ 
nents  with  only  simple  or  no  interference  corrections.  Another  important 
conclusion  was  that  the  effects  due  to  thickness,  twist,  camber,  and  flap 
deflection  could  be  handled  through  modification  of  spanwise  loading  so 
that  the  wing  could  be  assumed  to  be  a  "thin  airfoil".  These  conclusions 
allowed  treatment  of  wing  induced  velocities  and  angles  independent  of  other 
bodies  and  allowed  treatment  by  a  number  of  powerful  theoretical  techniques 
developed  for  thin  airfoils. 

Paul  F.  Yaggy  documented  in  Reference  1  a  technique  for  estimating  the 
upwash  angles  induced  by  wings  at  the  propeller  planes  of  aircraft  with  mul¬ 
tiple,  wing-mounted  engines.  Yaggy  used  basic  lifing-line  theory  to  esti¬ 
mate  the  wing  effects  for  straight  wings  in  incompressible  flow.  Using  a 
system  of  horseshoe  vortices  and  equations  from  H.  Glauert  in  Reference  15, 
Yaqqy  estimated  the  induced  angles  and  confirmed  by  comparison  with  experi¬ 
mental  data  that  the  technique  gave  good  results.  He  also  confirmed  that 
the  span  loading  of  the  wing  alone  could  be  used  in  the  technique  since  dif¬ 
ferences  in  span  loading  caused  by  nacelles  had  negligible  effect  on  upwash 
for  greatly  differing  plan  forms.  The  technique  worked  well  for  distances 
greater  than  60-percent  chord  ahead  of  the  leading-edge  of  the  wing  and  for 
unswept  wings  of  moderate  to  high  aspect  ratio.  However,  Yaggy  recognized 
deficiencies  in  the  technique  and  recommended  that  for  swept  wings  and  other 
configurations,  techniques  such  as  those  described  by  John  De  Young  in  Ref¬ 
erence  5  should  be  used. 

Vernon  L.  Rogallo  expanded  Yaggy* s  work  to  swept  wings  and  high  subson¬ 
ic  Mach  numbers  in  Reference  2.  Based  on  Yaggy' s  recommendation,  Rogallo 
incorporated  De  Young's  technique  for  estimating  wing-induced  upwash  angles. 
In  Reference  6  De  Young  investigated  three  possible  replacements  for  basic 
lifting-line  theory;  two  were  modified  lifting-line  techniques  and  the  third 
a  lifting  surface  technique.  As  expected,  the  Falkner  lifting  surface  tech¬ 
nique  was  most  accurate  but  most  difficult  and  time  consuming  for  making 
estimates  of  wing-induced  upwash  angles.  Of  the  tv o  modified  lifting-line 
techniques,  the  Weissinger  technique  proved  superior  to  the  Mutterperl  tech¬ 
nique  and  gave  only  sliqhtly  lesB  accurate  results  than  the  lifting  surface 
technique.  The  Weissinger  technique  was  considerably  less  time  consuming, 
easier,  and  proved  to  be  the  method  best  suited  to  an  overall  study  of  ef¬ 
fects  of  geometry  and  wing  loading  on  wing-induced  upwash  angle.  The  com- 
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parison  of  the  three  methods  was  done  by  comparing  the  three  methods  with 
experimental  data  of  aerodynamic  center,  lift  curve  slope,  and  other  aero¬ 
dynamic  parameters.  Because  the  techniques  did  accurately  predict  other 
aerodynamic  characteristics,  equations  developed  by  De  Young  in  Reference 
4  for  induced  downwash  were  considered  adequate.  Rogallo  adapted  the  equa¬ 
tions  to  upwash  prediction  by  changing  the  influence  coefficients  and  applied 
the  equations  to  prediction  of  upwash  angle  at  the  propeller  planes  of  wing- 
mounted  engines.  He  compared  the  results  with  experimental  data  at  Mach  num¬ 
bers  of  0.31  and  0.90  and  at  wing  sweeps  of  0.0  degrees  and  40.0  degrees  with 
good  results. 

The  equations  developed  by  De  Young  from  Weissinger  and  reiterated  by 
Rogallo  in  Reference  2  were  not  well  suited  for  flight  test  application.  The 
estimates  from  these  equations  were  confined  to  distances  close  ahead  of  the 
wing  where  propeller  planes  would  normally  lie.  Also,  calculation  of  upwash 
from  the  known  influence  coefficients  required  inputting  estimates  of  span 
loading  which  often  were  not  readily  available.  Investigation  by  Don  Bellman 
of  NASA/DFRC  led  to  a  presentation  of  the  data  in  a  concise  form  which  was 
easily  used  in  making  upwash  estimates  for  flight  test  application. 

MEISSINOER  METHOD 

The  Weissinger  method  of  prediction  of  wing-induced  upwash  is  a  modified 
lifting-line  theory  considered  by  De  Young  and  Rogallo  to  be  the  best  tech¬ 
nique.  The  method  uses  a  lifting-line  located  at  the  quarter-chord  line  and 
continuous  sheet  of  trailing  vortices  as  shown  in  Figure  Bl.  The  strength 
and  variation  of  the  vortex  filament  is  specified  by  a  boundary  condition  at 
a  finite  number  of  control  points  alonq  the  three-quarter-chord  line.  Flow 
perpendicular  to  the  surface  of  the  "thin  airfoil"  is  not  allowed  at  these 
control  points.  Weissinqer  investigated  the  use  of  7,  15,  and  31  control 
points  to  accurately  represent  the  wing.  Although  there  were  small  gains  in 
accuracy  with  the  greater  number  of  points,  the  gains  were  negligible  and 
seven  points  were  used  in  all  referenced  analyses.  The  spanwise  location  of 
the  seven  points  was  specified  in  terms  of  the  dimensionless  lateral  coordi- 
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nate,  n ,  which  was  defined  by  equation  (Bl)  where  y  is  the  lateral  coordinate 
and  b  is  the  wing  span.  The  spanwise  location  of  the  nth  control  point  is 
defined  by  equation  (B2) .  The  location  of  the  seven  points  is  shown  in  Fig¬ 
ure  B2. 

n  “  TB/TT  (B1) 

n  ■  cos  (n"/8)  (B2 ) 

n 

The  equations  for  calculation  of  induced  upwash  which  Rogallo  adapted 
from  De  Young's  development  in  Reference  4  were  the  result  of  Weissinger 
analysis.  The  upwash  can  be  calculated  at  points  with  specific  spanwise 
coordinates,  n  ,  corresponding  to  Weissinger' s  control  points  as  defined  in 
equation  ( B 2 )  ind  ahead  of  the  leading  edge  of  the  wing.  To  simplify  the 
analysis,  De  Young  assumed  symmetric  loading  so  that  the  analysis  could  ba 
confined  to  four  points  (v  ■  1,  2,  3,  4)  and  intermediate  results,  if  neces¬ 
sary,  could  be  determined  from  interpolation  as  discussed  in  Reference  4. 

The  equation  for  induced  velocity,  W,  over  freestream  velocity,  VT,  is  equa¬ 
tion  (B3) . 


<w/vT> 


I 

n» 


a  0 

vn  n 


(B3) 
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The  value  of  V/VT  is  a  function  of  influence  coefficients,  a  ,  which 
control  the  inf luence 1  that  the  dimensionless  circulation,  G  ,  at  ^fle  nth 
semispan  station  has  on  the  spanwise  station  v  where  upwashnestimates  are 
required.  The  dimensionless  circulation  is  defined  by  G  =  T  /bV„  where  the 
circulation  at  spanwise  station  n,  T  ,  is  dependent  on  Mich  nEtmber  and  geo¬ 
metry.  Equation  (B3)  is  awkward  to  work  with  frr  quick,  easy  estimates  of 
upwash.  The  equation  was  rearranged  and  rewritten  in  Reference  2  to  put 
upwash  in  terms  of  unit  lift  coefficient  and  other  more  easily  understood 
terms.  Equation  <B 4 )  is  the  rearranged  equation. 


where 

M  =  Mach  number 

B  —  /l— M* 

Cj  =  total  lift  coefficient 

k  *  average  of  ratio  of  actual  sectional  lift- 

curve  slopes  to  theoretical  values 

S  =  total  wing  area 

=  section  lift  coefficient 

c  =  section  chord  lenqth 

c  mean  aerodynamic  chord  length 

The  parameter  [C  c/C,cj  is  the  span  loading  coefficient  which  is  pre¬ 
sented  in  references  |  ana  7  for  numerous,  varied  configurations.  Results 
in  these  references  indicate  that  the  span  loading  coefficient  is  not  a  func¬ 
tion  of  Mach  number,  but  is  only  a  function  of  effective  geometry.  The  in¬ 
fluence  coefficients  are  functions  of  effective  wing  sweep  defined  by  equa¬ 
tion  (B5)  and  dimensionless  distance  ahead  of  the  quarter-chord  line  defined 
by  equation  (B6).  The  influence  coefficients  used  throughout  this  analysis 
are  given  in  Figure  17,  p.  35-46  of  Reference  2. 


(B5) 


distance  ahead  of  C/4 
(b/2 ) 


(B6) 
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BELLMAN  PARAMETERIZATION 


Don  Bellman  of  NASA/DFRC  simplified  application  of  the  Weissinger  tech¬ 
nique  in  making  upwash  estimates  for  flight  test  use.  Through  the  U3e  of 
normalizing  parameters,  the  Weissinger  technique  was  adapted  to  general 
flight  test  application  and  the  range  of  t/B  values  extended  to  values  which 
would  include  nosebooms.  Initially  equation  (B4)  was  simplified  and  conver¬ 
ted  to  parameters  more  generally  used  flight  test  personnel.  Since  the 
upwash  angle,  e,  is  defined  by  z  =  tan  (W/V  )  =  (W/V„)  and  aspect  ratio,  AR, 
is  defined  by  AR  =  b*/s,  equation  (B4)  was  rewritten  as  equation  (B7) .  The 


value  of  v  can  be  set 
line  (midspan)  of  the 
unnecessary. 


to  4 
wing 


CAR 

CT 


by  assuming  that  the  noseboom  is  on  the  center- 
(n  =  0)  or  close  enough  that  interpolation  is 


(B7) 


n=l 


Research  in  many  of  the  references  indicated  that  span  loading  coeffi¬ 
cients  varied  very  little  directly  with  Mach  number  but  had  major  variation 
with  A.  with  minor  variations  due  to  taper  ratio,  A.  Further  analysis 
showed^that  the  minor  variations  of  span  loading  coefficient  due  to  taper 
ratio  had  negligible  effect  on  upwash  estimates.  Values  of  influence  coef¬ 
ficient  are  a  function  only  of  t/B  and  A„  as  shown  in  Reference  2.  From 
this  functional  analysis  it  was  concluded  that  the  parameter  eAR/C  could 
be  plotted  versus  t/B  with  lines  of  A„.  Bellman  generated  a  curve  which  was 
disseminated  to  numerous  projects  and^used  successfully  both  at  NASA/DFRC 
and  AFFTC.  Several  curves  derived  from  the  original  were  obtained  from  diff¬ 
erent  sources.  Because  the  curves  from  different  sources  differed  slightly 
in  detail  and  there  were  questions  concerning  the  range  of  independent  para¬ 
meters  over  which  the  curves  were  valid,  the  original  curve  could  not  be 
determined.  Since  the  original  analysis  was  no  longer  available,  independ¬ 
ent  curves  were  developed  which  could  be  substantiated  with  data  and  docu¬ 
mented. 

DATA  ANALYSIS  AND  FAIR (NO : 

Redevelopment  of  the  Bellman  curve  was  done  with  theoretical  and  exper¬ 
imental  data  from  three  references  and  included  twenty-two  span  loadings  sub¬ 
stantiated  by  numerous  tests  and  studies.  Data  with  seven  geometric  wing 
sweeps  were  used  and  in  every  case  possible,  the  data  were  used  at  more  than 
one  effective  wing  sweep  to  insure  that  Ag  was  normalizing  Mach  effects  on 
span  loading.  Taper  ratios  from  0.0  to  1.0  were  used  to  validate  the  assump¬ 
tion  that  data  with  all  taper  ratios  could  be  combined.  Aspect  ratios  from 
2.0  to  10.0  were  used  to  insure  that  cAR/C  normalized  aspect  ratios  effect¬ 
ively.  A  summary  of  data  included  in  the  analysis  is  given  in  Table  Bl. 

Values  of  eAR/CL  were  calculated  for  effective  wing  sweep  values  0.0, 
31.0,  45.0,  and  60.0Ldeqrees  and  for  t/B  values  from  0.15  to  0.70  using  equa¬ 
tion  (B7 ) .  The  data  were  plotted  and  faired  as  shown  in  Figure  o3.  The 
data  appear  nonlinear  at  lower  t/B  values,  but  become  linear  above  t/B  of 
0.4  when  plotted  on  a  log-log  plot.  Some  minor  disagreement  between  the 
data  and  fairings  can  be  seen  at  t/3  of  0.7  and  A„  of  45.0  and  60.0.  These 
discrepancies  were  attributed  more  to  inability  to  accurately  read  the  in¬ 
fluence  coefficient  curves  than  to  trends  in  the  data.  As  stated  in  the 
note  in  Figure  B3,  there  is  some  ordered  effect  within  A„  groupings  attrib¬ 
utable  to  taper  ratio.  Generally  lower  taper  ratios  are^at  the  top  of  the 
grouping  and  higher  values  at  the  bottom  of  the  grouping,  but  the  effects 
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TABLE  BU  SUMMARY  OF  INCLUOEO  OATA 


POINT 

NUMBER 

-A-C/4 

TAPER 

RATIO 

ASPECT 

RATIO 

SPAN  LOAOINO  COEFFICIENT 
REFERENCE  AN0  FI0URE 

1 

0.0 

0.0 

0.0 

NACA  TN  1491.  FI0  10(C) 

2 

0.0 

0.0 

— 

ELLIPTIC  LOAOINO 

3 

o.o 

0.0 

0.4 

NACA  TN  2795.  FI0  7(B) 

4 

0.0 

0.0 

0.5 

NACA  TN  1491.  FI0  10(0 

5 

o.o 

0.0 

1 .0 

NACA  TN  1491.  FI0  10(0 

6 

31.0 

31.0 

0.4 

NACA  TN  1476.  FI0  2(0) 

7 

45.0 

31-0 

0.4 

NACA  TN  1476.  FI0  2(0) 

8 

45.0 

35.0 

0.5 

6.0 

NACA  TN  2795.  FIO  3(A) 

9 

45.0 

40.0 

0.4 

10.0 

NACA  TN  2795.  FIO  7(A) 

10 

45.0 

45.0 

o.o 

3.0 

NACA  TN  1491.  FIO  10(B) 

11 

45.0 

45.0 

0.5 

3.0 

NACA  TN  1491.  FIO  10(B) 

12 

45.0 

45.0 

1.0 

3.0 

NACA  TN  1491.  FIO  10(B) 

13 

60.0 

31.0 

0.4 

4.7 

NACA  TN  1476.  FIO  2(0) 

14 

60.0 

35-0 

0.5 

6.0 

NACA  TN  2795.  FIO  3(A) 

15 

60.0 

40.0 

0.4 

10.0 

NACA  TN  2795.  FIO  7(A) 

16 

60.0 

45.0 

0.0 

3.0 

NACA  TN  1491.  FIO  10(B) 

17 

60.0 

45.0 

0-5 

3.0 

NACA  TN  1491.  FIO  10(B) 

18 

60.0 

45.0 

1 .0 

3.0 

NACA  TN  1491.  FIO  10(B) 

19 

60.0 

56.3 

0.0 

2.0 

NACA  TN  1491.  FIO  7 

20 

60.0 

60-0 

0.0 

3.0 

NACA  TN  1491.  FIO  10(A) 

21 

60.0 

60.0 

0.5 

3.0 

NACR  TN  1491.  FIO  10(A) 

22 

60.0 

60.0 

1.0 

3.0 

NACA  TN  1491.  FIO  10(A) 

were  not  well  defined  and  are  insignificant  for  almost  all  applications. 
Althouqh  there  were  perceivable  differences  between  the  data  fairings  in 
Figure  B3  and  Bellman  curves  from  other  sources,  the  differences  are  within 
limits  that  could  arise  from  interpretation  of  the  data  and  were  ignored. 

The  new  fairings  were  used  in  all  analyses  in  this  memorandum  and  in  develop¬ 
ment  of  the  production  software  presented  in  Appendices  D  and  E. 

CURVE  FIT: 

The  data  fairings  in  Figure  B3  were  fitted  with  an  equation  to  allevi¬ 
ate  curve  reading  in  future  applications.  The  equation  ( B 8 )  is  valid  for 
i/3  values  above  0.4  and  is  the  recommended  procedure  to  extrapolate  to 
t/P  values  higher  than  0.7.  For  t/3  lower  than  0.4,  the  data  were  faired 
at  the  four  wing  sweeps  where  data  were  available  and  interpolated  to  obtain 
values  at  intermediate  wing  sweeps.  The  results  are  given  in  Figure  B4. 

The  assumptions  in  the  technique  deteriorate  rapidly  as  t/3  goes  below  0.15 
and  no  confidence  is  possible  in  extrapolation  below  0.15. 


loq  (eAR/C.  )  =*  Aloq  (  t/3  )  +  BA Q  ♦  C 
L  p 


T/3  >  0.4 


where 

A  » 
B  *= 
C  = 


-1.48897301 

-0.008447868 

-1.099368684 


(B8) 
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PROCEDURAL  OUTLINE  FOR  THE  YAGGY-ROGALLO 
UPHASH  ESTIMATION  TECHNIQUE 
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The  outline  contained  in  this  appendix  is  a  concise  presentation  of 
the  procedures  necessary  to  generate  upwash  angle  estimates  usinq  the 
Yaggy-Rogallo  technique.  It  is  intended  as  a  quick  reference  and  checklist 
for  users  who  are  familiar  with  the  vaqqy-Rooallo  technique  and  dots  not 
replace  the  detailed  discussion  elsewhere  in  the  memorandum. 
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PROCEDURAL  OUTLINE  FOR  THE  YAOGY-ROOALLO 
UPWASI!  ESTIMATION  TECHNIQUE 

I.  Perform  Yaqqy-Roqallo  analysis 

A.  Define  geometric  confiquration  (s)  to  be  analyzed 

B.  Divide  aircraft  into  two  t'^oes  of  comoonents 

1.  Non-lifting  bodies 

2.  Lifting  surfaces 

3.  Other  (see  discussion  for  classifvino  comoonents) 

II.  Define  tipwash  Estimates  for  Individual  Comoonents 


A.  Cylindrical  Components 

1.  Designate  individual  comoonents 

a)  Each  component  with  non-concurrent  centerline 

b)  Each  component  where  a  separate  estimate  is  desired 

2.  Define  cross  sectional  area  distribution  of  each  component 

a)  Obtain  definition  of  cross  sectional  area  versus  fuselage 
station  for  aircraft 

i)  Plot 

ii)  mabulation 

iii)  Scale  drawinq 

b)  Eliminate  area  due  to  liftinq  surfaces 
or  canards 

c)  Generate  final  plot  of  cross  sectional 
station 

3.  Select  points  from  curve  and  assiqn  numbers 

a)  Points  selected  so  that  straight,  line  seqments  connecting 
them  adequately  define  curve 

i)  Each  data  ooint.  consists  of  two  values 

(1)  Fuselage  station 

(2)  oorresnondinn  cross  sectional  area 

ii)  At  least  two  ooints  must  be  selected 

(1000  points  is  maximum  number  for  use  of  production 
software) 

iii)  Accuracv  depends  on  how  w®ll  straiqht-line  seqments 
represent  curve;  not  on  soacinq  of  points 
h)  selected  ooints  are  nur^ered  from  1  to  N 


such  as  wings,  tails, 
area  versus  fuselage 


^  H 


is  a*  leadirn  e^ae  o*  bodv 


i)  Point 

ii)  Point  "’I"  is  at  trailing  edqe  of  body 
iii)  Points  designated  at  ith  point 
Define  location  of  angle-of-attack  van® 

a)  Longitudinal  location  of  anole-of-attack  vane 

i)  Designated  *>v  center  of  pressure  of  vane(s) 
ii)  Puselage  station  must  po  :Ln  same  r»,*r»nc'>  svstem 
as  rest  of  data 

b)  Angular  and  radial  location 

i)  nesinnated  bv  center  of  pressure  of  vane(s) 
ii)  Known  geometrv 

r  =  distance  from  centerline  to  center 


of  pressure 

annular  location  of  center  of  pressure  with  0C 


iii) 


being  vertically  above 
Unknown  georetrv 
r  »  distance  from  plan®  of 
pressure 
2  =  RO.O0 


centerline  (degrees) 


svmmetrv  to  center  of 
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Select  Mach  numbers  where  unwash  estimates  are  desired 

a)  Restrictions 

i)  Mach  numbers  must  be  greater  than  or  ecrual  to  0.0 

ii)  Mach  numbers  must  be  less  than  1.0 

b)  Recommended  range  and  spacing 

i)  0.0  to  1.0  in  steps  of  0.10  is  usuallv  accentable 

ii)  0.0  to  1.0  in  steps  of  0.025  is  more  than  adequate 

iii)  Range  and  snacinn  should  b°  adjusted  to  meet  project 
needs 

Select  "production  software"  or  "hand  calculation"  approach 
(Results  will  be  almost  identical  but  "hand  calculation"  is 
long,  tedious,  and  time  consuminq) 

a)  "Production  software"  approach  (consult  User’s  Cuide  for 
Yaggv-Rogallo  Production  Software) 

i)  Develop  HEAP1  card 

ii)  Develop  HEAD2  card 

iii)  Develop  VANES  card 

iv)  Develoo  MACHV  card(s) 

v)  Develon  AREAV  cards 

vi)  Run  BODY  proqram 

vii)  Skip  to  sten  II,  A,  13.  of  this  outline 

b)  Hand  calculation  aonroach 

i)  Establish  a  work  sheet  for  actual  geometric  calcula¬ 
tions 

ii)  Establish  a  work  sheet  for  effective  geometric  calcu¬ 
lations  and  upwash  estimates  for  each  Mach  number 
where  an.  unwash  estimate  is  desired 

iii)  Calculate  compressibility  factor,  8,  for  each  Mach 
number  and  nlac»  a*  ton  o *  work  sheet 
°  «  >a-M* 

Perform  coordinate  transformation  so  that  distances  fore  and 
aft  are  measured  from  the  c“nfer  of  pressure  of  the  angl*  of 
attack  vane  (distances  are  positive  forward) 

i)  Equation  uses  fuselage  stations  previously  determined 

d.  -  ES  -FS, 
i  v  i 

ii)  Record  values  on  actual  q^ometrv  worksheet 


8.  Calculate  eouivalent  radius  a*-  each  selected  fuselage  station 


i  -  1,N 


8.  Calculate  effective  distances  corresponding  to  each  geometric 
distance.  Must  be  repeated  for  each  Mach  value  and  placed  in 
individual  worksheets. 

r?i  «  d^*0  d  >  D  (forward) 


ai  -  V’ 


d  <  o  (a't) 


i  -  1,  N 


10.  Calculate  effective  angles  corresponding  to  each  geometric 

distance.  Must  be  repeated  for  each  Mach  value  and  placed  on 
individual  worksheets. 

n,  ■  cot-3 (d./| r! )  i  ■  1  ,H 


IS 


11.  Calculate  the  contribution  o*  each  segment,  A  (c/a) ,  to  upwash 
angle  estimates 


A(o/a)  =  [K1(ain8i-sin9i_1)  +  K2 (cos9i_1-cos9i) ] 


'  2  2 

where  Kx  =  (R  i_1"Ri)  /  (cot9i_1~cot9i) 


K~  =  -X, cot9 .  +  Rf 
2  111 


12.  Sum  contributions  of  all  segments  to  obtain  an  upwash  angle 
estimate  for  the  body  at  each  Mach  number. 

—  =  z  A  (c/a )  (for  each  Mach) 

a  2r  i=2 

13.  Plot  upwash  angle  estimates  as  c/a  versus  Mach  number  and 
check  for  apnropriate  characteristics 

B.  Lifting  Surface  Components 


1.  Designate  individual  comnonentn 

2.  Define  an  "eguivalent  nlanform"  for  each  lifting  surface  having 
straight  leading  and  trailing  edges  and  tins  narallal  to  free- 
stream  flow. 

a)  n-fine  geometrv  of  lifting  surface 

i)  Fuselao-  station  of  guart»r-chord  line  at  midsnan, 

ii)  SpS^  of  lifting  surface,  b 

iii)  Aspect  ratio  of  lifting  surface,  AR 

iv)  Sweep  of  the  quarter-chord  line, 

b)  Oefin®  fuselage  station  of  cent»r  of  nressure  of  angle- 
of-att.ack  vanes,  FS 

3.  Select  Mach  numbers  where  unwash  estimates  are  desired 

a)  Restrictions 

i)  Mach  numbers  must  b«  greater  t^an  O.o 

ii)  Mach  numbers  must  be  less  than  1.0 

b)  Recommended  ranqe  and  snacina 

i)  0.0  to  1.0  in  s«-ens  of  0.10  is  usuallv  acceptable 

ii)  o.n  to  1.0  in  steps  of  o.o?s  is  more  than  adequate 

iii)  Range  and  snacing  should  b»  adjusted  to  meet  project 
needs 

4.  Select  "production  software"  or  "hand  calculation"  annroach 

a)  "Production  software"  annroach  (consult  Pser's  Guide  for 
Yaggv-rogallo  Production  Softv'are) 

i)  Pevelon  HFAOl  card 

ii)  n-velon  HFAD2  card 

iii)  Pevelon  Wivgp  card 

iv)  Develon  MA.CH'f  card(s) 

v)  Run  WING  program 

vi)  Skin  to  step  ii,  b,  7. 

b)  "Hand  calculation"  annroach 

i)  Establish  a  worksheet  for  each  lifting  surface  to 
be  analv7.ed 


ii)  Calculate  compressibility  factor,  6,  for  each  Mach 
number  on  place  adjacent  to  the  Mach  value  on  work¬ 
sheet  _ 

8  =  /PBF 

5.  Calculate  dimensionless  parameters 

a)  Calculate  dimensionless  distance,  t,  of  the  center  of  pres¬ 
sure  ahead  of  quarter-chord  of  wing 

t  =  (FSv-FSc/4)/(b/2) 

b)  Calculate  effective  wing  sweep,  A0 

P 

Ag  =  tan-1  (tanA^/8) 

c)  Calculate  dimensionless  parameter,  t/B 

6.  Calculate  upwash  parameter,  eAR/C^ 

a)  t/8  >  0.4 

£ AR/C.  =  10tA1^<T/e>  +  BA8  +  C] 

where  A  =  -1.488973010 
B  =  -0.008447868 
C  «  -1.099368684 

b)  0.15  <  T/8  <  0.4 

Curve  lookup 

c)  t/8  <  0.15 

Assumptions  violated  (see  text) 

7.  Calculate  e/Cl 

t/CT  =  (eAR/CT)/AR 

U  J-J 

8.  Plot  upwash  angle  estimates  as  e/C.  versus  Mach  number  and 
check  for  appropriate  characteristics 

HI.  Combine  upwash  estimates  from  various  components 

A.  ylindrical  components 

1.  Account  for  incidence  angles  such  as  noseboom  incidence,  iNfi, 
and  wing  incidence,  i 

w 

2.  Combine  upwash  angle  estimates 

e  =  (c/Ql)NB<0t”iNB)  +  {e/Qt)F(a-iW) 

3.  If  interference  corrections  are  required,  see  text. 

B.  Lifting-surface  components 
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1. 


Develop  lift  coefficient  for  each  component 
i)  Divide  lift  among  component 
ii)  Develop  CT  for  each  component  based  on  lift  of  com¬ 
ponent,  area  of  component,  and  local  dynamic  pressure 

2.  Calculate  total  upwash  estimate 


i 

i 


3. 


e 


(e/CL>W 


+ 


(£/cl>c 


If  required,  change  to  e/a  format  (not  recommended) 


i 


£ 


( (e/CT )  *CT  a) 


((£/cl) 


CL  “° 
a 


> 


i)  e/C  is  a  function  of  Mach  number 
ii)  CL  Land  are  a  function  of  Mach  number  and  lift 
a 

coefficient 


IV.  Apply  Results 


A.  Use  for  preliminary  estimate 

B.  Adjust  level  of  curve  based  on  flight  test  data  to  obtain  final 
curves. 


1 


M 


APPENDIX  D 

USER'S  GUIDE  FOR  THE 
YAGGY-ROGALLO  PRODUCTION  SOFTWARE 


INTRODUCTION 

The  equation*  and  procedures  developed  for  apDlication  of  the 
Yaggy-Rogallo  upwash  estimation  technique  were  incorporated  in  two  pro¬ 
duction  computer  programs.  Equations  developed  in  Appendices  A  and  B 
were  implemented  in  a  manner  to  insure  maximum  capability  for  users  with 
a  minimum  requirement  of  time  and  effort.  The  language  used  was  American 
National  Standards  Institute  (ANSI)  FORTRAN  77,  commonly  called  FORTRAN 
Version  5  (or  FORTRAN  V) ,  as  described  in  Reference  20.  This  newest 
version  of  FORTRAN  was  u  ,ed  to  insure  the  longest  possible  useful  life 
for  the  software  but,  more  importantly,  FORTRAN  V  supposedly  has  greater 
commonality  among  computers  than  FORTRAN  iv.  Ev»ry  effort  was  made  to 
create  software  which  would  move  eaai.lv  between  computers.  Only  a  sin¬ 
gle  card,  the  file  declaration  or  "PROGRAM"  card,  of  each  proqram  is  non- 
ANSI.  No  suggested  control-card  setups  vr®re  included  in  the  Pser's 
Guide  since  the  programs  are  very  small  and  easily  adapted  to  numerous 
application  technioues. 

The  two  computer  programs  correspond  to  the  two  tvpes  of  components 
into  which  an  aircraft  would  he  divided  in  the  preliminary  stages  of  a 
Yaggy-Rogallo  analysis.  Program  ponv  i«?  intended  for  makinq  unwash  est¬ 
imates  for  components  which  are  designated  an  bodies  of  revolution  and 
program  WING  is  intended  for  makinq  uowash  estimates  for  components 
designated  as  liftinq  surfaces.  Althouqh  the  two  programs  are  completely 
independent,  there  is  great  similarity  in  their  int»rna)  development, 
input  and  output,  and  t^ere  is  a  large  deor““  of  commons] itv  in  the 
input  data  cards  which  reduces  the  amount  of  work  to  prepare  input 
cards.  Both  programs  extensively  check  the  input  data  to  determine 
any  inputs  which  would  cause  abnormal  termination  of  the  program,  and, 
to  the  maximum  extent  possible,  check,  for  errors  which  would  giv»  un¬ 
reasonable  results.  Messag-3  are  written  bv  both  programs  to  the  user 
if  unexpected  inputs  are  encountered  to  aid  in  locating  the  erroneous 
card(s).  The  following  sections  describe  the  input  data,  the  output, 
and  error  checking  for  both  programs  in  detail.  If  more  detail  is 
needed,  the  user  should  consult  the  Programmer's  Guide  or  program  list¬ 
ings  in  Appendix  F. 

PROORRN  BOOT  -  B0DY-0F-REV0LUTI0N  UPHR8H  E8TIKRT0R 

Program  BODY  makes  uowash  estimates  for  bodies  of  revolution  using 
equations  recommended  in  Appendix  A.  The  program  does  not  require  the 
user  to  develop  the  equivalent  hody-of -revolution,  angular  relationships, 
and  "effective"  geometry  required  by  the  final  equations,  Instead,  the 
input  is  in  more  familiar  terms  such  as  lonqitudinal  (fuselage)  station, 
longitudinal  (fuselage)  cross  sectional  area  and  Mach  number,  with  the 
program  calculating  other  needed  parameters.  The  proqram  can  calculate 
upwash  estimates  for  as  manv  bodies  as  the  user  desires  during  a  single 
run.  The  normal  output  is  a  comprehensive,  linenrinter  listing  of  the 
input  data,  the  calculated  geometric  data,  th®  effective  geometric  data, 
and  the  upwash  estimates.  Should  the  user  improperly  prepare  the  input 
data,  the  output  data  will  be  greatlv  abbreviated  and  consist  of  as  much 
information  as  possible  or.  the  location  and  nature  of  the  error.  The 
error  messages  are  output  by  error-checking  portions  of  code  and  depend 

n 


The  input,  output,  and  error  checking 


on  the  location  and  type  of  error, 
are  described  more  fully  below. 

BOOT  INPUT  CMOS: 

The  input  cards  for  program  body  consist  of  six  tvnes  of  cards  which 
are  arranged  into  data  sets  as  shown  in  Finure  01.  Each  card  of  the  input 
data  sets  contains  one  of  six  card  identifiers  entered  as  alphanumeric 
data  in  the  first  five  columns  of  the  card.  The  card  identifiers  are 

used  as  card  titles  in  discussing  the  cards  as  well  as  serving  as  flags 

in  the  error  checkinq  process.  The  remainder  of  the  card  contains  input 
data  to  be  used  in  calculating  unwash  estimates  or  labeling  output.  There 
are  limitations  on  the  number  and  range  of  some  parameters  which  must  be 
observed.  The  content  of  each  type  of  card  and  limitations  on  specific 
variables  will  be  discussed  in  detail  below,  but  the  critical  information 
is  summarized  for  quick  reference  in  Fiaure  D2. 

HE ADI  Card.  A  single  HFAD1  card  is  the  first  card  of  each  input 

data  sed  'The  card  identifier  is  the  alphanumeric  characters  "HEAbl"  in 

the  first  five  columns.  Columns  6  through  RO  contain  75  alphanumeric 
characters  of  user-supplied  heading.  This  heading  is  the  first  of  two 
lines  of  heading  printed  on  each  page  of  output  to  identifv  and  clarify 
the  Output. 

HEAD2  Card.  A  sinole  HFAD2  card  is  th“  second  card  or  each  inout 
data  set.  The  card  identifier  is  the  alphanumeric  characters  "HEAD2" 
in  the  first  five  columns.  Columns  6  throuah  BO  contain  7^  alphanumeric 
characters  of  user-supplied  h3adinq  comprisinn  th“  second  heading  line 
to  be  printed  on  each  paoe. 

VANES  Card.  A  single  VANES  card  is  the  third  card  of  each  inout 
data  set.  The  card  identifier  is  the  alphanumeric  charcters  ",TANES"  in 
the  first  five  columns.  The  information  on  tMs  card  specifies  the  loca¬ 
tion  of  the  center  of  pressure  of  the  anqle-of -attack  vane  relative  to 
the  centerline  of  the  bodv  of  revolution  and  the  longitudinal  (fuselage) 
station  reference  plane.  The  first  parameter  is  a  two-character  abbre¬ 
viation  for  the  units  of  length  in  which  other  parameters  will  be  input. 
Two  alphanumeric  characters  such  as  "IN"  for  inches  or  "cm"  for  centime¬ 
ters  are  entered  in  columns  9  and  lb  and  stored  in  variable  HNim.  The 
only  use  within  the  proqram  for  this  variable  is  as  a  column  heading  for 
all  columns  having  units  of  length;  no  units  checking  or  conversion  is 
done.  Nothing  adverse  will  hannen  within  the  program  if  the  abbreviation 
doesn't  agree  with  the  actual  units  or  even  if  the  field  is  blank.  It 
is  important,  however,  that  all  values  b3  entered  in  consistent  units  of 
length.  If  the  vane  radial  distance  is  ent°red  in  inches,  for  instance, 
the  longitudinal  stations  must  be  in  inches  and  th“  areas  must  he  in 
square  inches  to  obtain  valid  results.  Accurate  entrv  of  the  units  abbre¬ 
viation  will  act  as  a  reminder  of  this  fact  and  is  highlv  recommended. 

The  next  parameter  is  the  lonqitudinal  station  of  the  center  of  pressure 
of  the  anqle-of-attack  vane.  The  parameter  is  entered  in  columns  11 
through  20  in  FlO.b  format  and  stored  in  variable  PS'rANF.  The  next  two 
parameters  are  the  radial  distance  from  tKe  centerline  and  angular  loca¬ 
tion  around  the  centerline  of  the  center  of  pressure  of  the  angle-of- 
attack  vane  which  are  stored  in  variables  RVANE  and  ANVANF  respectively. 
The  radial  distance  is  measured  from  the  centerline  to  the  center  of 
pressure;  it  must  be  non-zero  and  should  be  entered  as  a  positive  value 
in  F10.0  format  in  columns  21  through  30.  The  anqular  location  of  the 
center  of  pressure  is  measured  positive  clockwise  when  facing  forward 
and  is  zero  when  vertically  above  the  centerline.  The  value  of  the  angle 
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NO 

CARD  OF 

TYPE  CARDS  VARIABLE  FORMAT  COLUMNS  DESCRIPTION 


HEADl 


VANES 


1  HEAD 


FSVANB 


RVANE 


AN  VANE 


1-  5  Card  identifier  -  Must  be  "HEADl" 


5  Card  identifier  -  Must  be  "ITEAD2" 


23Us£3 


Card  identifier  -  Must  be  "VANES" _ 

Blank 

Two  alphanumeric  characters  as  abbreviat¬ 
ions  for  units  of  innut  (eg;  IN,  CM,  FT) 
Longitudinal  station  of  center  of  pressure 

of  angle-of-attack  vane _ 

Radial  distance  from  body  centerline  to 
center  of  pressure  of  angle-of-attack  vane 
Annular  location  of  angle  of  attack  vane  - 
Must  be  in  degrees 


MACHV  1  HEAD 
to 

14 

(as  NMACH 
req'd) 


AREAV  2  HEAD 
to 
1000 

as  FS ( I) 
req ' d )  AREA  I 


Card  identifier  -  Must  be  "MACHV"  on  each 
card _ 

Blank 

Number  of  Mach  number  values  to  be  read 
from  all  MACHV  cards  (Note:  This  variable 

is  active  on  first  card  only) _ 

Seven  Mach  number  values  (Note:  Number  of 
MACHV  cards  must  match  exactly  the  number 
of  cards  to  read  in  "NMACH"  values  at  7/ 
card)  _ 


Card  Identifier  -  Must  be  "AREAV"  on  eac 

card  _ 

Blank  . 

A  single  longitudinal  station  value  ~ 

A  sinnle  station  area  corresponding  to  the 
FS  value  on  the  s.».me  card 


1  HEAD 


1-  5  Card  identifier  -  Must  be  "FND 
This  terminates  a  data  set)  _ 


(Note : 


F I  (HIRE  02  •  8UTOWKT  OF  ORTH  CMOS  FOU  PROMftH  BOOT 


must  be  entered  in  decrees  in  FlO.o  format  in  columns  31  through  40. 

MACHV  Cards.  One  or  more  MACHV  cards  follow  the  VANES  card  in  each 
input  data  set. Each  MACHV  card  has  the  alohanumeric  characters  "MACHV" 
in  the  first  five  columns  as  the  card  identifier.  The  first  MACHV  card 
only  has  variable  NMACH  in  12  format  in  columns  0  and  10.  NMACH  contains 
the  number  of  Mach  numbers  where  upwash  estimates  are  desired.  The  num¬ 
ber  of  Mach  numbers  is  limited  to  98  by  program  dimensions  and  efforts 
to  enter  more  values  will  abnormal lv  terminate  the  program.  Each  MACHV 
card  will  hold  up  to  seven  values  of  Mach  number  in  7F10.0  format  in 
columns  11  through  80.  Since  seven  values  per  field  for  fourteen  cards 
yields  ninety-eight  values,  the  program  expects  no  more  than  fourteen 
MACHV  cards.  The  program  exoects  all  Mach  number  values  to  be  entered 
in  the  first  NMACH  fields  and  only  enouqh  cards  to  contain  NMACH  fields 
to  be  in  a  data  set.  Thus,  if  seventeen  Mach  number  values  are  to  be 
read  in,  it  would  require  exactly  three  MACHV  card?  the  first  two  would 
contain  full  seven  fields  each  and  the  third  would  contain  three  Mach 
number  values  in  the  three  left-most  fields.  All  Mach  number  values 
must  be  greater  than  or  equal  to  0.0  and  less  than  or  eoual  to  1.0  as 
assumed  in  deriving  the  basic  equations. 

AREAV  Cards.  At  least  two  AREAV  cards  must  follow  the  MACHV  card(s) 
in  each  input  data  set.  Each  AREAV  card  has  the  alohanumeric  characters 
"AREAV"  in  the  first  five  columns  as  the  card  identifier.  The  AREAV 
cards  contain  the  longitudinal  station/body  cross  sectional  area  oairs 
necessary  to  describe  the  bodv.  Values  are  entered  with  one  pair  per 
card  to  allow  easy  insert  or  removal  if  it  is  found  that  more  or  fewer 
points  are  needed  to  describe  the  body  accurately.  The  lonqitudinal 
station  value  of  a  pair  is  entered  in  F10.0  format  in  columns  11  through 
20  and  stored  in  the  FS  variable  arrav.  The  corresnonding  body  cross 
sectional  area  of  a  pair  is  stored  in  F10.0  format  in  columns  21  through 
30  and  stored  in  the  same  element  of  the  AREA  variable  array.  The  arrays 
are  dimensioned  to  1000  in  the  program  and  exceeding  1000  AREAV  cards 
will  cause  abnormal  termination.  The  cards  must  be  ordered  within  the 
AREAV  cards  in  ascending  value  of  longitudinal  station.  Two  consecutive 
cards  with  equal  values  of  longitudinal  station  are  allowed  to  easily 
account  for  sudden  increases  in  fuselage  area  such  as  at  inlets.  In  this 
case  a  zero  value  will  be  assigned  to  the  increment  uov’ash  for  the  segment. 

END  Card.  A  single  ENP  card  must  follow  the  AREAV  cards  to  end 
reading  of  AREAV  cards  and  to  terminate  the  input  data  set.  The  card 
identifier  is  the  alohanumeric  characters  "FNn  "  in  the  first  five 
columns.  The  rest  of  the  card  is  blank. 

B60T  BUTfUTi 

The  output  of  proqram  BODY  consists  of  two  tvpes  of  lineprinter 
listings  as  shown  in  Figure  D3.  The  first  type  of  listing  shown  in 
Figure  D3(a)  lists  the  actual  geometric  data  for  the  bodv  as  input  by 
the  user  and  the  calculated  geometrv  for  the  eauivelent-body-of-revol- 
ution.  The  second  type  output  shown  in  Figure  03(b)  lists  the  effective 
geometric  data  and  upwash  estimates  calculated  for  a  specific  Mach 
number.  The  format  of  both  tvoes  of  pages  wer»  planned  such  that  the 
actual  printing  takes  only  7.5  inches  in  width  and  about  8.5  in  length 
with  40  data  ooints  oer  nage  and  6  lin»s  per  inch.  This  leaves  adequate 
margins,  when  trimmed  to  8.5  bv  11. n,  to  olace  the  paqes  in  notebooks 
and  reports.  Both  types  of  pages  are  repeated  as  reauired  to  accommodate 
the  number  of  points  up  to  1000.  Further  details  o *  both  pages  are  in- 
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YAGGY-ROGALLO  UPwASH  CALCULATION 
GEOMETRIC  DATA  FOR  CYLINDRICAL  COMPONENTS 


PAGE 
1/  1 


FUSELAGE 


F-11IA/TACT  AIRCRAFT 


ANCLE-OF-ATTACK  VANE  LOCATION: 

RAOIAL  DISTANCE  FROM 
ANGULAR  LOCATION  .  . . 

BODY  CtNItRLINfc 

LONGI TLDINAL  STATION 
COMPONENT  AREA  DISTRIBUTION: 

OF  VANES  ....... 

CROSS  StCTIGNAL 

LONGI T GDI NAL 

PT 

LONG  I T  UOINAL 

AREA  Gf  BUOY 

DISTANCE  FROM 

NO 

STATION 

AT  STATION 

VANE  TC  STATION 

(  IM 

(SO  IN) 

(  IN) 

1 

0.0000 

0. 0000 

-68. *500 

2 

25.0000 

2CO.OOOO 

-93.5500 

3 

125.0000 

2000.0000 

-193. *500 

* 

175.0000 

2650. COOO 

-2*3. *500 

5 

2CO.OOOO 

3150. COOO 

-268. *500 

6 

250.0000 

*500.0000 

-318. *500 

7 

275.0000 

*925.0000 

-3*3. *500 

8 

ACC. 0000 

6275. COOO 

-A68. *500 

9 

**0.0000 

68  75  »  COOO 

-508. *500 

10 

*60.0000 

7575.0000 

-528. *500 

11 

*75.0000 

7600. OOUO 

-5*3. *500 

12 

525.0000 

685C.0000 

-593. *500 

13 

6  CO . 0000 

6200. COOO 

-668. *500 

1* 

65C.0000 

5375.0000 

-713. *500 

15 

7 CO. 0000 

*725.0000 

-768. *500 

16 

75C.OOOO 

3950. 0000 

-818. *500 

17 

775.0000 

3*25.0000 

-8*3. *500 

18 

78I.OOOO 

1325.0000 

-8*9. *500 

19 

800.0000 

1025.0000 

-  868 . *500 

20 

850.0000 

650. COOO 

-918. *500 

21 

8  75.0000 

350.0000 

-9*3. *500 

R)  OCCHCTRIC  ORTR 

FIpt«£  01 t  PROGRAM  BODY  SAMPLE  OUTPUT 
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21  POINTS 


7.8750  IN 
90.G000  DEG 
-68. *500  IN 

EQUIVALENT 
RAOIUS  AT 
STATION 
(IN) 


0.0000 
7.9788 
25.2313 
29.0*3* 
31.6651 
37.8*70 
39.5939 
**.6922 
*6.7601 
*9.1039 
*9.16*9 
*6.6950 
**.*2*3 
*1.3632 
38.7816 
35. *588 
33.0183 
2C.5368 
18.0629 
I*. 38*1 
10.5550 


YACGY-kOGALlG  '  PWASH  CALCULATION 
UP  *  ASH  ESTIMATE  FOk  CYLINDRICAL  COMPONENTS 
MACH  NUMBLk  -  .8000 


PAGE 
1/  1 


F-lllA/TAoT  AIRCRAFT 

FUSELAGE  21  POINTS 


LONGITUDINAL 

EFFECT IVL 

EQUIVALENT 

INCREMENT 

PT 

OISIANCE  FROM 

DISTANCE  FROM 

RADIUS  AT 

IN  EPSILON 

Mj 

VANE  TU  STATION 

VANE  TO  iTATIUN 

THETA 

STATION 

OVER  ALPHA 

(INI 

(IN) 

(RAO) 

(IN) 

1 

-68.A500 

-1 1A.0833 

3.0727 

O.OOOC 

0. 

2 

-9  3. 9500 

-1SS. 7S00 

3.0911 

7.9768 

.2385E-03 

3 

-19 J.A500 

-322. A167 

3.1172 

25.2313 

. 197AE-02 

<i 

-2A  l.ASOO 

-AOS. 7500 

3.1222 

29. 0 A  3  A 

• t A51E-03 

5 

-2  6  M • A5  00 

-AA  7 . A 1 6 7 

3 • 1 2 AO 

11.6651 

•2A78E-03 

t 

-  3 1  d .  ASOu 

-SJO. 7500 

3.1260 

37.0A70 

. A332E-03 

7 

-  3 A  3 • ASOO 

-S72.A167 

3.1270 

39.5939 

•  186AE-0  3 

8 

-AtA. ASOO 

-7e 0.7500 

3.131S 

AA .6922 

•  6 172E-C3 

9 

-5C8. ASOO 

-dA  7 . A  lb  7 

3.1323 

A6.78C1 

• 1295E-03 

1C 

-52d.A500 

-ddi,.  7SOO 

3.132/ 

A9.1039 

• 5939E-0 A 

11 

-5A1.A6U0 

-70S. 7500 

3.1329 

A9.18A9 

. A237E-0A 

12 

-59  3.  ASOO 

-989.0033 

3. 1 336 

Ab. 69  50 

. 113AE-03 

13 

-660. ASOO 

-  1 11 A.0MJ3 

3.13AS 

AA.A2A3 

.1127E-03 

1  A 

-/Id. ASOO 

-1197. Alb/ 

3.13S0 

A1.3632 

• A998E-0 A 

IS 

- /6H.A5T0 

-1280.  7  SOO 

3.1 3SA 

38.7816 

. 3536E-0A 

16 

-81 3. ASOO 

-  1  30A.CS1J3 

3.13SO 

35 • AS  8  8 

. 2AS9E-0A 

17 

-8A3.A5UO 

-1ACS.7S00 

3.1 JbO 

33.0183 

•  9220E-05 

18 

-mA9. ASOO 

—  1  a l S . 7  SOO 

3.1360 

20.5368 

•1JA8E-05 

19 

-  8  6  H  .  AS  00 

- 1AA7. A167 

3.1 362 

18.0629 

. 2022E-05 

2C 

-9  Id. ASOO 

-1530.  ( SO  ) 

3 . 1 36 A 

1A.30A1 

•  3 J90t-C5 

21 

-  )<t  3 .  At  00 

-1 >72. A  16  7 

3.13b6 

10. 5550 

.8915E-06 

TOTAL  UPAASH  ESTIMATE  ( t P S I  LON/ AL PH A  I  » 

. A926E-02 

B)  EFFECTIVE  OATH  AND  UPHRSH  E8TIHRTE8 
F I DURE  OSt  PROORRM  BOOT  SAMPLE  OUTPUT  (CONCLUOED) 
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eluded  below 


Geometric  Data,  The  first  page(s)  output  by  BODY  for  each  input 
data  set  present  the  geometric  data  for  the  case  as  shown  in  Figure  D3(a). 
The  output  is  identified  by  the  two  lines  of  user-supplied  heading 
entered  on  the  HEABl  and  HEAD 2  cards  and  printed  between  the  dashed  lines 
directly  below  the  standard  headings.  The  next  information  is  the  loca¬ 
tion  of  the  angle-of-attack  vane  relative  to  the  centerline  of  the  enui- 
valent-body-of-revolution  as  read  from  the  VANES  card.  The  radial  dis¬ 
tance  from  the  body  centerline  and  the  longitudinal  station  of  the  vanes 
are  printed  and  labeled  with  the  dimensional  abbreviation  innut  in 
variable  UNIT.  It  is  reiterated  that  no  conversions  or  checks  are 
done  to  insure  consistency  of  units.  The  annular  location  is  printed 
and  labeled  in  its  mandatory  units  of  decrees. 

The  component  area  distribution  is  then  printed  along  with  prelim¬ 
inary  calculations  of  the  dimensions  for  the  eouivalent-body-of-revolution 
The  point  numbers  are  assigned  to  lonqitudinal  station  values  used  to 
define  the  actual  geometry  of  the  body  ir  a  consecutive  manner  as  the 
AREAV  cards  are  read.  This  is  the  means  ot  identifyino  the  points  and 
segments  of  the  geometric  data  listing  and  corresponding  points  on  the 
effective  data  listing.  The  longitudinal  stations  and  corresponding 
cross  sectional  area  of  the  bodv  at  the  station  are  printed  exactly  as 
read  from  the  AREAV  cards  in  the  innut  data  set-.  The  lonqitudinal  dis¬ 
tance  from  the  vane  to  the  appropriate  station  is  the  distance  used  to 
determine  upwash  and  in  calculated  by  the  algebraic  addition  of  the 
longitudinal  station  and  the  lonqitudinal  station  of  the  anqle-of-attack 
vane.  The  equivalent  radius  at  the  station  is  calculated  a*  the  radius 
of  a  circle  having  an  area  eciual  to  the  cross  s«c*-ional  area  of  the 
station.  All  values  on  the  geometric  data  pace  represent  purely  geometric 
relationships  which  do  not  varv  with  Mach  number.  The  geometric  data 
are  output  as  one  set  per  input  data  set. 

Effective  Geometry  and  Upwash  Estimates.  The  second  type  of  page(s) 
output  by  BODY  are  the? effective  q°ometry  anci  upwash  estimat»s  shown  in 
Figure  D3(b).  These,  values  are  Mach  number  d°oendent  and  there  is,  there¬ 
fore,  one  set  of  paqes  for  each  Mach  number  requested  on  the  f'ACUV  card(s) 
The  appropriate  Mach  number  is  printed  as  part  of  the  standard  heading 
on  each  page  directly  above,  the  two  lines  of  user-supplied  heading  from 
the  HEAD1  and  HFAD2  cards.  The  point  numbers  on  this  paqe  directly 
correspond  to  those  on  the  geometric  data.  The  lonqitudinal  distance 
from  the  angle  of-attack  vane  to  the  longitudinal  s*-a*ion  is  reprinted 
on  this  page  for  completeness.  Th“  actual  lonqitudinal  distances  are 
corrected  for  compressibility  which  decreases  distances  ahead  of  the 
vane  and  lengthens  distances  aft  of  the  vane  to  'Held  effective  geometric 
distances.  The  angles  between  th»  c°nterline  o*  the  body  and  lines 
connecting  the  center  of  the  eouivalent  circles  at  each  effective  station 
are  calculated  as  T1IFTA  and  printed  in  units  of  radians.  The  eouivalent 
radius  at  the  station,  which  is  reprinted  for  completeness,  and  'HIETA 
are  the  primary  innuts  to  determine  the  upwash  increment  for  the  seqment 
in  terms  of  upwash  angle,  eosilon,  over  anqle  of  attack,  aloha.  The 
upwash  increment  for  a  seqment  is  printed  at  the  lonqitudinal  station  for 
the  aft  station  of  the  segment.  On  the  last  page  of  listinq  for  each 
Mach  number  is  printed  the  sum  of  unwash  estimates  for  all  segments  to 
give  a  total  upwash  estimate  for  t*e  bodv  a*  the  given  ”ach  number  in 
terms  of  epsilon  over  aloha. 

BODY  ERROR  CHECK  I MO : 


Program  BODY  performs  numerous  error  checks  as  data  is  read  and  as 
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upwash  calculations  are  made.  The  main  purpose  of  the  error  checking 
is  to  avoid  abnormal  termination  of  the  proqram  or  to  give  the  user  as 
much  information  as  possible  on  the  location  and  nature  of  the  error  if 
termination  is  unavoidable.  Several  minor-error  checks  are  made  which 
change  internal  values  which  are  obviously  wrong  and  would  cease  execu¬ 
tion  of  the  program  immediately  if  uncorrected.  These  changes  which 
were  included  mainly  to  protect  the  user  from  gross  misapplication  or 
keypunch  errors  are  not  accompanied  by  error  messages.  Since  these 
checks  might  limit  the  experienced  user,  give  unexpected  results,  or 
cause  confusion  due  to  apparent  inconsistencies,  thev  are  discussed 
briefly.  The  main  body  of  the  checks  are  for  major  errors  with  input 
data  which  do  not  have  easily  determined  solutions  and  which  must  term¬ 
inate  program  execution.  These  major-error  checks  write  error  messages 
to  the  user  prior  to  termination  to  identify  where  and  why  termination 
occurred.  Both  the  minor-error  and  major-error  checks  are  discussed 
below. 


Minor-Error  Checking.  Errors  which  appear  obvious  and  usually  re¬ 
sult  from  keypunch  errors  or  gross  misapplication  of  the  technique  were 
corrected  without  accompanying  output  messages.  This  was  done  to  allow 
continued  execution  of  the  program  so  that  output  could  be  produced  to 
aid  error  analysis  and  so  that  subsequent  data  sets  could  be  processed. 

In  some  cases  this  results  in  invalid  or  unexpected  results  which  should 
be  investigated  and  fixed;  in  other  cases  it  simoly  allows  data  to  be 
entered  in  a  convenient  manner  when  it  would  normally  terminate  program 
execution. 

Mach  number  values  entered  on  MACHV  cards  were  checked  for  values 
less  than  0.0  and  greater  than  1.0.  Obviously  values  less  than  0.0  have 
no  meaning  and  usually  result  from  keypunch  errors.  Values  less  than 
0.0  were  changed  to  0.0  internally  but  printed  at  their  original  value 
to  allow  them  to  be  found  and  corrected.  Mach  number  values  greater  than 
1.0  are  equally  invalid  since  velocities  cannot  b«  induced  ahead  of  a 
body  in  supersonic  flow.  Upwash  estimates  for  Mach  numbers  greater  than 
1.0  were  set  to  zero  and  execution  continued. 

The  longitudinal  station  values  input  on  ARFAV  cards  to  define  the 
body  of  revolution  wore  checked  on  input  to  insure  there  were  no  descen¬ 
ding  data  values,  "‘his  was  necessary  to  insure  valid  upwash  estimates 
were  obtained.  Two  cards  with  the  same  longitudinal  station  were  allowed, 
however,  to  conveniently  represent  step  increases  in  cross  sectional 
area  such  as  at  inlets  without  resort  to  or.e  station  slightlv  less  and 
one  slightly  more  than  the  actual  station.  An  incremental  upwash  of 
0.0  was  assigned  to  the  "segment"  defined  by  the  two  equal  longitudinal 
stations. 

Major-Error  Checking,  rheckinq  for  major  errors  which  necessitated 
program  termination  was  accomplished  in  subroutine  RDDATA .  Any  major 
error  discovered  resulted  in  an  error  message  followed  bv  program  termin¬ 
ation. 


Normal  program  termination  resulted  from  reading  an  end-of-file 
designator  vrhile  attempting  to  read  a  HFAD1  card.  Unexpected  reading  of 
an  end-of-file  designator  by  any  other  r»ad  statement  results  in  the 
statement 

*  "  IS  IN  ERROR  -  END-OF-FILE  FNrotTNTFPFO 
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where  the  quotes  have  the  five  alphanumeric  characters  of  the  card  iden¬ 
tifier  being  sought  when  the  end-of-file  designator  was  encountered. 

The  input  data  set  should  be  reassembled  with  the  proper  number  and 
types  of  cards.  A  related  message  results  if  the  input  data  set  is 
improperly  arranged  or  a  card  identifier  is  misspelled.  The  statement 

"  "  IS  IN  ERROR  -  CARD  READS  " 

will  be  printed  if  an  unexpected  card  identifier  is  encountered.  The 
first  set  of  quotes  enclose  the  five  alphanumeric  characters  of  the  ex¬ 
pected  identifier  while  the  second  set  enclose  the  eiqhty  alphanumeric 
characters  on  the  card  actually  read.  Deck  must  be  checked  for  proper 
order  of  the  input  data  set. 

Other  major  errors  result  from  improper  number  or  arrangement  of  the 
AREAV  cards.  As  the  cards  are  read,  the  longitudinal  station  values  are 
checked  for  ascending  order.  If  a  descending  value  is  found,  the  state¬ 
ment 


LONGITUDINAL  STATION  IS  NOT  IN  ASCENDING  ORDER 

is  written  and  execution  terminated.  The  AREAV  cards  must  then  be  man¬ 
ually  checked  and  reordered.  If  reading  of  AREAV  cards  is  terminated 
normally  by  reading  an  END  card  and  only  one  ?REAV  card  has  been  read, 
the  statement 

NUMBER  OF  LONGITUDINAL  STATIONS  MUST  EXCEED  1 

is  written.  Since  two  cards  are  necessary  to  define  a  single  segment, 
at  least  one  additional  AREAV  card  must  be  added  to  adequately  define 
the  body  of  revolution. 

PROMAN  MIMA  -  LIFT  I  MO  SURFACE  UPMASH  ESTIMATOR 

Program  WING  makes  upwash  estimates  for  lifting  surfaces  using 
equations  developed  in  Appendix  B.  The  user  is  reauired  to  inout  the 
geometry  of  the  lifting  surface  to  be  analvzed  and  the  location  of  the 
center  of  pressure  of  the  angle-of-attack  vane.  In  most  cases  the  test 
aircraft  will  resemble  the  Weissinger  planform  discussed  in  Appendix  B 
closely  enough  that  geometric  data  can  be  obtained  directly  from  the 
physical  description  of  the  aircraft  and  input  into  the  program.  If, 
however,  the  test  aircraft  doesn't  match  the  assumed  planform  closely 
enough,  the  user  must  enter  the  geometry  of  an  "equivalent"  Weissinger 
planform  as  discussed  in  the  yaggy-Roqallo  Analysis  section  of  this  re¬ 
port.  The  program  calculates  all  "effective"  geometry  to  account  for 
compressibility  effects  and  calculates  unwash  estimates.  The  program 
can  calculate  upwash  estimates  for  as  many  lifting  surfaces  as  the  user 
desires  during  a  single  run.  '"he  normal  output  is  a  comprehensive,  line- 
printer  listinq  of  the  input  geometric  data,  the  effective  geometric  data, 
and  the  upwash  estimates.  Should  the  user  improperly  prepare  the  input 
data,  the  output  will  be  greatly  abbreviated  and  consist  of  as  much  in¬ 
formation  as  possible  on  the  location  and  nature  of  the  error.  The  error 
messages  are  output  by  error-checking  portions  of  cod®  and  depend  on 
the  location  and  type  of  error.  The  innut,  output,  and  error  checking 
are  described  more  fully  below. 


71 


* 

*$  I 

.■>  i 


Mjjlft  IjjfjlT  CUMt:  | 

The  input  cards  for  proqram  WING  consist  of  five  tynes  of  cards  f 

which  are  arranged  into  data  sets  as  shown  in  Fiqure  D4.  Each  card  of 

the  input  data  sets  contains  one  of  five  card  identifiers  entered  as  j  « 

alphanumeric  data  in  the  first  five  columns  of  the  card.  The  card  iden¬ 
tifiers  are  used  as  card  titles  in  discussinq  the  cards  as  well  as  ser¬ 
ving  as  flags  in  the  error  checking  process.  The  remainder  of  the  card 
contains  input  data  to  be  used  in  calculating  unwash  estimates  or  labeling 
output.  There  are  limitations  on  the  number  and  ranqe  of  some  parameters 
which  must  be  observed.  The  content  of  each  tvoe  of  card  and  limitations 
on  specific  variables  will  be  discussed  in  detail  below,  but  the  critical 
information  is  summarized  for  quick  reference  in  Fiqure  D5. 

HEAD1  Card.  A  single  HFAD1  card  is  the  first  card  of  each  input 

data  se€.  The  card  identifier  is  the  alphanumeric  characters  "HEADl"  in 

the  first  five  columns.  Columns  6  through  80  contain  78  alphanumeric 
characters  of  user-sunolied  heading.  This  headina  is  the  first  of  two 
lines  of  heading  Printed  on  each  naqe  of  output  to  identify  and  clarify 
the  output. 

HEAD2  Card.  A  single  HFAh2  card  is  the  second  card  of  each  input 
data  set.  fKe- card  identifier  is  the  alphanumeric  characters  "HEAD2"  in 
the  first  five  columns.  Columns  f  through  80  contain  78  alphanumeric 
characters  of  user-supplied  heading  comnrisinq  the  second  headinq  line 
to  be  printed  on  each  page. 

WINGS  Card.  A  sinqle  WINGS  card  is  the  third  card  of  each  input 

data  set.  The  card  identifier  is  the  alphanumeric  characters  "WINGS"  in 

the  first  five  columns,  mhe  information  on  this  card  specifies  the  loca¬ 
tion  of  the  center  of  pressure  of  the  anqle-of-attack  vane  relative  to 
the  location  of  the  lifting  surface  and  the  qeometrv  of  the  lifting  sur¬ 
face.  The  first  Parameter  is  a  two-character  abbreviation  for  the  units 
of  length  in  which  other  parameters  will  b»  inout.  Two  alphanumeric 
characters  such  as  "IN"  for  inches  or  "CM"  for  centimeters  are  entered 
in  columns  8  and  10  and  stored  in  variable  PNI"'.  The  only  use  within 
the  program  for  this  variable  is  as  a  column  headinq  for  all  columns 
having  units  of  length;  no  units  checkinq  or  conversion  is  done.  Nothing 
adverse  will  happen  within  the  program  if  the  abbreviation  doesn’t  agree 
with  the  actual  units  or  even  if  the  field  is  blank.  It  is  important, 
however,  that  all  values  be  entered  in  consistent  units  of  l*nath. 

If  the  longitudinal  station  of  the  center  of  pressure  of  the  angle-of- 
attack  vane  is  in  centimeters,  for  instance,  the  lonqitudinal  station 
of  the  quarter-chord  line  at  midspan  and  the  span  must  also  be  in  centi¬ 
meters  to  obtain  valid  results.  Accurate  entry  of  the  units  abbreviation 
will  act  as  a  reminder  of  this  fact  and  is  highly  recommended.  The  next 
parameter  is  the  longitudinal  station  of  t**  center  of  pressure  of  the 
angle-of-attack  vane.  The  parameter  is  entered  in  columns  11  throuqh  20 
in  F10.0  format  and  stored  in  variable  FSVANF.  ^he  n«xt  parameter  is 
the  longitudinal  station  of  the  quarter-chord  line  at  midsnan  and  is 
stored  in  variable  FSQGL.  FSOGL  is  entered  in  columns  21  throuqh  30  in 
F10.0  format.  The  last  three  parameters  which  give  the  planform  geometry 


71 


NO. 

CARD  OF 

TYPE  CARDS  VARIABLE  FORMAT  COLUMNS  DESCRIPTION 


UEADl  1  HEAD _ A5 _ 1-  5  Card  identifier  -  Must  be  "HEAP1" _ 

_ HEADl  A75  P-RQ  First  line  of  user  specified  heading 

HEAD2  1  III'.AD _ A5 _ 1-  5  Card  identifier  -  Must  be  "HFAD2" 

_ HEAD2  A75  6-RO  Second  line  of  user  specified  heading 


WINGS  1  HEAD  A5  1-  5  Card  identifier  -  Must  be  "WINGS" 

3X  6-8  Blank  ~~ 

OnTt  S2  9-15  Two  alphanumeric  characters  as  abbreviation 

for  units  of  irnut  (eg;  IN,  CM,  FT) _ 

FSVANF  fTO  11-20  Ijonqitudinal  station  of  center  of  pressure 

_ of  ancle-of-attack  vane _ _ 

FSOCL  FTO  5T^T5  Longitudinal  station  o<!  gurrter-ehord  line 

_ at  midspan _ 

SPAN  P10.O  31—4 0  Span  of  liftinq  surface 

AR  FlO.O  41-50  Asnect  ratio  of  lifting  surface 

WSQCL  Pitt .  0  51-65  Wing  sweep  of  quarter  chord  line 

27J3?  ST^17)  Manic 


MACHV  1  HEAD  A5  1-  5  Card  identifier  -  Must  be  "*V\CHV"  on  each 

to  _ card _ 

14  3X  6-  5  Blank 

(as  UMACH  T 2  o—l tt  'lumber  of  Mach  number  values  to  be  read 

req'd)  from  all  MACHV  cards  (Mote:  This  variable 

_ is  active  on  first  card  only)  _ 

XMACH  7Pl5. 5  11-R0  Seven  Mach  number  values  (Note :  Number  of 

MACHV  cards  must  match  exactly  the  number 
of  cards  to  read  in  "M*»ACH"  values  at  7/ 
card)  _ 


are  the  span  of  the  lifting  surface,  aspect  ratio  of  the  lifting  surface, 
and  sweep  of  the  quarter-chord  line.  The  parameters  are  stored  as  vari¬ 
ables  SPAN,  AR,  and  WSQCL  resDectively .  All  three  are  entered  in  ”10.  h 
format  in  columns  31  through  ^0,  41  through  50,  and  51  through  60  resnec- 
tively.  The  remainder  of  the  card  is  blank. 

MACH”  Cards.  One  or  more  MACHV  cards  follow  the  ”ANFF  card  in  each 
inout  data  set.  Fach  *IACH”  card  has  the  alphanumeric  characters  "MACHV" 
in  the  first  five  columns  as  the  card  identifier.  Th“  first  MACHV  card 
only  has  variable  NMACH  in  12  format  in  columns  o  and  10.  N^ACH  contains 
the  number  of  ”ach  numbers  where  uowash  estimates  are  desired.  The 
number  of  Mach  numbers  is  limited  to  OR  bv  program  dimensions  and  efforts 
to  enter  more  values  will  abnormally  terminate  th«  orooram.  Fach  ‘1ACH” 
card  will  hold  up  to  seven  values  of  *’AC!I  number  in  7vin.n  format  in 
columns  11  through  80.  Since  seven  values  ner  field  for  fourteen  cards 
yields  ninety-eight  values,  the  nroqram  expects  no  more  than  fourteen 
MACHV  cards.  The  program  expects  all  Mach  number  values  to  be  entered 
in  the  first  NMACH  fields  and  only  enough  cards  to  contain  NMACH  fields 
to  be  in  a  data  set.  Thus,  if  seventeen  Mach  number  values  are  to  be 
read  in,  it  would  require  exactly  three  MACH”  cards;  the  first  two  would 
contain  full  seven  fields  each  and  the  third  would  contain  three  Mach 
number  values  in  the  three  left-most,  fields.  Al  ’  Mach  number  values 
must  be  greater  than  or  eaual  to  0.0  and  less  than  or  eoual  to  l.o  as 
assumed  in  deriving  the  basic  aquations. 

FND  Card.  A  single  FND  card  must  follow  the  MACH”  cards  to  end 
reading  of  ARKAV  cards  and  to  terminate  the  innut  data  3»t.  "h®  card 

identifier  is  the  alphanumeric  characters  "FNn  "  in  tb®  first  five 
columns.  The  rest  of  the  card  is  blank. 

HI HP  QUTrUT: 

The  output  of  program  WIN”.  consists  of  a  sirale  page  of  line- 
printer  listing  as  shown  in  Flour®  D6.  mhe  outout  lists  the  actual 
geometric  data  for  the  lifting  surface  as  inout  hv  the  user  as  well 
as  the  calculated  and  effective  geometry  at  each  Mach  numb®r  renuested 
and  unwa3h  estimates  for  each  specific  "ach  nuirt®r.  Th*  format  of 
th®  listing  was  Planned  such  tha*-  the  actual  printing  takes  only 
7.5  inches  in  width  and  about  8.5  in  length  with  40  data  points  per 
naqe  and  6  lines  oer  inch,  ""his  leaves  adeouate  margins,  when  trimmad 
to  8.5  by  11.0,  to  Place  the  pages  in  notebooks  and  reports.  ’"he 
page  is  repeated  as  reouired  to  accommodat®  the  number  of  Mach  number 
values  requested. 

The  output  is  identified  bv  the  two  lines  of  user-supplied 
headinq  entered  on  the  HFAD1  and  HFAD2  cards  and  printed  between  the 
dashed  lines  directly  below  the  standard  h®adinos.  The  next  informa¬ 
tion  is  the  location  of  the  lifting  surface  relative  to  the  center 
of  pressur®  of  the  anqle-of-at*ack  vane  and  *-he  geometry  of  the 
planiform  as  read  from  the  wiNFF  card,  ’"he  distances  and  dimensions 
are  printed  and  labeled  with  *h®  dimensional  abbreviation  input  in 
variable  tOUT.  As  with  program  BODY,  no  conversions  or  checks  ar» 
done  to  insure  consistency  of  units. 
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YACGY-KUGALLU  UP. ASH  CALCULATION  PAGE 

UPwASH  ESTIMATE  FQk  LIFTING  SURFACE  COMPONENTS  1/  1 


TRAPEZOIDAL  WING 

F 

-UlA/TAOl 

A I  R  C  R  A  F  T 

26  WING  SWEEP 

LIFTING 

SURFACE 

GEOMETRIC  Detail: 

LCNCIUOINAl  SIATIUN 

(JF  V A Nt  j  ., 

LONGITUDINAL  STATION 

UF  APEX  uF 

QUAR  TER- CHORD 

LINE 

997.2800  IN 

CISTANCt 

ERCM  APIX  JF 

GUARTtR-CHORl)  TO  VANES 

(XV). 

929.7300  IN 

SPAN  (PI 

669.6200  IN 

tau  (xv/ij/zn  . 

1.9839  B/2 

ASPECT  RATIL  . 

9.0700 

Sweep  cf 

QUA  R  T  EK  -LHUR 1)  LINE  ...  , 

23.3900  DEG 

LIFTING 

surface 

UP. ASH  ESTIMATES: 

EFFECT  i  Vt 

EPSILON  AR 

M  AC  H 

beta 

TAU/HETA 

.ing  sweep 

/CL 

EPSILON /CL 

EPSILON/CL 

<0t  CJ 

(  RADI 

(KAO) 

(DEG) 

c.coco 

l.OCCC 

1 .98  39 

2  i  •  J  9  3  J 

.029982 

.009026 

.287  76 

.  1CCC 

.09  SC 

i  .o<m 

7  3  M  <1 A  9 

. 029291 

.009978 

.289297 

./occ 

.9  798 

1  .fa  lb  1 

2  3 .  7bH9 

.079919 

. 009b  39 

.277033 

.  300 C 

.9929 

1  .0599 

29.  1  JC  7 

.02  379(3 

.009999 

.263289 

.  90c  c 

.olfci 

1./2/7 

^  1 1 1 

.0219(31 

.009257 

.293880 

•  6  00  C 

.  a66  C 

1 .8  2a  9 

2t . 989  7 

.01  yJ99 

.003816 

.218663 

.fcOoC 

.-3  0Cl 

1 .97  n 

2 1 .  :  9  1  i 

.0 lu‘>89 

.303271 

.187918 

.  /ccc 

.EMI 

7.71/1 

ll.Mlo 

.013263 

.002616 

.199880 

.HOCC 

.6  OC  L 

2  •  6  39 1 

19  .72/2 

.009 J61 

.001896 

.109786 

.scco 

.9  359 

3  . '  3/7 

AN.  709  8 

.009(389 

,0009t3 

.059191 

.so  c  c 

.Mil 

l  1  ./'/SO 

7  »  .  r.  /cl 

.000936 

.000106 

.006063 

FI OUSE  06 i  PROOMN  NINO  SAMPLE  OUTPUT 
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The  Mach  number  data,  effective  geometry,  and  upwash  estimates 
are  then  printed  with  one  line  ner  Mach  number  value  requested. 

The  first  and  second  values  are  the  Mach  number  and  corresoonding 
compressibility  factor  labeled  MACH  and  BFTA  respectively.  The 
effective  distance  of  the  angle-of-attack  vane  ahead  of  the  quarter 
chord  line  in  terms  of  r/°  and  the  effective  wing  sweep  for  the  par¬ 
ticular  Mach  number  value  are  then  nrinted.  The  upwash  parameter 
eAR/C  calculated  for  the  Mach  number  value  is  then  printed  as 
EPSILON  AR/CL  in  radians.  The  upwash  angle  estimate  ner  unit  lift 
coefficient  in  radians  and  then  deqrees  are  printed  as  EPSILON/OL 
as  the  last  two  values. 

HIHO  ERR  Oft  CHECK  I  NO-: 

Program  WING  performs  numerous  error  checks  as  data  is  read  and 
as  upwash  calculations  are  made.  The  main  purpose  of  the  error 
checking  is  to  avoid  abnormal  termination  of  the  program  or  to  give 
the  user  as  much  information  as  possible  on  the  location  and  nature 
of  the  error  if  termination  is  unavoidable.  Several  minor-error 
checks  are  made  which  change  internal  values  which  are  obviouslv 
wrong  and  would  cease  execution  of  the  program  immediatelv  if 
uncorrected,  '’’hese  changes  which  were  included  mainlv  to  protect 
the  user  from  gross  misapplication  or  keypunch  errors  are  not  accom¬ 
panied  by  error  messages.  Since  ♦hese  checks  might  limit  t*>®  exper¬ 
ienced  user,  give  unexpected  results,  or  caus°  confusion  due  to  apparent 
inconsistencies,  thev  are  discussed  briefly.  The  main  body  of  the 
checks  are  for  maior  errors  wi*h  input  data  which  do  not  have  easily 
determined  solutions  and  which  must  terminate  program  execution. 

These  maior-error  checks  write  error  messaqes  to  the  us»r  prior  to 
termination  to  identify  where  and  why  termination  occurred.  Both 
the  minor-error  and  major-error  checks  are  discussed  b“low. 

Minor-Error  rhecklnn.  Frrors  which  ano“ar  obvious  and  usually 
result  from  keypunch  errors  or  gross  misapplication  of  the  t°chniaue 
were  corrected  without  accompanvina  output  messaqes.  This  was  done 
to  allow  continued  execution  of  the  program  so  that  outout  could 
be  produced  to  aid  error  analysis  and  so  that  subsenuent  data  sets 
could  be  processed.  In  some  cases  this  results  in  invalid  or  unex¬ 
pected  results  which  should  be  investigated  and  fixed;  in  other  cases 
it  simply  allows  data  to  be  entered  in  a  convenient  manner  when  it 
would  normally  terminate  urogram  execution. 

Mach  number  values  entered  on  MACHV  cards  were  checked  for 
values  less  than  o.n  urd  greater  than  l.o.  Obviously  values  l*es 
than  h.O  have  no  meaning  and  usuall”  result  from  keypunch  errors. 

Values  less  *-han  o.O  were  cianged  to  0.0  irternallv  hut  printed  at 
their  oriqinai  value  to  allow  th»m  to  he  found  and  corrected.  Mach 
number  values  greater  than  i.o  ar*  eouallv  invalid  since  velocities 
cannot  be  induced  ahead  of  a  hodv  i"  supersonic  flow,  t’nwash  esti¬ 
mates  for  Mach  numbers  qreat<»r  t^nn  l.o  were  «=«t  to  zero  and  execu¬ 
tion  continued. 
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Major-F.rror  Checking,  Checking  for  maior  errors  which  neces¬ 
sitated  program  termination  was  accomplished  in  subroutine  DATARD. 

Any  major  error  discovered  resulted  in  an  error  message  followed 
by  program  termination. 

Normal  program  termination  resulted  from  readinq  an  end-of- 
. file  designator  while  attemotinq  to  read  a  HFAD1  card.  Unexpected 

reading  of  an  end-of-file  desiqnator  by  any  other  read  statement  results 
in  the  statement 

;  ’* 

:<*  "  "  IS  IN  ERROR  -  END-OF-FILE  ENCOUNTERED 

lr. 

where  the  quotes  have  the  five  alphanumeric  characters  of  the  card 
r'4  identifier  being  sought  when  the  end-of-file  designator  was  encountered. 

The  input  data  set  should  be  reassembled  with  the  proper  number  and 
types  of  cards.  A  related  message  results  i*  the  inout  data  set  is 
improperly  arranged  or  a  card  identifier  is  missoelled.  The  statement 

"  "  IS  IN  ERROR  -  CARD  READS  " 

will  be  printed  if  an  unexpected  card  identifier  is  encountered. 

The  first  set  of  quotes  encloses  the  five  alphanumeric  characters 
of  the  exoected  identifier  while  the  second  set  enclose  the  eiahtv 
alphanumeric  characters  on  the  card  actually  read.  Deck  must  be 
checked  for  proper  order  of  the  input  data  s«t. 

other  major  errors  result  from  errors  ir  input  data  values 
which  would  abnormally  terminate  execution.  If  one  or  more  errors 
are  found,  a  statement 

INPUT  DATA  FURORS: 

is  printed  followed  by  an  explanation  of  the  innut  data  errors.  The 
center  of  pressure  of  th°  an^le-of-a  :tack  vanes  must  be  well  ahead 
of  the  quarterchord  line  of  the  liftinq  surface.  If  th»  distance  is 
found  to  be  zero  a  messaqe 

DISTANCE  FROM  WINO  TO  VANF  MUST  BF  NON-ZFRO 

will  he  printed.  If  the  value  for  snan  of  the  lifting  surface  is 
found  to  be  neqative  or  7«ro,  a  message 

SPAN  MUST  BE  NON-ZrRO,  POSITIVE  VALUE 

will  he  printed.  If  the  value  for  aspect  ratio  is  found  to  be 
negative  or  zero,  a  messaoe 

ASPECT  RATIO  MUS*  BF  NON-ZFRO,  POSITIVE  VALUE 

will  be  printed.  The  wino  sweep  anqle  of  the  ouarter-chord  line 
must  be  greater  than  zero  (aft  sweep)  because  no  forward  sweep  data 
was  analyzed.  If  the  sweep  is  less  than  O.o  a  messaoe 

SWEEP  ANOLFS  LESS  THAN  0.0  ARE  Invalid 
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is  written.  If  the  sweep  is  greater  than  90.0  a  message 

SWEEP  ANGLES  GREATER  THAN  90.0  ARE  INVALID 

is  written.  If  no  Mach  number  values  are  soecified,  the  program  need 
not  be  executed;  therefore,  if  the  number  of  Mach  number  values 
soecified  by  variable  MACHS  is  zero  a  message 

NUMBER  OP  MACH  NUMBERS  MUST  BE  GREATER  THAN  0 
is  written. 

Another  message  is  written  if  the  center  of  pressure  of  the 
angle-of-attack  vanes  doesn't  lie  far  enouoh  ahead  of  the  quarter- 
chord  line.  The  production  software  uses  the  eouation  described  in 
Appendix  B  to  calculate  the  uovash  parameter.  As  explained  in 
the  derivation,  the  eouation  is  valid  only  for  x/S  values  greater 
than  0.4.  Therefore,  if  x/P  is  less  than  0.4  a  message 

'T'AU/BETA  VALUES  LFSS  THAN  0.4  ARE  INVALID 


is  written. 


APPENDIX  E 

PROGRAMMER'S  GUIDE  AND  PROGRAM  LISTINGS 
FOR  THE  YAGGY-ROGALLO  PRODUCTION  SOFTWARE 


INTRODUCTION 

Programmer  information  and  program  listings  for  the  tvc  programs 
developed  to  implement  the  Yaggy-Rogallo  upwash  estimation  technique  are 
included  in  this  appendix.  Both  programs  are  ANSI-standard  FORTRAN  V 
with  the  exception  of  a  single  file-declaration  list  on  the  first  card 
of  each  program.  Practically  no  user  or  programmer  guidance  is  supplied 
as  comments  within  the  code  because  of  the  extreme  simplicity  of  the 
code  and  the  production  nature  of  the  software.  A  user's  guide  to  input, 
output,  and  error  checking  for  both  programs  is  included  in  Appendix  D. 

Program  BODY  is  the  program  which  calculates  upwash  estimates  for 
bodies  of  revolution  using  equations  developed  in  Appendix  A  and  program 
WING  is  the  program  which  calculates  upwash  estimates  for  lifting  sur¬ 
faces  using  equations  developed  in  Appendix  B.  The  two  programs  are 
completely  independent,  but  there  is  great  similiarity  in  their  internal 
development,  input,  and  output.  Information  needed  to  understand  the 
working  or  limitations  of  these  programs  or  to  modify  them  is  included 
in  the  following  sections. 
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PROGRAM  BODY; 

Purpose.  Program  BODY  implements  the  Yaggy-Rogallo  technique  of 
making  upwash  angle  estimates  for  bodies  of  revolution.  The  body  is 
entered  as  a  series  of  longitudinal  station/cross  sectional  area  pairs 
which  describe  the  distribution  of  area  along  the  body  length.  At  each 
station  a  radius  is  calculated  which  is  the  radius  of  a  circle  having  an 
area  equal  to  the  actual  area  at  the  station.  The  upwash  angle  estimate 
for  a  segment  between  two  adjacent  stations  is  calculated  by  assuming  the 
area  of  intermediate  sections  varies  linearily  as  described  in  Appendix  A 
and  implemented  in  equation  (A40). 

Conditions  of  Validity. 

(1)  The  specified  location  of  the  center  of  pressure  of  the 
angle-of-attack  vane  must  not  lie  on  the  centerline  of  the  equivalent 
body  since  this  is  a  singularity  in  the  original  equations.  Mathe¬ 
matically  the  radial  distance  from  the  centerline  of  the  equivalent  body 
to  the  center  of  pressure  of  the  angle-of-attack  vane,  r  (variable  RVANE) , 
cannot  be  equal  to  zero. 

(2)  Program  dimensions  limit  the  body  to  999  segments  (1000 
longitudinal  station/cross  sectional  area  pairs). 

(3)  Program  dimensions  limit  the  number  of  Mach  number  values 
per  input  data  set  to  98. 

(4)  Mach  number  values  requested  must  be  greater  than  or  equal 
to  0.0  and  less  than  or  equal  to  1.0  since  other  values  have  no  meaning 
in  upwash  calculations. 

Storage  Required. 


PROGRAM 

COMMON 

TOTAL 


Subprograms  Used . 

ABS,  ATAN,  EPOAL,  RDDATA,  SQRT,  WRITE2 

COMMOM  Inputs.  Only  the  common  inputs  used  by  this  program  are 
described  be  1 ow , 


Octal  Words 

— ns — 

15701 

16051 


Decimal  Words 

— m - 

7105 

7209 
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COMMON/BBB/FSVANE , RVANE , ANVANE , NMACH , XMACH (98), FS (1000), AREA( 1000 ) 


INPUT 

DESCRIPTION 

UNITS 

r 

RVANE 

Radial  distance  from  the  centerline  of  the  equiv¬ 
alent  body  to  the  center  of  pressure  of  angle- 
of-attack  vane 

t  3 

a 

ANVANE 

Angular  location  of  the  center  of  pressure  of  angle- 
of-attack  vane  around  the  equivalent  body  center- 
line  with  0  being  vertically  above  the  centerline 

deg 

NMACH 

Number  of  Mach  number  values  where  upwash  angle 
estimates  are  desired 

N-D 

XMACH 

Mach  number  values  where  upwash  angle  estimates 
are  desired 

N-D 

COMMON/CCC/NPTS, NPGS, X( 1000 )  ,  R(  1000) 

INPUT 

DESCRIPTION 

UNITS 

NPTS 

Number  of  points  into  longitudinal  station/long- 
itudal  station  pairs  used  to  describe  the  equi¬ 
valent  body  (Number  of  AREAV  cards  read  by  RD- 
DATA) 

N-D 

X 

X 

Longitudinal  distances  from  center  of  pressure  of 
angle-of-attack  vane  to  longitudinal  station 

t 

R 

R 

Equivalent-circle  radii  at  each  longitudinal 
station 

i 

COMMOM  Outputs.  Only  the  common  outputs  used  by  this  program  are 
described  below. 


COMMON/DDD/XE( 1000 ) , THETA( 1000 ) , DELEOAf 1000 ) , SUM 
OUTPUT  _ _ DESCRIPTION _ 

x  XE  Effective  longitudinal  distances  from  center  of 

pressure  of  angle-of-attack  vane  to  longitudinal 
stations 

0  '’'HETA  Angle  between  body  centerline  and  line  connecting 

vane  and  center  of  equivalent  circles 

A( c/u) DELEOA  Increment  in  e/a  for  a  segment  of  tht  body 

I A ( c/a) SUM  The  sum  of  A(e/a)  values  to  give  the  e/a  for  the 
whole  body 

^Units  specified  as  units  of  length,  t,  are  not  in  any  particular  units  of 
length  but  all  parameters  having  units  of  t  must  be  in  the  same  unit  of 
length,  e.g.,  all  in  inches  or  all  in  centimeters. 


UNITS 

t 

rad 

N-D 

N-D 


r>r>  ->r»  n 


Program  Description.  Program  BODY  calculates  upwash  estimates  for 
bodies  of  revolution.  BODY  uses  subroutine  RDDATA  to  read  a  single  input 
data  set,  check  the  input  for  errors,  and  call  WRITE1  to  print  the  actual 
geometric  data  for  the  body.  The  effective  body  geometry  is  then  calcu¬ 
lated  from  the  actual  geometry  by  BODY.  Increments  in  upwash  angle  due 
to  each  segment  of  the  body  are  calculated  using  function  EPOAL  and 
summed  by  BODY  to  obtain  an  upwash  angle  estimate  for  the  whole  body. 
WRITE2  is  called  to  print  the  effective  geometry  and  upwash  estimates. 
RDDATA  is  then  called  to  read  another  input  data  set  or  terminate  exe¬ 
cution  . 

Program  Listing. 


RRHfiRAM  iHHITj J^TRUTtT APE 5 «INPUT#TAPF6 •OUTPUT! 


ERF  pFfl  L] 


NON-ANSI 


THIS 


PROGRAM  «nnv  IS  «  FORTRAN  V  RR06RAH 
fop  nnnTFS  OF  RFvni  yTION  USING  THE  YA 

OFPTVATint  nr  F"U»TI0MS  IJSF  ‘  - 

APE  CONTAIN*"  TV  "A  METHOD  OF 
NOUNTFO  VANFS-, AFFTC-TIH-il- 
OUESTTONS,  rnMMFNTS*  OR  REOO 
OTRECTEO  Tn  KEN  RAWLINGS*  65? 

TEST  rFNTFfl, cnuAPOS  REP,CA  93523  OR  TeL 
PFVIsjnN  pcrnant 
WRITTEN  NY  »AT  JUEN^NANN 


!TNG 
' D  IN 

DF  E _  _ 

*1  iJUNE  1981 
JESTS  FOR  TH 
J?0TG/E  NOT  I 


FgR^MAKING.  UPWASH  ..ESTIMATES 


RPGALLO  TECHNIOIIFj  THF 
PROGRAM  AND  A  USER'S  GUIDE 


ESTIMATING  UPWASH  ANGLE 


BODY 

BODY 

•BODY 

BODY 

BODY 

BODY 


AT  NOSFBnPM-BODY 
BODY 

USER'S  GUIDE  SHOULD  BE 
IP  239  It  AIR  FORCF  FLIGHT 
‘  f«05)  27Y-S7TO. 


TO 
•PHONE 


REVISFn  RV 


REN 


!NGS 


APR 

1900  1 

JUN 

1081 

10 


COMMON /MB  A*s  VANF>p  yAME? AMV  ANE  *  MMACH*  RMACHC  9B  )»  ES  ( 1000 )  »  A  RE  A(  1000 

F0A(1000)t  SUM 


20 


30 


*0 

50 


Common  /nnn  /v* i 
CALL  POOATA 
no  50  1 -1  •  nn a f h 

sun-0. 

AMACH-VNAnur  .1  i 

IFIANAOH  , l T ,  0.01  AMACM-0.0 
I  F (  AMA* H  .  C.T,  1,01  AMACH-  1.0 
BETA-SORTM.O-IAMACHRAMAChI  t 

on  30  t-i,n#ts 

Tf»xm  .1*.  0.0)  GO  TO  20 

XEm-XrilRBFTA 

THFTAfT 1-3,1 41SO>f/»,o 

IFfXEfT)  .NC.  0.01  THETA(IMATANCABS»RVANE)/XE(I)  1 
GO  TO  30 

XE ( I ) --999090080, 0999 

IFIRFTA  .NC,  n ,0)  XE II 1 -X( I) /BETA 

THETA(T1-3.141B9?a/3,c 

I*»XE(i)  .nc,  0.0)  THFTA(IMATANfAPS(RVAMF)/XE(I1»33.1B150?6 
CONTTNII* 

nn  40  T  •  3. N RTS 
DELFOAf T1-0.D 
JMPVANc  ,  *  0 , 


0.01  GO  TO 


I F  f THF  T  A  f  I  1  ,F0.  THETAfl-Ii?  GO  76  40 

OELFOAUI-cRnAL  fRII),  R»  1-1), THE  TAUIt  THETA!  1-1)  »R  VANE  tANVANF) 
SUM.SUM*"M  *"411 1 


f  DNT INIIF 

CALI  WRITE? 
CONTINUE 
GO  JP  10 
FND 


rxNArumi 


BODY 

BODY 

BODY 

BODY 

BODY 

BrpY 

BODY 

1BCDY 

BODY 

BODY 

BCDY 

BODY 

PPPV 

BODY 

BODY 

BODY 

BODY 

BODY 

BODY 

BODY 

BODY 

BODY 

BODY 

BCDY 

BODY 

BODY 

BODY 

BODY 

BODY 

BODY 
BODY 
BODY 
BODY 
BODY 
B  ODY 
BODY 
BODY 


FUNCTION  EPOAL; 


Purpose.  Function  EPOAL  calculates  the  upwash  estimate  for  a  single 
segment  o£  an  equivalent  body  of  revolution  using  the  linear-area  method 
recommended  in  Appendix  A  and  implemented  in  equation  (A40). 

Conditions  of  Validity. 

(1)  The  angle  at  the  leading  edge  of  the  segment,  ®i-l,  and 
the  angle  at  the  trailing  edge  of  the  segment,  °i,  cannot  be  equal. 

(2)  The  radial  distance  from  the  body  centerline  to  the  center 
of  pressure  of  the  angle-of-attack  vane,  r,  cannot  be  zero. 

Storage  Required. 


Octal  Words 


Decimal  Words 


70 


56 


Subprograms 

Used . 

COS,  SIN 

,  TAN 

Calling  Statement. 

FUNCTION 

EPOAL  ( RI,RIM1, THI, THIM1, RVANE, ANVANE) 

Calling  Argument  Inputs. 

__INPUT _ 

DESCRIPTION 

UNITS 

RI 

Radius  of  equivalent  circle  at  trailing  edge  of 
segment 

t 

1  RIM1 

P.adius  of  equivalent  circle  at  leading  edge  of 
segment 

l 

THI 

Angle  for  trailing  edge  of  segment 

rad 

1  THIM1 

Angle  for  leading  edge  of  segment 

rad 

RVANE 

Radial  distance  from  equivalent  body  centerline 
to  center  of  pressure  of  angle-of-attack  vane 

t 

ANVANE 

Radial  location  of  center  of  pressure  of  angle- 
of-attack  vani-  around  body  centerline 

deg 

Funct iona 1 

Output . 

Output 

DESCRIPTION 

UNITS 

:/a)  EPOAL 

Upwash  increment  for  a  single  segment  in  terms  of 
upwash  angle,  c,  per  angle  of  attack,  a 

N-D 

n 


Program  Listing 


rUNCTi™  e»n*l  ( «>T,  PIP1  ,THI  ,  THI  HI,  P  V*  NF  ,A*VANF) 

OPAL 

«l.f  (PTHl*oTni1-(Pi*pi))M(t./TAN(THIHl))-(l./TAN(THI)  )) 
*?«(9T*PT1-(M  /  t*N(  TNT)  > 

STNfTMT )-SIN(THTHl) ) )4(K  ?*( C05( TH IM U -r PS ( THT ) ) 1 
p  4CTHP  *  ( f  I  ptn  f  AWV»HF/'57.29*>P  »)**?)-  ((CHS  (  AN  VANE  /5  20*5"  >)*  +  ’)  ) 

1  /f  ?,*ovamp*PVANF ) 

FPnAl«P»Ctn'>*»!lt*M 

pun 


FPOAL 

FPOAt 

EPPAl 

F°OAl 

FPPAl 

PPOAl 

F»0*L 

FPflAl 

FPOAl 


\ 

3 

A 

5 

6 
7 

A 

9 


II 


SUBROUTINE  KQDATB s 


Purpose.  Subroutine  RDDATA  reads  one  input  data  set  per  call, 
checks  data  tor  input  errors,  and  passes  information  out  through  common. 


Conditions  of  Validity. 

(1)  Data  must  be  formated  in  80-column  card  images  as  described 
in  Appendix  D,  Users  Guide  for  the  Yaggy-Rogallo  Production  Software. 


Storage  Required. 

PROGRAM 

COMMON 

TOTAL 


Subprograms  Used . 

ABS,  SQRT,  WRITE1 
Calling  Statement. 
CALL  RDDATA 


Input  File  Inputs. 


INPUT 

DESCRIPTION 

UNITS 

HEADl 

First  line  of  user-supplied  heading  (75  alpha¬ 
numeric  characters) 

N-D 

HEAD2 

Second  line  of  user-supplied  heading  (75  alpha¬ 
numeric  characters) 

N-D 

UNIT 

Abbreviation  for  units  of  length  for  input  values 
(2  alphanumeric  characters) 

N-D 

FSV 

FSVANE 

Longitudial  (fuselage)  station  of  center  of  pres¬ 
sure  of  angle-of-attack  vane 

l 

r 

RVANE 

Radial  distance  from  the  body  centerline  to  the 
center  of  pressure  of  the  angle-of-attack  vane 

t 

A 

ANVANE 

Angular  location  of  the  center  of  pressure  of 
angle-of-attack  vane  around  the  body  centerline 
with  0  being  vertically  above  centerline 

deg 

Octal  Words  Decimal  Words 

— m —  - m — 

10030  4120 

10436  4382 


II 


1 


& 

n 


area  corresponding  to  stations  in  array  FS 


NMACH 

Number  of  Mach  number  values  where  upwash  angle 
estimates  are  desired 

N-D 

XMACH 

Mach  number  values  where  upwash  angle  estimates 
are  desired 

N-D 

• 

FS 

FS 

Longitudinal  (fuselage)  stations  used  to  describe 
equivalent  body 

l 

• 

AREA 

Longitudinal  (fuselage)  station  cross  sectional 
areas  corresponding  to  stations  in  array  FS 

2 

t 

COMMON  Outputs.  Only  those  common  values  output  by  this 

subroutine 

are  discussed  below. 

COMMON/AAA/HEADI , HEAD2 .UNIT 

OUTPUT 

DESCRIPTION 

UNITS 

HEAD1 

First  line  of  user-supplied  heading  (75  alpha¬ 
numeric  characters) 

N-D 

HEAD2 

Second  line  of  user-supplied  heading  (75  alpha¬ 
numeric  characters) 

N-D 

UNIT 

Abbreviation  for  units  of  length  for  input  values 
(2  alphanumeric  characters) 

N-D 

• 

COMMON/BBB/F  S VANE , RVANE , ANVANE , NMACH , XMACH ( 98 ) , FS ( 1 000 ) , AREA (1000) 

OUTPUT 

DESCRIPTION 

UNITS 

FSy 

FSVANE 

Longitudinal  (fuselage)  station  of  center  of  pres¬ 
sure  of  angle-of-attack  vane 

t 

r 

RVANE 

Radial  distance  from  the  body  centerline  to  the 
center  of  pressure  of  the  angle-of-attack  vane 

t 

il 

ANVANE 

Angular  location  of  the  center  of  presure  of  the 
angle-of-attack  vane  around  the  body  centerline 
with  0  being  vertically  above  centerline 

deg 

NMACH 

Number  of  Mach  number  values  where  upwash  a.'glo 
estimates  are  desired 

N-D 

XMACH 

Mach  number  values  where  upwash  angle  estimates 
are  desired 

N-D 

FS 

FS 

Longitudinal  (fuselage)  stations  used  to  describe 
equivalent  body 

t 

AREA 

Longitudinal  (fuselage)  station  cross  sectional 

2 

t 

COMMON/C CC/NPTS , NPGS, X ( 1000 ) , R( 1000 ) 


OUTPUT  DESCRIPTION  UNITS 


NPTS  Number  of  points  into  longitudinal  station  pairs  N-D 

used  to  describe  the  equivalent  body  (Number  of 
AREAV  cards  read  by  RDDATA) 

NPGS  Number  of  pages  required  for  NPTS  at  40  points  per  N-D 
page 


x  X  Longitudinal  distances  from  center  of  pressure  of  t 

angle-of-attack  vane  to  longitudinal  stations 

R  R  Equivalent-circle  radii  at  each  longitudinal  t 

station 

Messages, 

(1)  "  "  IS  IN  ERROR  -  CARD  READS  "  " 

Error  message  printed  when  RDDATA  encounters  an  unexpected 
card  identifier.  First  quotation  marks  contain  card  identifier  of  card 
expected  and  second  prints  80  alphanumeric  characters  on  card  actually 
read . 


(2)  LONGITUDINAL  STATION  IS  NOT  IN  ASCENDING  ORDER 

Error  message  printed  when  the  longitudinal  station  FS(I)  on 
an  AREAV  card  is  less  than  that  on  the  preceding  card. 

(3)  "  "  IS  IN  ERROR  -  END-OF-FILE  ENCOUNTERED 

Error  message  printed  when  an  end-of-file  is  encountered 
during  reading  of  any  card  other  than  a  HEAD1  card.  Quotation  marks 
contain  card  identifier  of  card  expected  when  end-of-file  was  encountered. 

(4)  NUMBER  OF  LONGITUDINAL  STATIONS  MUST  EXCEED  1 

Error  message  printed  when  only  one  AREAV  card  is  read 
and  NPTS  equals  one. 


IS 


Program  Listing. 


SUBROUTINE  RDDATA 

£oMMON/fi$$/£fvWl#lv^I*AftvA)*E#NMACH#XMACH(9B)#FS(l©|De)#AREA( 

^HARACTiSCHEAoii?5fHtAoil8?ii^EAi3*75»HEAD*5»HEAOER*5»UNIT*2 


READ  *HEAD1*  CARO  - 

HEADER-'HEADl* 

RE AD(9#100#END«90)  H6AD#HFA01 
IF  ( HE  AO  .NE.I  HEADER)  GO  TO  5*. 
Rt  AD  ' HE AD2 '  CARO - 


1/  »¥  »  fifcmw/iii. 

HEADER)  GO_TO_5t_  _ 


HE  A0ER»fHEA02' 

RE  AOO*1C  #END"80)  HE AD#  HEAD2 

itw*uw  - 

HEADER"*  VANES' 

RF  AD (5* 110# END" 80)  HE  AD# UNIT, FS WANE #R VANE# AN VANE 
IF (HEAD  .NE.  HEADER)  GO  TO  9C 

READ  ‘HACHVT  CARO  - 

HEADER"*  MACHV* 

03  1»>  J"l#  92#  7 

READ(5*120*EN0-8O)  HE AD# MACH# ( XNACH ( I ) *  I "J* ( J  +  6 ) ) 
IF (HEAD  «NE«  HEADER)  GO  TO  it 
IF  t  J  • EQ«  1)  NHACH-NACH 
IF (NMACH  «LE •  (J+6))  GO  TO  2 C 
1/  CONTINUE 

READ  *AREAV*  CARDS  - 

HE  ADER"' AREA  V' 

NPTS"0 

RE k-  ?5 0"B0)  HEAO#F  S (I >  #  AREA  < I ) 

IF  (HEAD  •E0«>'  END  •  )  GO  TO  AC 


IF (HEAD  ,NE,  HE 
NPT$"NPTS*:. 

X(l)  •  FSVANt  - 


HEADER)  GO  TO  )(, 


K ( 1 1  •  F S VANE  -  FS1I) 

R< I )  -  SQRT( ABS ( ARE  A ( I ) )  /  3. 
If  II  *E0«  li  GO  TO  30 
IF(FS(I).LT.iFS(I-ll  )  GO  TO  60 
3  *  CONTINUE 
4j  HPGS"( NPTS /AO ) ♦! 

IF(NPTS  .EG.  1)  GO  TO  7l 

CALL  WRITE! 

RETURN 

j  •  BACKSPACE  5 

Re  AD  (9# 100)  HE AD*HF  AD3 


3.1A15926) 


RP  IT F ( 6# 1 AO)  HEADl 

STOP 

63  WR I TF (fa# 190 ) 

STOP 

n  V* I TE ( 6# 160) 

STOP 

4)  WR I TE ( 6# 170)  HEADER 
V3  STOP 

FORHAT  STATEMENTS 
FJRMAT(A5#A75) 


R  »  HE  A  D,  HE  A  D3 


FORMAT 

FORMAT 


5  IS  IN  ERROR  -  CARD  READS  "'#A5#A7 

IV,  HgBiIJiSi;:)i8»fIiu8INfWI)IiyS*is5mii»siS6i¥siSfE?1 

17  »  F3«MAT11H1#*,,*#A5#'"  IS  IN  ERROR  -  ENC-OF-FILE  (NCOUNT 

EHD 


IS  IN  ERROR  -  CARD  READS  " *#A5# A75 * • •• ) 

NtWi)ISfS.Psm.5lls‘I6Sf5{5!E?5I,£S:i 

IS  IN  ERROR  -  ENC-OF-FILE  (NC OUNTE RED' ) 


RDDATA  1 
RDDATA  2 
iaoa)RODATA  3 
RDDATA  A 
RODATA  9 

- RDDATA  6 

RDDATA  7 
RDDATA  8 
RDDATA  9 
- RD0ATA10 

5881' it) 
- SBttliH 

RD0ATA15 

RDDATA16 

X.DDATA17 

- RDDATA18 

ROD  AT  A1 

rddat  a|o 

RD0ATA21 
RDDATA22 
RDDATA23 
RDDATA2A 
KOOAT A25 
- RDDATA26 

SSBiliiJ 

RDDATA29 
KDDATA30 
KDDATA3I 
RDQATA32 
RDDAT  A33 
RDDATA3A 
RDDATA39 
RD0ATA36 
RDDATA37 
RDDAT  A36 
RDDATA39 
RDOATAAO 
kDDATAAI 
RDDATAAZ 
KDDATAA3 

rdoataaA 

RDDATAA9 

RCDATAAb 

RDDATAA7 

RDDATAAB 

RDDATAA9 

RDDATA93 

RDDATA91 

RCDATA92 

- kODATA93 

ROD AT  A3A 
RDDATA99 
KD0ATA96 

E88itt?« 

!88ifi?» 

R00ATA61 

K0DATA62 


f 

m 


$ 

2 

8 


8UW0UTIME  WHITE  I,  i 

Purpose.  Subroutine  WRITE1  prints  the  actual  geometric  data  for 
the  body  with  40  points  per  page  and  up  to  25  pages  as  required.  One 
set  of  pages  is  printed  for  each  input  data  set  of  the  input  file. 


Storage  Required. 


PROGRAM 

COMMON 

TOTAL 


Octal  Words 

— m - 

10030 

10441 


Decimal  Words 

- 51? - 

4120 

4385 


Calling  Statement. 

CALL  WRITE1 

COMMOM  Inputs.  Only  the  COMMON  inputs  used  by  this  subroutine  are 
described  below. 


COMMON/AAA/HEADl , HEAD2 , UNIT 


INPUT 


DESCRIPTION 


HEADl 

First  line  of  user-supplied  heading 
numeric  characters) 

(75  alpha- 

HEAD2 

Second  line  of  user-supplied  heading 
numeric  characters) 

(75  alpha- 

UNIT 

Abbreviation  for  units  of  length  for 
(2  alphanumeric  characters) 

input  values 

UNITS 

N-D 

N-D 

N-D 


COMMON/BBB/FSVANE , RVANE , ANVANE, NMACH , XMACH( 98) , FS( 1000) , AREA( 1000) 


INPUT 


DESCRIPTION 


UNITS 


FSV  FSVANE  Longitudinal  (fuselage)  station  of  center  of  pre¬ 
ssure  of  angle-of-attack  vane 

r  RVANE  Radial  distance  from  the  body  centerline  to  the 

center  of  pressure  of  the  angle-of-attack  vane 


ANVANE  Angular  location  of  the  center  of  pressure  of 

angle-of-attack  vane  around  the  body  centerline 
with  0  being  vertically  above  centerline 

NMACH  Number  of  Mach  number  values  where  upwash  angle 
estimates  are  desired 


deg 


N-D 


•7 


1 


X 


XHACH 

Mach  number  values  where  upwash  angle  estimates 
are  desired 

N-D 

FS 

FS 

Longitudinal  (fuselage)  stations  used  to  describe 

Z 

equivalent  body 

AREA 

Longitudinal  (fuselage)  station  cross  sectional 
areas  corresponding  to  stations  in  array  FS 

£» 

% 

COMMON/CCC/NPTS , NPGS ,  X  ( 1 0  00  ) ,  R ( 1 000  ) 

INPUT 

DESCRIPTION 

UNITS 

NPTS 

Number  of  points  into  longitudinal  station 
longitudinal  station  pairs  used  to  describe  the 
equivalent  body  (Number  of  AREAV  cards  read  by 
RDDATA) 

N-D 

NPGS 

Number  of  pages  required  for  NPTS  at  40  points 
per  page 

N-D 

X 

X 

Longitudinal  distances  from  center  of  pressure  of 
angle-of-attack  vane  to  longitudinal  stations 

Z 

R 

R 

Equivalent-circle  radii  at  each  longitudinal 
station 

Z 

Printer  Outputs. 

INPUT 

DESCRIPTION 

UNITS 

NPGS 

Number  of  pages  required  for  NPTS  at  40  points 
per  page 

N-D 

HEADl 

First  line  of  user-supplied  heading  (75  alpha¬ 
numeric  characters) 

N-D 

HEAD2 

Second  line  of  user-supplied  heading  (75  alpha¬ 
numeric  characters) 

N-D 

UNIT 

Abbreviation  for  units  of  length  for  input  values 
(2  alphanumeric  characters) 

N-D 

r 

RVANE 

Radial  distance  from  the  body  centerline  to  the 
center  of  pressure  of  the  angle-of-attack  vane 

l 

ANVANE 

Angular  location  of  the  center  of  pressure  of 
angle-of-attack  vane  around  the  body  centerline 
with  0  being  vertically  above  centerline 

deg 

N 


**« 


FSVANE  Longitudinal  (fuselage)  station  of  center  of 
pressure  of  angle-o£-attack  vane 

FS  Longitudinal  (fuselage)  stations  used  to  describe 

equivalent  body 

AREA  Longitudinal  (fuselage)  station  cross  sectional 
areas  correrponding  to  station  in  array  FS 

X  Longitudinal  distances  from  center  of  pressure 

of  angle-of-attack  vane  to  longitudinal  stations 

r  Equivalent-circle  radii  at  each  longitudinal 

station 


*1 

e 

j*r 


Program  Listing. 


SUBROUTINE  WRITEI 

(!0Mm8n/BB6 /f10X&£>  ^VA$£»  AlSV  ANE,  NHAC  H>  XHACH(  98) »  FS  ( 1006) t ARE  A ( 1000 
CONMON/CCC /NPTS, NPGS»X  CloQQ ) ,R (lO&O ) 

01  MENTION  l INE6 ( 2  I » LINE7(Z1 »LI NE8(2 I 

CHAF.ACTfc R  HEADi*75,HEAD2*75#UNIT*2,TITLEl*32.TITlE2*41 
CMA'-.ACTER  LINe!*30»lIne2*5©»LINE3*30#l!nE4*50, LINE 3*28,  LINE 6*34 
CM  AW  A^TE  R_tINE7*34»  L!NE8*34*  LINE9*5 2»DASH4*4,DASH1 3*13*  DASH 13* 15 

DATA  Tir 


WRITEI 

>m\ 

WRA6E1 

write: 


TER  LINE?*34,LINE8*34»LINE9*52,DASH4*4.DASH13*13, 
ITLE1/ ' YA6GY-R06ALLC  UPWASH  CALCULATION'/ 

ITLE*/ ' GEOMETRIC  DATA  FOR  CYLINDRICAL  COMPONENTS' 

DATA  DASH4,pASH13,DASH15/' - ',' - - 

DATA  LINEI/' ANGLE-OF-ATTACK  VANE  LOCATION!*/ 

DATA  LINE2/ 'RADIAL  DISTANCE  FROM  BODY  CENTERLINE  . 

DATA  LINE3/' ANGULAR  LOCATION  . . 

-  LINE4/' LONGITUDINAL  STATION  CF  VANES  ••••• 

L INE5 /' C ONPONENT  AREA  OISTRIBUT ION» • / 

CROSS  SECTIONA* 
LONGITUDINAL  ECU* VALENT  • 

LONGITUDINAL 
DISTANCE  FROM 


wm  * 


STATION 
VANE  TO  STATICN 
)  (SQ 


AREA  OF  BODY 
RADIUS  AT 
AT  STATION 
STATION 
I 


•/ 

(  ) 


FUR  GEOMETRIC  DATA 


DATA 
DATA 

DATA  LINE6/' 

1  'L 

DATA  LINE7/' PT 
1  • 

DATA  LINEB/'NO 
.  • 

DATA  LINED/' (  ... 

PRINT  GEOMETRIC  DATA  FUR  CYLINDERICAL  COMPONENTS  - 
LPT-0 

DO  20  Ial,23 

^EfflEfi 

*>  TITLE2,I,NPGS 

>)  DASH13,DASH15,DASH'.3,DASH13,DASH15 
II  HEAD! 

head! 

DASH15,DASH1«,DA$Ha5»DASH15,DASH15 
LINEI 

LINE2»RVANE»UNIT 
L INE3, AN VANE 
LINE4»FSVANE»UNIT 
LIME5 
LINE6 
LINE7 
LINES 
LINE9 

UNIT, UNIT, UNIT#  UNIT 
DASH4,DASH13,DASH13,DASH13, DASH13 

PAGE  OF  GEOMETRIC  DATA  - 

,40 


•••••••'/ 

. '/ 

. . 


WRl 

WRITE! 
WRITE! 

write!: 

'/WRITE 
WRITE 
WRITE 
WRITE.. 
WRITEI! 

write!! 


WR 

WR 

WR 

(  )'/WR 

W- 


PRINT  PAG. 
WR I TE (6, It 
WRITE(6,‘‘ 
WR I TE (6, . 
WRITE(6,13l 
WR  I  T  E  ( 6, 13i 
WRITE  (6,121 
WRITE  6,  141 
WR ITE (6, Isl 
WRITE(6,16{ 

mim 

URITE(6,18( 
WRITE (6,If 
W*ITE(6,: 
WRIT6(6, 


WR  ITE(6, 

W? ITE  < 6, 

PRINT  ON 
OJ  10  J- 
LpT*LPT* 

IF (LPT  • 

WR ITE<6,2 
1"  CONTINUE 
20  CONTINUE 
RE  TURN 

FORMAT  STATEMENTS - - — - — 

F  JR  NAT (jH1,T24»A32, T73, 'PAGE'  I 

U9^A41,T72,I2,'/',I2) 


>T  fNPTS  I  RETURN 
1201  LPT,FS(LPTI, 


AREA(LPT1,XCLPTI,R(LPT) 


WR'GE14j 
WRlTul4< 
WRITEf 
WRITE] 

- WRI0E143 

WRITE146 
WRITEJ.R 

m 

write: 
write: 


I'JD 


I1J  F JRMAT (IH 
1/5  FORMAT  *Im 


FORMAT 
13  J  F3RMAK 
14J 


formaUIh 


ll*» 

lli 


.  H 


FORMA 

193  FORMATC. 
JjO  F3RNAT(!h» 
F OR  MAT (IH 
FORMAT (IH 
HD 


■»  »»  w 

ft* 


5eg>  i 


SAlIl 

f  aIoI 

FORMAT (iH  ,T10,A30,?X,F10.4,1X,A2) 

FOR  MAT (IH  ,TIo,A30,3X,f!o.4,1X, 'OE 

IP’ 

T  ),A2,16X,A2,13X,A2,13X, A2) 

T09, A4,2X, A1 3# 2X,A13,2X, A15,2X,A13) 

7X, !4,2X,F13«4,2V,F13.4,2X,f!?.4,2X,F13*4] 


: 


8UM0UT1MC  HHITC2 : 


Purpose.  Subroutine  WRITE2  prints  the  effective  geometric  data  for 
the  body  with  40  points  per  page  and  up  to  25  pages  as  required.  One 
set  of  pages  is  printed  for  each  Mach  number  where  upwash  estimates  are 
requested . 

Storage  Required. 

Octal  Words  Decimal  Words 
PROGRAM  Tfl  TSl 

COMMON  11633  5019 

TOTAL  12202  5250 

Calling  Statement. 

CALL  WRITE2(AMCT) 

COMMON  Inputs.  Only  the  COMMON  inputs  used  by  this  subbroutine 

are  described below. 


COMMON/ AAA/HEADI , HEAD2 , UNIT 

INPUT 

DESCRIPTION 

UNITS 

HEAD1 

First  line  of  user-supplied  heading 
numeric  characters) 

(75  alpha- 

N-D 

HEAD2 

Second  line  of  user-supplied  heading 
numeric  characters 

(75  alpha- 

N-D 

UNIT 

Abbreviation  for  units  of  length  for 
values  (2  alphanumeric  characters) 

input 

N-D 

COMMON/CCC/NPTS, NPGS, X ( 1000 ) , R( 1000 ) 

INPUT 

DESCRIPTION 

UNITS 

NPTS 

Number  of  longitudinal  station  longitudinal  cross 
sectional  area  pairs  used  to  descirbe  the  equival¬ 
ent  boly  (number  of  AREAV  jards  read  by  RDDATA) 

N-D 

NPGS  Number  of  pages  required  for  NPTS  at  40  points 
per  page 

X  Longitudinal  distances  from  center  of  pressure  of 

angle-of-attacfc  vane  to  longitudinal  stations 


N-D 


COMMON/DDD/XE( 1000 )  ,THETA( 1000 ) , DELEOA( 1000 ) , SUM 


1 


Kl 


OUTPUT 

DESCRIPTION 

UNITS 

X 

XE 

Effective  longitudinal  distances  from  center  of 
pressure  of  angle-of-attacfc  vane  to  longitudinal 
stations 

l 

e 

THETA 

Angle  between  body  centerline  and  line  connecting 
vane  and  center  of  equivalent-circle 

rad 

A( e/a) 

DELEOA 

Increment  in  e/a  for  a  segment  of  the  body 

N-D 

I A ( e/a) 

SUM 

The  sum  of  A (e/a)  values  to  give  the  e/a  for  the 
whole  body 

N-D 

m 


Program  Listinc 


SlIPROUTIwe  W»TTF?UNC7) 

comnon/aaa  /hfaoi,heao2*unit 

COMMON /Cf  C /N«*TS, N-Pp$.  X(  1000  l»R(  1000) 

r  nuHnv  f  ftn-) /vc  n  r>oo). TMf  T*  (looo)  iOfiFOAj  1000  )»sw 

niMENSTOM  i  TNEl  (?T»llNE?(?)»lINf3(2J»tINM(2> 
CHARACTER  TTTm*22,TITLE2*42.HEA0I*75»HEAD2*75i 
CHAP  if  tcb  l^Mci*i«,LlNE2*35,LINE3*35»LINM*V 
f  HAP  AC  TCP  r)ACMA*4,OASW6*6,0  4SH10*10*OASHlf*15 
n  ATA  TTTLel/»VAfRV-pn6ALL0  UPWASH  CALCULATION*/ 
maTA  TIT(e?f,t'°W ASH  ESTIMATE  FOP  CYLINDRICAL  COI 


IINIT*2 


TTT(  e?/t(iowAMH 


NATE  EO*  CYLINDRICAL  COMPONENTS • / 


DATA  PASHA, n»SH/., PASHIO/*— '  »♦ - - '  ,l  — - - 


fASMU/t - 


PITA  I TNC1 /* 


LONGITUDINAL 


OINAL  EFFECTIVE 
EQUIVALENT  I  NCR*  ME  NT 


OATA  ITNC»/*PT  "TSTANCE  FROM  DISTANCE  EPPM*, 
l  *  RADIUS  AT  IN  EPSILON  •/ 

PATA  ITVC1/*NP  VANE  TO  STATION  VANE  TO  STATIS 
L  *n>j  tufta  STATION  nvER  ALP«A  »/ 

DATA  I TNCfc /i  f  )  (  )  *, 

L  *  (PAD!  f  1  »/ 

PRINT  IIPVASM  c?TT(.ATFs  COP  CYLINOFRICAL  COMPONENTS 
IPT-O 

no  PO  T«1,P« 

PPTNT  PAfc  *c*0Tkif.5  CfJP  IIPWASH  ESTIMATES  — - - 

VPTTE(A,!f>Ot  TTTLE1 
VPITE(A,1V.1»  TTTirJ,  l,N»r,S 
VPTTF(A,1»0*  AMCT 

WPTTF( 4,1*0*  n»CHl5,PASHl5» 0ASH15, PASH1 5, PASHi* 
VPITE (A,l*n)  ue#r>T 
vRTTEM.Un)  mc A n P 

WPTTF( 4,130)  PASH15,PASH15»PASM1?,  PASm15,PAShi*; 
WPITF(A,1?0*  LTMCi 
WPITF(A,1S0*  iTNr* 

WPTTF(4,1*0*  ITwss 
VPITCfA.HO*  11**4 


VRITEfA.lVrt) 
VPTTF (4,1*0* 
WRTTF(4,130* 
VPITE(A,1*'M 

WRTTEM.IAO) 
WPTTE ( A ,t SO* 
WPITF(A,1?0* 
WPITF(A,1S0* 
WPITFM.HO* 
VRITCfA.lSO* 


WPTTE(A,l«.OI  "HTT,UNIT»UNTT 

WPITF (4, 170)  PASua, PASH1S,0ASN15,0ASHE,DA 

PPJNT  "N*  PASC  nc  UPWASH  ESTIMATES  - 


SHI 0,0  ASH10 


no  10  *«I,4D 
| PT.IPT*! 

TC(I  PT  ,<“T,  M»T*1  op  TP  TO 

WRITE l PT, V( LPT), XE( LPT)* THETA (L PT I, R(LPT) , DELE OA( L »T) 
rnwTiwMc 

CnHTTNUC 

WPITF  *, 10*>t  PA^uto 
WPTTE ( A, ?0Q)  SUM 
RETURN 

FORMAT  STATCMCWT*' - - - - - - - - - . - - — — — — — 

FOPMSTMH|,T*A,  ASP,  T73*  pp  AG  E* ) 
eORMATflM  , T1  0, A A 2. T7?.  I P  ,  •  /*,!?) 

E0PMAT*1M  , T*0, ' MAC  H  NUMACP  ■  *,E6,A) 

CPOM  AT  (  1H  ,*A1M 
EOPMATdH  ,at«) 

FORMAT ( in  ,T*.|*SAV) 

FnRNAT(lH»,TPO,AP,l5X»A2»20X»A?( 

FOPMATMH  ,T7,«a,px*A1«,pX,  A1S.2X,  AF»  ?X#  A  10.2X,  At 0> 

F0PNAT(1H  ,TTLTA,?X,Fl5.A,2X,Ell,A,2X,E6. A,2X,F10,A,PX,C10,4) 

FPPMATMH  .TaSAIO) 

FORMAT  ( lu  ,TPT,  ’TOTAL  UPWASH  FSTIMATE  (EPSI  LON/ AL  P**A)  •  *,F10.4) 

FNO 


FOPMATMH  ,T7,»a,?X*A1«,PX,  A1S.2X,  A6, 2X, A 10.2X, At 0> 

F 0®M AT ( 1 H  ,TTLT4,?X,Fl5.4,2X,Ell,4,2X,E6.4,2X,E10.4,PX,E10,4) 
COPMATMU  .TA7.A10) 


SSItEi  i 

VRTTF2  I 
VRITFI  4 
VRITeI  9 
WRITE?  6 
WRITE?  7 
WRTTF2  8 
WPITF 2  0 
WRITE 210 
WRITE 211 
WRITE 
WRITF213 
WPTTF214 
WRITE219 
WPITE216 
WR  TE217 
WRITF210 
WRITE?!* 
VRITF220 
— VRJTF2?! 
WRITE??? 
WRITE  ?23 
—WRITE??* 
WRITE225 
WRTTF  226 
WRITEP27 
WRTTF  ??§ 
WRITE??* 
WPITE230 
WPITF  231 
WPITE23| 
WRITE 233 
WRJTF234 
WRITF239 
WRTTEPSft 
WPITE237 
—WRITE  230 
WRITE  239 
WRITF2A0 
WPITF241 
WR  TF24? 
WR  TF 243 
WP  TF244 
WRTTF24? 
WR1TE240 
WRITE ?47 
— WRITE240 
WR!TF|a9 

WRITF250 

WPITEfJl 
WRITFil? 
WPTTE253 
WRTTF234 
WRITF253 
WRTTE256 
WRITF297 
WRITE??* 
WR!  TF259 
WRTTF 260 


riTOMftH  MI  NO 

PWOOKBM  NINO: 


Purpose.  Program  WING  implements  the  Yaggy-Rogallo  technique  of 
making  upwash  angle  estimates  for  lifting  surfaces.  The  lifting  surface 
is  estimated  by  a  thin  airfoil  with  trapezoidal  planform  as  described  in 
Appendix  B  and  the  upwash  angle  is  calculated  using  equation  (B8). 

Condition  of  Validity. 

(1)  The  lateral  location  of  the  angle-of-attack  vane  must  be 
very  near  the  midspan  of  the  lifting  surface. 

(2)  The  longitudinal  location  of  the  angle-of-attack  vane  must 
be  ahead  of  the  quarter-chord  line  of  the  lifting  surface.  Mathematical 
calculations  require  only  that  the  value  be  non-zero  to  continue  program 
execution  but  the  basic  assumptions  of  the  equations  require  the  effec¬ 
tive  distance  to  be  at  least  0.4  of  the  span  of  the  lifting  surface. 

(3)  The  specified  span  and  aspect  ratio  of  the  lifting  surface 
must  be  a  non-zero,  positive  values. 

(4)  The  specified  quarter-chord  sweep  angle  must  be  greater 
than  or  equal  0.0  ana  less  than  or  equal  90.0. 

(5)  Program  dimensions  limit  the  number  of  Mach  number  values 
per  input  data  set  to  98. 

(6)  Mach  number  values  requested  must  be  greater  than  or  equal 
to  0.0  and  less  than  or  equal  to  1.0  since  other  values  have  no  meaning 
in  upwash  calculations. 

Storage  Required ■ 

Octal  Words 
PROGRAM  H75 

COMMON  1266 

TOTAL  1426 

Subprograms  Used. 

ALOG10 ,  ATAN,  DATARD,  SORT,  TAN,  WRITE 3 

COMMON  Inputs.  Only  the  COilMON  inputs  used  by  this  program  are 

described  below. 

COMMON/EEE/FSVANE, FSQCL, XV, SPAN,TAU, AR,WSQCL, NMACH , XMACH( 98) 


Decimal  Words 

56 

694 

790 


m 


INPUT 

DESCRIPTION 

UNITS 

FSV 

t 

FSVANE 

Longitudinal  (fuselage)  station  of  center  of 
pressure  of  angle-of-attack  vane 

l 

psc/4 

FSQCL 

Longitudinal  (fuselage)  station  of  quarter- 
chord  line  at  midspan 

i 

xv 

XV 

Longitudinal  distance  from  midspan  of  quarter- 
chord  line  to  longitudinal  station  of  center  of 
pressure  of  angle-of-attack  vane  Xv  *  FSc/4~FSv 

l 

b 

SPAN 

Span  of  the  lifting  surface 

t 

1 

T 

1 

TAU 

Dimensionless  distance  from  quarter  chord  to 
center  of  pressure  of  angle-of-attack  vane 

N-D 

AR 

AR 

Aspect  ratio  of  the  lifting  surface,  AR  «  b2/S 

N-D 

Ac/4 

( 

WSOCL 

Sweep  of  the  quarter-chord  line  of  the  lifting 
surface 

deg 

NMACH 

Number  of  Mach  number  values  where  upwash  angle 
estimates  are  desired 

N-D 

XMACH 

Mach  number  values  where  upwash  angle  estimates 
are  desired 

N-D 

COMMON/FFF/BETA( 98 ) , TOB { 98 ) , WSB{ 98 ) , EARCL( 98 ) , ECLR{ 98 ) , ECLD( 98 ) 

0 

i 

BETA 

Compressibility  parameter,  8  «/l-M2 

l  -D 

r/B 

TOB 

Distance  ahead  of  quarter-chord  line  corrected 
for  compressibility 

N-D 

a3 

j 

WSB 

Sweep  of  quarter-chord  corrected  for  compressi¬ 
bility  tan  Ag  ■  tan  Ac  /4/B 

deg 

eAR/C L 

EARCL 

Lifting  surface  upwash  parameter 

t*-D 

!  e/cL 

ECLR 

Upwash  angle  per  unit  lift  coefficient  in  radians 

1  ad 

e/rL 

ECLD 

Upwash  angle  per  unit  lift  coefficient  in  degrees 

deg 

Program  Description.  Program  WING  calculates  upwash  estimates  for 
lifting  surfaces.  WING  uses  subroutine  DATARD  to  read  a  single  input 
data  set  and  check  th»'  input  for  errors.  It  then  calculates  upwash 
estimates  for  all  Mach  numbers  requested  and  calls  WRITE3  to  print  the 
effective  geometry  and  upwash  estimates.  DATARD  is  then  called  to  re*d 
another  input  data  set  or  terminate  execution. 


U  1 


Program  Listing. 


PRPGPAN  WTNG(TNPHT, OUTPUT# TAPF 5 ‘INPUT  a TAPF6-0UT PUT) 

PILE  PeCL  4p  ATTON  LIST  NON-ANSI 

PROGRAM  WTHG  h  7”fORTR AN  V  PROGRAM  FOR  MAKING  UPWASH~fSTIMATFS 
pnp  LIFTING  Giipr  AOf  s  USING  THF  YAGGY-POGALLO  TECHNIQUE.  THE 
PEPIVATing  nr  fomatIONS  USEP  IN  THIS  PPOGPAH  AND  A  USER'S  GUTDF 
A»F  Ff'NTATNcn  yw  "A  METHOD  OF  ESTIMATING  UPWASH  ANGLE  AT  WOSF RO0H- 
NnUNTFn  VANPP",  ACFTC-TIM-fll-1#  JUNE  19  Ml. 

OIJFS  TI TNP,  rnuMCNTS#  OR  REOUESTS  FOP  THE  USFP*S  GUIDE  SHOULD  RF 
PIPFCTFn  tn  KPN  PAULINGS*  6520TG/ENDT  (STOP  ?39)»ATP  FOOCF  CL TGWT 
TPPT  'NTEP.FPUAPns  AFB»C A  <53  5  2  3  0®  TfLfPHONC  (POST  277-3779. 
PEVISI  iN  pc-'npnt 

WRITTEN  P Y  KFn  RAWLINGS  J UN  1981 

rnMMON/ccP /PFVAMc7FSOri 7*v7sPAN,TAU»AR, WS QC l » NM AC H , XMAC H( 9R ) 
COMMON  zfcf  /ncTA  f  «R)  » TOP (98)  »WSR  (98  )  *F  AR  CL  (9  8)  *  FCL  R  fRfl  )  .  FT L  D C 9R ) 
PATA  A.“»C /-l.APPO7301»-0.008AA7R6P»-I.0993t«6RAA 
10  CALL  n *TA»n 

TF(NMACH  .FO.  0)  GO  TO  10 

TAt'»X\'/ (RB  AN /?,  1 

IFopnp.o 

no  50  T-l.NMATM 

T  F ( X  MAGHf  T 1  ,GF.  0.0)  GO  TO  20 
PETAfT)  »V.O 
g n  Tn  4 o 

?0  I F ( XHA r  Hf  < )  . L 7 .  1,01  GO  TO  30 
nETA ( I )  •0.0 
TORfI)  »Qonoo.Po«g 
wGPfi)  «on,n 
F  ARCLfT  )-7.0 
FCI R  f I  )  *5.0 

rCLOm  -n.n 
GO  TO  "0 

7 0  RFTA(n-7DPTfl.-fXMAC*<(I)*XMACH(I)  )) 

AO  TOPfl)  .TA'I  /BCTA  f  T  ' 

Tc(TnofTi  .it.  o,4)  i FPPno.i 

W  FR  (  1  )  ■  ( A  T  AN  f  TA  A'  (  W$  OCL  /  57  ,?  Q5  A  )  /RF  T  A  (  I)  ))  *5  7.  20  58 
F  ARC  L  (  T  )■)  n,0**f  (A*ALOG10(T0Rmi)*(«RWSR(I  ))*C  ) 

FCI  R  (T  )  »f A 0 I  f  T  )  /AP 

rciAfTliF'i  *(’**<7, ?Q5R 
■0  rpNTTNMF 

PALI  V9  TTr  A 

Ic(IPPRnp  ,rz>,  i)  WPT  TF  f  5*  1 00  ) 
on  TO  m 

100  FOPMATf  l«n,*TA"/oPTA  ''ALOES  LFSJ  THAN  O.A  A  PF  INVAIIDM 
F  NO 


WING 

WING 

•WING 

1 

3 

WING 

A 

WING 

5 

WTNG 

6 

■WING 

7 

WING 

R 

WTNG 

9 

WING 

10 

WING 

11 

WING 

1? 

WING 

13 

•WING 

1A 

WING 

15 

WTNG 

16 

WING 

17 

WING 

1R 

WING 

10 

WING 

20 

WTNG 

21 

WING 

22 

WING 

23 

WTNG 

2  A 

WTNG 

25 

WING 

26 

WI»'G 

27 

WING 

?* 

WING 

?9 

WING 

30 

w  JNG 

31 

WING 

3? 

WING 

33 

WING 

3  A 

WING 

35 

WING 

36 

WING 

37 

WING 

3R 

WING 

39 

WTNG 

AO 

WING 

*1 

WING 

A2 

WING 

A3 

WING 

AA 

WING 

A3 

WING 

A6 

SUMOUTIHC  g*TB*D: 


Purpose.  Subroutine  DATARD  reads  one  input  data  set  per  call, 
checks  data  for  input  errors,  and  passes  information  out  through  common. 


Conditions  of  Validity. 

(1)  Data  must  be  formated  in  80-column  card  images  as  described 
in  Appendix  C,  User's  Guide  for  the  Yaggy-Rogallo  Production  Software. 


Storage  Required. 


Octal  Words 

Decimal  Words 

PROG RAM t 

450 

296 

COMMON 

172 

122 

TOTAL 

642 

418 

Subprograms  Used ■ 
ABS,  WRITE3 
Call ing  Statement . 
CALL  DATARD 


Input  File  Inputs. 

INPUT  DESCRIPTION  UNITS 


HEAD1 

First  line  of  user-supplied  heading  (75  alpha¬ 
numeric  characters) 

N-D 

HEAD  2 

Second  line  of  user-supplied  heading  (75  alpha¬ 
numeric  characters) 

N-D 

UNIT 

Abbreviation  for  units  of  length  for  input 
values  (2  alphanumeric  characters) 

N-D 

FSV 

FSVANE 

Longitudinal  (fuselage)  station  of  center  of 
pressure  of  angle-of-attack  vane 

t 

FSc/4 

FSQCL 

Longitudinal  (fuselage)  station  of  quarter-chord 
line  at  midspan 

t 

SPAN 

Spar  of  the  lifting  surface 

1 

AR 

AR 

Aspect  ratio  of  the  lifting  surface,  AR  ■  b2/S 

N-D 

ac/4 

WSQCL 

Sweep  of  the  quarter-chord  line  of  the  lifting 
surface 

deg 

MW 


NMACH  Number  of  Mach  number  values  where  upwash  N-D 

angle  estimates  are  desired 

XMACH  Mach  number  values  where  upwash  angle  esti-  N-D 

mates  are  desired 

COMMON  Outputs.  Only  those  common  values  output  by  this  subroutine 
are  discussed  below. 


COMMON/ V 

’.A/HEAD1 ,  HEAD2 ,  UNIT 

OUTPUT 

DESCRIPTION 

UNITS 

HEAD1 

First  line  of  user-supplied  heading 
numeric  characters) 

(75  alpha- 

N-D 

HEAD2 

Second  line  of  user-supplied  heading 
numeric  characters) 

(75  alpha- 

N-D 

UNIT 

Abbreviation  for  units  of  length  for 
(2  alpanumeric  characters) 

input  values 

N-D 

COMMON/EEE/FSVANE,FSQCL, XV, SPAN, TAU, AR,WSQCL, NMACH, XMACH( 98) 


OUTPUT  DESCRIPTION  UNITS 


FSV 

FSVANE 

Longitudinal  (fuselage)  station  of  center  of 
pressure  of  angle-of-attack  vane 

l 

FSc/4 

FSQCL 

Longitudinal  (fuselage)  station  of  quarter-chord 
line  at  midspan 

t 

Xv 

XV 

Longitudinal  distance  from  midspan  of  quarter- 
chord  line  to  longitudinal  station  of  center  of 
pressure  of  angle-of-attack  vane,  Xv  »  Fsc/4-FSv 

t 

b 

SPAN 

Span  of  lifting  surface 

t 

T 

TAU 

Dimensionless  distance  from  quarter  chord  to 
center  of  pressure  of  angle-of-attack  vane 

N-D 

AR 

AR 

Aspect  ratio  of  the  lifting  surface,  AR  *  b2/S 

N-D 

Ac/4 

WSQCL 

Sweep  of  the  quarter-chord  line  of  the  lifting 
surface 

deg 

NMACH 

Number  of  Mach  number  values  where  upwash  angle 
estimates  are  desired 

N-D 

XMACH 

Mach  number  values  where  upwash  angle  estimates 
are  desired 

N-D 

l« 


Messages. 

(1)  "  "  IS  IN  ERRROR  -  CARD  READS  "  " 

Error  message  printed  when  DATARD  encounters  an  unexpected 
card  identifier.  First  quotation  marks  contain  card  identifier  of 
card  expected  and  second  prints  80  alphanumeric  characters  on  card 
actually  read. 

(2)  "  "  IS  IN  ERROR  -  END-OF-FILE  ENCOUNTERED 

Error  message  printed  when  an  end-of-file  in  encountered 
during  reading  of  any  card  other  than  a  HEAD1  card.  Quotation  marks 
contain  card  identifier  of  card  expected  when  end-of-file  was 
encountered. 

(3)  INPUT  DATA  ERROR: 

Error  message  printed  when  an  error  is  encountered  in  the 
input  data.  It  will  be  followed  by  one  or  more  of  the  following 
messages : 

(a)  DISTANCE  FROM  WING  TO  VANE  MUST  BE  NON-ZERO 

Error  message  printed  when  the  fuselage  station  of  the  center 
of  pressure  of  the  angle-of-attack  vane  ( FSVANE)  exactly  equals  the 
fuselage  station  of  quarter-chord  line  at  midspan  ( FSQCL) . 

(b)  SPAN  MUST  BE  NON-ZERO,  POSITIVE  VALUE 

Error  message  printed  when  the  span  of  the  lifting  surface 
(SPAN)  is  zero  or  negative. 

(c)  ASPECT  RATIO  MUST  BE  NON-ZERO,  POSITIVE  VALUE 

Error  message  printed  when  the  aspect  ratio  of  the  lifting 
surface  (AR)  is  zero  or  negative. 

(d)  SWEEP  ANGLES  LESS  THAN  0.0  ARE  INVALID 

Error  message  printed  when  the  sweep  of  the  quarter-chord 
line  of  the  lifting  surface  (WSQCL)  is  less  than  0.0. 

(e)  SWEEP  ANGLES  GREATER  THAN  90.0  ARE  INVALID 

Error  messaqe  printed  when  the  sweep  of  the  quarter-chord 
line  of  the  lifting  surface  (WSQCL)  is  greater  than  90.0. 

(f)  NUMBER  OF  MACH  NUMBERS  MUST  BE  GREATER  THAN  0 

Error  message  printed  when  the  number  of  Mach  number  values 
where  upwash  angle  estimates  are  desired  ( NMACH )  is  equal  to  0. 


Program  Listinc 


SUPPOUTTne 

common  maa 

CnMMON/Fcc 

CHARACTER 

rwARACTFO 

tat*  epri/ 
DATA  E PR?/ 
r>ATA  F°P7/ 
DATA  F  °RA 1 
PATA  FOR'/ 


DATA  fppa/ 
OEAP  IHFAP 
HFAOFP.tmf 
P  FAD ( 5* 100 
T  E  f  H  E  A  0  •  M 
PFAP  * UCAD 
MFAPFPa'ME 
READ  (5 » 1°P 
I  F  (  HE  A  P  .VI 
PFAP  'MINI 
HEAPFP«*WT 
PEAD(K,m 
I  F ( MF  A P  . M 
PFAD  'M*rM 
MFAOCp.lMA 

pn  io  ?«i. 

BFAPU.I’O 

1F(MF4P  ,N 
IMJ  .c0. 
tejmmaph  . 

CONTINUE 


PATAR 
/4FAP1 
/EEVAM 
ME  API  * 
FRP1*« 
• OTFT  A 
«S°AN 

*  AEpcr 
»<?UFFP 

•  PWFF  n 
• MHMAC 

I  '  PAP 
API  • 

,end«« 
F.  mca 
■>«  PAP 
A07  • 

.  E  N  0«  A 

E.  UFA 
E*  PAP 
MP<  I 

,  curv«  A 

F.  MCA 
V*  pap 

PMW  * 
09,7 

.Enin,  a 
c.  MCA 

II  MM  A 
IF.  1. 1 


»HF AD?, UNIT 
P»FSOCL*XV,SRAN,TAy.A 
7*»  HF  AD?*  75  #HE  A  M3*  73, 
5*ERP?*A5,FRRT*A«, FRR 
NPf  FROM  WING  TO  V  ANF 
MUST  BE  NON-7FPO,»OSI 
T  PATIO  HUS T  BE  NON-7 
ANIGFS  LEST  THAN  0.0 
ANLGES  GREATER  THAN 
p  OF  MACH  NUMBERS  MUS 

0»  HEAP*MFA01 
PER  1  GO  TP  10 


R» WS«Cl»  NMACH,XMACH(9B) 
MEAD*5»HEADFP*5,UNIT*? 


5»HEADFP*5,UNIT*? 
A*A5»E RP5*A3,EPR6*A5 
MUST  BE  N0N-7FR-0  •/ 

TIVE  valuf  •/ 

ERO»FOSITIVE  valuf  •/ 
APE  INVALID  '/ 

90.0  A PF  INVALID  •/ 


BE  GREATFR  THAN  0* / 


0)  heap, head? 
per  »  GO  TP  10 

0’  HE AO, UNTT.ES VAN F »F SOCL »SRA N, AR .WSOCL 
ncR  j  pn  TO  30 


pi  he  ad .machs, (>mach< n,i «j»( j*mi 
per  1  go  TO  30 
CM. MACHS 

♦t.)  )  Sn  to  ?o 


RFAD  'Fun  * 

r  a  p  r>  — 

.... 

mm. 

ME APFP. TEN  0 

* 

READ (*» 1P0.EMD.AP1  MFAD 

TEfMFAO  ,Nc, 

ME aotr  ) 

GP 

TP 

30 

XV. A  ns ( rsoct - 

FI  V ANF  ) 

IFJNMiru  ,r\ 

0  ? 

CO 

TO 

60 

IF(XV  , F  0 . 

0.0) 

GP 

TO 

70 

I F ( E  PA  M  ,(E, 

0.0) 

GP 

TP 

70 

TC(AP  « l r . 

0.0) 

GP 

Tn 

70 

I F ( W S CP L  .IT. 

P.O) 

GO 

TO 

70 

IFIVSOOI  ,CT, 
p  F  Til  PM 

00.0? 

GP 

TO 

70 

backspace  * 
PFAD  <?. IPO  i 
VPTTF1A.no? 
ETPO 

WPITC(A.TAP1 

STD® 

MACHS-P 
NMAPu.o 
PALI  V*  TTC  7 
VRTTF(6.1«0? 
TF (XV  ,co 

I  F ( *  PAN  .1  = 
T  F( AP  ,ic 

TF(W$ori  ,lT 
I c ( u  EOF  L  ,0* 
TF  ( M  AP  ME  ,  E<v 
RETURN 

FOOMAT  ETATC 
EORMATI A*. A  7 
FORMAT?  A*. 7* 
FORM  AT  f A*. 7V 
FORMAT? TH1 . * 
FORMAT? 1U1  .* 
FORM  AT  ? 1u  ,« 
FORMAT  ?TU  ,T 
FNO 


?  mfap.meaOI 
0?ufapfr,mfad»hfapi 

grAPicp 


WRITE  (6 » l^O  ?  EPR1 
WRTTF (6,160 I  ERR? 
WRITE (6.160)  E  R  R1 
WRITE (6,1601  pppA 
WRITE (6,160)  E  P B* 
WRITE (6.1601  f  RR6 


TCMEMTE - - - 

A’5  ) 

7», 17, «E1 0.01 
7V, I?»  7T10.0) 
,'"'.AS"  TS  TN  FPR OR 
."".AS"  IS  IN  ERROR 
.'  TNPI'T  DATA  FRRORSt  *  ) 
A  T1  P»  A  A*  ) 


-  CARD  R  FAD*  " A5. A75, • "* ) 

-  FND-rF-FUF  ENCOUNTERED*  ) 


r AT  APD  1 
DATARD  2 
DATARD  3 
DATARD  A 
OATARP  5 
PATAPO  6 
DATARD  7 
DATARD  8 
DATARD  9 
DATAPD10 
DATARD11 
— — PATAP012 
D ATARD13 
DATARD16 
TATAPD15 

- PATARD16 

DATAPD17 
PATARD10 
D ATAPD19 
— — DATARD20 
D ATAPD21 
PATARD22 
PATARD23 
DATAPD2  A 
DATARD25 
DAT  APD26 
DATARD27 
0ATARD2B 
DATAPD29 
PATAPD30 
DATARD11 
— — PATAPD32 
DATARD33 
DATAPD16 
DATAPD35 
PATARD36 
DATAPD37 
DATAPD3B 
0 ATARD39 
PATARD60 
DATARD61 
DATAR0A2 
0 ATAPPA3 

patapdaa 
PATAPDA5 
DAT  APDA6 
DATARD*? 
DATARDAM 
PATAPDA9 
DATARD50 
PATARD51 
0ATARDJ2 
PATARD6I 
DATARD5A 
DATARP55 
DATAR056 
DATAPD57 
DATAR058 
DATARD59 
DATARD60 

- PATARD61 

DATARD6? 
DATARP63 
D ATARD6 A 
DATAPD65 
PATAPD66 
PATAR067 
PATAR06* 
DATARD69 
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msmtiSH.  1 


r 


8UIW0UT IWE  WITE3: 

Purpose.  Subroutine  WRITE3  prints  the  actual  geometric  data, 
effective  geometric  data,  and  upwash-angle  estimates  with  40  points  per 
page  and  up  to  3  pages  as  required.  One  set  of  pages  is  printed  for  each 
input  data  set  of  the  input  file. 


Storage  Required. 


Calling  Statement, 
CALL  WRITE3 


PROGRAM 

COMMON 

TOTAL 


Octal  Words 
472 
1306 
2000 


Decimal  Words 
314 
710 
1024 


COMMOM  Inputs.  Only  the  COMMON  inputs  used  by  this  subroutine  are 
described  below. 

COMMON/AAA/H EADl , HEAD2 , UN IT 


FSV 

FSVANE 

Longitudinal  (fuselage)  station  of  center  of 
pressure  of  angle-of-attack  vane 

l 

fSc/4 

I 

FSQCL 

Longitudinal  ( fuselage (station  of  quarter-chord 
line  at  midspan 

l 

1  X 

Av 

XV 

Longitudinal  distance  from  midspan  of  quarter- 
chord  line  to  longitudinal  station  of  center  of 
pressure  of  angle-of-attack  vane,  Xv  *  FSC/4~FSV 

t 

b 

SPAN 

Span  of  lifting  surface 

t 

T 

TAU 

Dimensionless  distance  from  quarter  chord  to 
center  of  pressure  of  angle-of-attack  vane 

N-D 

?  : 
j  i 


INPUT 

DESCRIPTION 

UNITS 

HEAD1 

First  line  of  user-supplied  heading  (75  alpha¬ 
numeric  characters) 

N-D 

HEAD2 

Second  line  of  user-supplied  heading  (75  alpha¬ 
numeric  characters) 

N-D  J 

l 

UNIT 

Abbreviation  for  units  of  length  for  input  values 
(2  alphanumeric  characters) 

N-D 

COMMON/E EE/FSVANE, FSQCL, XV, SPAN, TAU, AR, WSQCL, NMACH , XMACH( 98) 

i 

INPUT 

DESCRIPTION 

UNITS 

i  I 


i 


AR 

AR 

Aspect  ratio  of  the  lifting  surface,  AR  =  b^/S 

N-D 

Ac/ 4 

WSQCL 

Sweep  of  the  quarter-chord  line  of  the  lifting 
surface 

deg 

NMACH 

Number  of  Mach  number  values  where  upwash  angle 
estimates  are  desired 

N-D 

XMACH 

Mach  number  values  where  upwash  angle  estimates 
are  desired 

N-D 

COMMON/FFF/BETA( 98 ) , TOB( 98 ) , WSB( 98 ) , EARCL( 98 ) , ECLR( 98 ) , ECLD( 98 ) 

B 

BETA 

Compressibility  parameter,  6  =/l-M^ 

N-D 

t/6 

TOB 

Distance  ahead  of  quarter-chord  line  corrected 
for  compressibility 

N-D 

aB 

WSB 

Sweep  of  quarter-chord  corrected  for  compres¬ 
sibility,  tan  A g  =  tan  Ac/4/B 

deg 

eAR/C L  EARCL 

Lifting  surface  upwash  parameter 

N-D 

e/CL 

ECLR 

Upwash  angle  per  unit  lift  coefficient  in  radians 

rad 

E/CL 

ECLD 

Upwash  angle  per  unit  lift  coefficient  in  degrees 

deg 

Printer  < 

Outputs. 

OUTPUT 

DESCRIPTION 

UNITS 

HEADl 

First  line  of  user-supplied  heading  (75  alpha¬ 
numeric  characters) 

N-D 

HEAD2 

Second  line  of  user-supplied  heading  (75  alpha¬ 
numeric  characters) 

N-D 

UNIT 

Abbreviation  for  units  of  length  for  input  values 
(2  alphanumeric  characters) 

N-D 

FSV 

FSVANE 

Longitudinal  (fuselage)  station  of  center  of 
pressure  of  angle-of-attack  vane 

l 

FSC/ 

4  FSQCL 

Longitudinal  (fuselage)  station  of  quarter-chord 
line  at  midspan 

l 

xv 

XV 

Longitudinal  distance  f-om  midspan  of  quarter- 
chord  line  to  longitudinal  station  of  center  of 
pressure  of  angle-of-attack  vane,  Xv  =  FSC/4~FSV 

t 

b 

SPAN 

Span  of  lifting  surface 

i 

I 

TAU 

Dimensionless  distance  from  quarter  chord  to 
center  of  pressure  of  angle-of-attack  vane 

N-D 
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AR 

AR 

Aspect  ratio  of  the  lifting  surface,  AR  =  b^/S 

N-D 

ac/4 

WSQCL 

Sweep  of  the  quarter-chord  line  of  the  lifting 
surface 

deg 

NMACH 

Number  of  Mach  number  values  where  upwash  angle 
estimates  are  desired 

N-D 

XMACH 

Mach  number  values  where  upwash  angle  estimates 
are  desired 

N-D 

8 

BETA 

Compressibility  Parameter,  8  =/l-M^ 

N-D 

t/b 

TOB 

Distance  ahead  of  quarter-chord  line  corrected 
for  compressibility 

N-D 

a8 

WSB 

Sweep  of  quarter-chord  corrected  for  compres¬ 
sibility,  tan  A g  =  tan  Ac/4/8 

deg 

eAR/CL 

EARCL 

Lifting  surface  upwash  parameter 

N-D 

e/CL 

ECLP. 

Upwash  angle  per  unit  lift  coefficient  in  radians 

rad 

e/cL 

ECLD 

Upwash  angle  pet  unit  lift  coefficient  in  degrees 

deg 

113 


V 


Program  Listing, 


uca  nl*7*,HEAD2*75,UNTT*2#TITLEl*32»TITLE2*A6 
LTMFi***, LINF2*50»LINF3*5C> LINEA*50»LINE5*50»LINF6*50 


WRITF3 
WR I TF  3 
WRITF3 
WPITF3 
WPITE3 
VRITF3 
UPITF3 


PATA 
PAT* 
HAT* 
DATA 
o  ATA 
PATA 
PATA 
PATA 
PATA 
l 

PATA 


F  (IPR  fll'TTNF  WRTTC9 
r  HHH  nN  A  A A  A /WF  AO  1 »  ME  A0?»  UN  I f 

COMMON /ccc /CFV ANC,e$OCL#XV,  SPAN,TAU»AR»  WS CCL»NMACH,XNACW(98) 
rnwMON/ccc  far  TAfOP)  ,  T03  ( 9  ft)  ,WS  8  (98  )  *E  AR  Ct  C9P)  ,  E  CL  R  f98 )  »  FC 1 0  (OR  ) 

DTMFNSTnu  LT‘|FiP'?),LINFll(2) 

r uaractco  —  -  - 

character 

PMiPifTFO  !.TNC7*«i0j  LINE**50.LINF9*33»  LINE  10*’7#  LlNEll*37, |  TNE1?*A1WRITE3 
CHARACTCP  PAFHf#6*DA$N10*10»PAEV]5*15  WRITE3 

PATA  TTT(  FT  M  YAP^>'-PPCALLO  UOWASH  CALCULATION*/ 

PATA  TT  Tl  c  2  /  •  HR V  A  SH  FSTTMATE  FOR  LIFTING  SURFACE  CnNpONeN  TS  * / 

PATA  PAFU$,pa*mi  0»OASH15/» - *,  » - »,« — - - •  / 

PATA  l  TNF 1  /  •  L  Ir  TTNG  SURFACE  GEPKFTRIC  DETAIL)*/ 

ITNF7/*LPNPTTU0TNAL  stattpn  of  vanfs  . •/ 

I TNC3/M  HNOTTUriNAL  STATTON  OF  A PF X  OF  OUA PTER-CHORO  tINF*/ 
l TNF4 /• PISTAWCE  from  APFX  OF  OUA RT€R -C HOPP  TO  VANe  S  (XV).'/ 

|TNF<S/*S»AN  (A)  . ., . . . ....*/ 

(TNC«,/*TAI)  (W/*B/2n  . . . . / 

[TNFT/tKPfAT  PATIO  . . . . . ..........*/ 

ITNFA/*fucpp  OF  CUARTFR-CHORD  LINE  . . . . ......*/ 

L I NF  9 / •  |  TFT  Twp  SUPFACF  UPWASH  ESTIMATES-'/ 
l  TNF  1  PM  FF  FF  CT TVF  •  » 

*  CRS  ILON  AP  '/ 

LTNC1  1  ' • M AC  h  RFT A  TAU/RCTA  WING  SWFE*, 

»P  /PL  EPSTLON/CL  FPSILPN/CL*/ 

DATA  l TNel ’/ * {nee  1  (RAO)  (RAOI  ( DEG ) • / 

PRTNT  GFOMCjRf r  o  AT  A  FOR  CYLINPEPICAL  COHPDNFNTS 
NPGJ  .  (mh«'-m/*p)+t 

I  PT-0 

pn  20  T.1,9 

PPTNT  DAPC  U.F4PTNCS  FOP  GEOMETRIC  OAT  A 
WRTTF(*S,10O*  TTTIFl 

.  .  TTTtF2,I,NPGS 

PASMVJ,0AFH15.  0APH15,  0A5H15.  DA  EMII 
mcapt 
hf  a  n  2 

PAPMIS.OAFHI?,  DASHI5,  0ASH15, DASH)} 

I  TNC1 

l  T*  FP»  FSVANF.ONIT 
l TWC7.fSCCL.UNIT 
I  TNC*. XV.UNTT 
1 Twc«, SPAN, UNIT 
l  T  N  C  A,  TAU 
I  TNF 7.  AP 
I  TNC*. W50CL 
0)  PE  TURN 
l  T  M  F  9 
IIWCTO 
I  T  *'  c  T  1 
lTNC-12 


V0TTE(6,T10) 
VRITF(6.1’0» 
WRTTF ( A . 1 9  0  > 
woTTFM.MO) 
WRTTF  IA.OOI 
uPTTFf A, lAO* 
VRITF(/».l«ft) 
urjtf(6,VS0) 
vrttf( p.ipoa 
VRTTEM.MO) 
WRITE(6,1 AOi 
WRTTMA.lTo' 
u°TTE(#.,T«0) 
TCfHMAPN  .cp. 
VRTTF(6,I  40* 
W'PTTCf  A.T901 
W»TT  F  (  A  ,1  90  T 


WR  T  T  c  (‘-.POP* 

VRITF(a,?10)  nASNF,PASHft.OASH10»OASWie.OASHIO»DASHiO»OASH10 

PRINT  pnc  ok  nr  nc  r-FOMFTRTC  DATA - - - 

po  10  J  «1  ,  A  9 


l  PT.|  t>T.1 
T  c ( l °T  ,PT  . 
U  R I Tc I 6,7’0* 


1 


f  PNT  T  N  l  IF 
rONTTNllc 

p  FTIIRw 

FORMAT  ^tatcmckits 


NNAPU)  RFTtJRN 

X**apu(  LpT).BFTA(LPT).  TOP  (LPT  *»  WS"(LpTT.f-APrL  HPT), 
ED  P  H  PT).FCLO  (LPT  ) 


l  0 
0 

I  3C 
1  AO 
i  *0 
160 
)  TO 
1  PO 
1  90 
?pp 
■>  i  o 

>?0 


r  [)PM  ATf  1^1 
c  non  AT  >  lu 
cnPMATMM 

CORNATflM 
FOPMATMW 
F  PPM  AT  Mu 
CPPWATflM 
format  mu 
FORM  AT  Mu 
FHRWAT  Mu 
FORMATflu 
FOR  WAT (IN 
F  OPMAT  MU 
r  NO 


TPA,  AO?,  T73»  '  PAGF' ) 

T1 7*  A  A  6, T72, I?,*  /«, I?) 

•5A1*  ) 

A  70  * 

ATM 

TTO,  AOO, 3X»  F10.A  »1X, A2) 
TTO, A 00»  3X,F10.A>1X. »«/?') 
TTO.A'O* 3X.F10.A.1X,  •  *) 

T10.AOO.3X, CIO, A, IX. *  OF  0*) 

n,A’7,  A37) 

T  P  2,  A*?) 

T9.A6 
1 V.2 ( c^ 


.  ’X. A6,  2X,A 10»  ?X, A 10» ?X,A10,2X, A10, ?X» A10  > 
cp.  A,  ?X).?(  F10.A,?X),3(  F10.F-.2X  n 
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WRTTF310 
VRITF311 
WPITF3I2 
WRITF313 
WRITF31A 
WPITF315 
WRITE  316 
VRTTF317 
WPITF31B 
WRI TF  319 
WRTTF 320 
WRITF  321 
WPTTF322 
WRTTF  323 
WPITF32A 
WRITF  325 
WRITE326 
VPTTF327 
WRITF32B 
WRITC329 
WRITF 330 
VRITE331 
WRTTE332 
WRITF  333 
WPTTE33A 
WR I TE  335 
WRITF  336 
WRITF337 
WRITF33B 
WRTTE339 
WRITF3A0 
WRTTF3A1 
WRTTE3A2 
WRITF3A3 
WRITF3AA 
WRITE3A5 
WPTTF3A6 
WRITF3A7 
WRITE  TAB 
WRITF3A9 
WRITF  350 
WR JTE351 
WRITF352 
WR I TF  35  3 
WPITF35A 
WRTTE355 
WRTTF  356 
WRTTF  3J7 
WRITF3^A 
WRTTF  359 
WRITF 360 
•WR  I  TF  361 
WRITF362 
WPTTC363 
WRITE  36 A 
WRITF 363 
WRTTF 366 
WRITF367 
WRITF76F 
WRITF  369 
WPTTF370 
WR I TF  371 
WRTTF372 
WR I TF  373 
WRTTC37A 
WRTTF  375 


